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We theoretically investigate the backaction of a sensor quantum dot with strong local Coulomb 
repulsion on the transient dynamics of a qubit that is probed capacitively. We show that the mea¬ 
surement backaction induced by the noise of electron cotunneling through the sensor is surprisingly 
mitigated by the recently identified coherent backaction [M. Hell, M. R. Wegewijs, and D. P. Di- 
Vincenzo, Phys. Rev. B 89, 195405 (2014)] arising from quantum fluctuations. This indicates that 
a sensor with quantized states may be switched off better than naively expected. This renormaliza¬ 
tion effect is missing in semiclassical stochastic fluctuator models and typically also in Born-Markov 
approaches, which try to avoid the calculation of the nonstationary, nonequilibrium state of the 
qubit plus sensor. Technically, we integrate out the current-carrying electrodes to obtain kinetic 
equations for the joint, nonequilibrium detector-qubit dynamics. We show that the sensor current 
response, level renormalization, cotunneling broadening, and leading non-Markovian corrections al¬ 
ways appear together and cannot be turned off individually in an experiment or ignored theoretically. 
We analyze the backaction on the reduced qubit state - capturing the full non-Markovian effects 
imposed by the sensor quantum dot on the qubit - by applying a Liouville-space decomposition 
into quasistationary and rapidly decaying modes. Importantly, the sensor cannot be eliminated 
completely even in the simplest high-temperature, weak-measurement limit since the qubit state 
experiences an initial slip depending on the initial preparation of qubit plus sensor quantum dot. 
The slip persists over many qubit cycles, i.e., also on the time scale of the qubit decoherence in¬ 
duced by the backaction. A quantum- dot sensor can thus not be modeled as usual as a “black box” 
without accounting for its dynamical variables; it is part of the quantum circuit. We furthermore 
find that the Bloch vector relaxes (rate l/Ti) along an axis that is not orthogonal to the plane 
in which the Bloch vector dephases (rate I/T 2 ), blurring the notions of relaxation and dephasing 
times. Moreover, the precessional motion of the Bloch vector is distorted into an ellipse in the tilted 
dephasing plane. 

PACS numbers: 73.63.Kv, 73.63.-b, 03.65.Yz 


I. INTRODUCTION 

The ongoing effort to mitigate the qubit decoherence 
due to environmental noise detrimental to quantum com¬ 
puting has made substantial progress by identifying well- 
isolated two-level systems [lHj] and developing efficient 
decoupling techniques However, active readout el¬ 

ements must be integrated into any quantum computer 
and “noise” from such sensors may soon become a rel¬ 
evant source of errors. Therefore, the unavoidable dis¬ 
turbance of the qubit evolution during a readout process 
gains importance, both for single-shot qubit readout 
as well as continuous qubit monitoring 0- Another 
question of central importance for qubit manipulations 
is which parameters should be varied to switch off a sen¬ 
sor most effectively when no measurement is intended to 
be made. 

In a single-shot measurement, the goal is to achieve a 
strong, projective measurement that dephases the qubit 
state as quickly as possible. However, any measurement 
still takes a finite time and relaxation processes @, exci¬ 
tation processes [Tlj] . and incoherent detector dynamics 
limit the detector efficiency p~2l4l5l | . For continuously 


monitoring the qubit evolution, by contrast, one has to 
realize a weak measurement. The aim is here to disturb 
the qubit evolution as weakly as possible to retain the 
(partial) purity of the quantum state [0, while avoid¬ 
ing quantum jumps and the quantum Zeno effect 01- 
Understanding the backaction of a weak measurement is 
therefore of great interest. 

In this paper, we focus on the noninvasive, weak- 
measurement backaction exerted on a qubit b y a capaci¬ 
tively coupled sensor quantum dot (SQD) @,[0. SQDs 
are attractive qubit detectors due to their strong tun- 
ability and higher signal-to-noise ratio as compared to 
quantum point contacts 0, 0] and dispersive readout 
schemes [0. This derives from the fact that an SQD is 
an interacting quantum system. For smaller quantum 
dots (QDs) than typically used for readout, the elec¬ 
trons may even occupy discrete energy levels. The re¬ 
quired readout current easily leads to a strong nonequi¬ 
librium and nonstationary sensor state. This altogether 
makes the description of the backaction arising from an 
SQD on a qubit challenging. Thus, typical weak-coupling 
approaches to decoherence assuming an environment in 
equilibrium with a continuous spectrum 10 - 10 cannot 
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be applied here and naive extensions are prone to er¬ 
rors as we will illustrate. Understanding the intertwined 
evolution of SQD and qubit is of key importance to un¬ 
derstand the measurement backaction jlj, [20l U2l . 

Out of these challenges arises the question which phys¬ 
ical effects have to be included for a consistent descrip¬ 
tion of the measurement backaction of a sensor QD on 
a qubit. For the qubit, we consider the simplest model 
that involves capacitive readout, a charge qubit, realized 
as a double quantum dot. An important aspect lies in 
the type of setup considered, namely that of indirect de¬ 
tection: one measures the conductance of a sensor QD in 
the attached electrodes, while only the SQD capacitively 
interacts with the qubit, see Fig. [T] To maximize the 
sensitivity, the SQD is operated at the threshold to the 
Coulomb blockade regime. Here, the conductance of the 
SQD shows the strongest response to small qubit-induced 
level shifts. Previous studies addressing the backaction 
of an SQD on a qubit in this regime focused mainly on 
the lowest-order approximation in the tunnel coupling T 
of the SQD to the attached electrodes [HIIEI2 §mI 3> 
strictly valid when operated in the single-electron tunnel¬ 
ing (SET) regime. There are two main reasons for going 
beyond this approximation. 

The first reason relates to the backaction on the qubit 
and its dependence on the level position, experimentally 
controlled by gate voltages. This is important since the 
level position is one of the key experimental control pa¬ 
rameters by which one can try to switch off the sen¬ 
sor backaction. In this approximation, the leading-order 
rates (SET oc T) become exponentially small when the 
level position (e) of the SQD is tuned away from the 
electrochemical potentials of the electrodes (/ i r ). Thus, 
when one is interested in the backaction at the onset of 
Coulomb blockade, next-to-leading order oc T 2 cotunnel¬ 
ing processes should also be accounted for because they 
are only algebraically suppressed, scaling as l/(e — p r ). 
One would thus naively expect that the backaction is 
suppressed only inversely proportional to the detuning 
from resonance. Yet, level-renormalization effects should 
also be considered; they lead to level shifts that depend 
logarithmically on the level position e and therefore the 
response of the level renormalization to the measurement 
perturbation also scales as l/(e — p r ). A central finding 
of our study is that level-renormalization effects in fact 
mitigate the naively expected cotunneling decoherence. 

A second reason for going beyond the lowest-order ap¬ 
proximation comes in view when one takes into account 
the experimentally used sensor signal: if one accounts for 
the first nonvanishing contributions oc XT/T that pro¬ 
duce a nonzero sensor signal, one has to include also 
renormalization effects since they appear in the same or¬ 
der. Basically, if one has time to measure, one has time to 
fluctuate as well. As already emphasized in Ref. [32} , in¬ 
corporating terms oc AF jT is another reason that forces 
us to keep also cotunneling processes oc T 2 /T since in 
the weak-measurement limit ACT the latter are larger. 
Only when one is not interested in the sensor current, 


one can consistently neglect cotunneling and renormal¬ 
ization effects by taking the high-temperature limit: as 
we show, they must either be kept or neglected together. 

The above-mentioned processes combine in a nontriv¬ 
ial way to give three types of backaction on the qubit 
[32} . First, both SET and cotunneling processes con¬ 
tribute to a stochastic switching of the SQD charge state. 
This switching generates a randomly fluctuating effective 
ma gnet ic field acting on the qubit Bloch or isospin vector 
|26l |27| . This “noise” - called here the stochastic back- 
action - leads to a shrinking of the Bloch vector, i.e., 
to decoherence. In addition, there is a dissipative back- 
action , which is the flip-side of the measurement action: 
it arises whenever one accounts for a nonzero response 
of the sensor SET tunnel rates to the qubit state and 
therefore a nonzero sensor signal. Finally, there is a co¬ 
herent backaction , the most striking finding of Ref. {32j. 
It arises from the above-mentioned level-renormalization 
response and translates into torque terms involving the 
qubit Bloch vector. These torques and related precession 
effects are similar to those emerging in various other QD 
transport setups: It is well-known that tunneling pro¬ 
cesses can produce exchange fields leadin g to an (iso)spin 
precession in the context of spintronics J34l [35 ] , double 
dots [Hi], molecular quantum dots H3, HH, and super¬ 
conducting devices [391 . All these level-renormalization 
effects arise from quantum fluctuations of electrons by 
tunneling into the attached electrodes. In this respect, a 
qubit coupled to a sensor QD is not different. 

An interesting question is how these different types 
of backaction relate to the information gained during 
the measurement process. Clearly, the quasistationary 
time-dependent current through the sensor contains in¬ 
formation about the qubit state and at the same time 
causes decoherence of the qubit. However, besides this 
fundamentally unavoidable backaction, the decoherence 
induced by the detector can be stronger. This can be for¬ 
mulated in terms of general inequalities relating the noise 
of the measured operator (here the position of the qubit 
electron) to the noise of the measurement signal (here 
the current) and additional noise cross terms [77}. The 
situation considered in this paper is far away from the 
quantum limit (meaning the above-mentioned inequality 
is far from being satisfied with equality). 

A part of our work actually focuses on a simple limit 
when only the stochastic backaction is accounted for 
but dissipative and coherent backaction are neglected. 
The SQD then acts rather as an “ordinary” environment 
and no information is obtained during the operation - a 
situation relevant when the detector is supposed to be 
switched off. Yet, even in this simple limit, there are ef¬ 
fects beyond the scope of the picture developed in Ref. [77|: 
the fast relaxation due to switching on the sensor may 
affect the qubit state also in a way that depends on the 
initial dynamical variables of the sensor. This is not cap¬ 
tured at all by the noise inequalities mentioned above. 
The explicit time evolution of the sensor, which is not 
considered in Ref. IB is thus also crucial to understand 
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the backaction on the qubit. An interesting question aris¬ 
ing from this insight is how this backaction effect has to 
be assessed in view of the information gain. Our work 
could thus spur new activity on the topic of information 
gain during the measurement. 

The impact of the three backaction effects on the full 
transient dynamics of the qubit so far remained an out¬ 
standing question that we address in this article (the 
analysis in Ref. ![32j was restricted to the stationary 
state). Our analysis is divided into two parts. 

Part 1. The main point of this paper is to eliminate 
the electrodes’ degrees of freedom [33 and to analyze 
the transient dynamics within the resulting physical pic¬ 
ture of the coupled SQD-qubit dynamics. In this way, 
we can deal with the nontrivial interplay of the SQD- 
qubit coherence, strong local Coulomb interaction in the 
SQD, nonequilibrium conditions imposed by the attached 
electrodes, as well as both leading (SET) and next-to- 
leading order effects in the tunneling (cotunneling). The 
necessary inclusion of the latter furthermore forces us 
to extend Ref. [32j] by including also the leading mem¬ 
ory effects on the sensor-qubit system due to the tunnel¬ 
ing to the electrodes, which is necessary for the study of 
the transient qubit dynamics. (For the stationary state, 
which we studied in Ref. [32j, they can be be ignored 
without making additional approximations.) The impor¬ 
tance of memory effects for the dynamics when g oing 
beyond weak coupling^ is known especially since Ref. [78], 
see also Refs. Jdol . Ell f79l - l82[ and progress for strong 
backaction \EHm . Non-Markovian corrections have also 
been studied in related contexts, such as the backaction 
of a quantum point contact on a double dot [gl], HI] or 
quantum-feedback control based on quantum measure¬ 
ments [HI, El]. The various effects of non-Markovian 
processes remain, however, an uncharted territory (87). 
The central result in our case are the kinetic equations 
(126)) for the system of qubit plus quantum-dot sensor. 
The equations reveal the above three-fold nature of the 
backaction of the sensor QD on the qubit; importantly, 
the relevant energy scale for this backaction is not sim¬ 
ply the internal capacitive interaction A (SET-induced 
stochastic backaction) but additionally involves the en¬ 
ergy scale TX/T (dissipative and coherent backaction in¬ 
volving transport processes). 

Our kinetic equation furthermore allows us to iden¬ 
tify slowly evolving quasistationary modes - containing 
the qubit evolution - and faster evolving decay modes re¬ 
flecting the dissipative SQD dynamics due to its coupling 
to the electrodes. The coupling between these modes 
generates the total backaction on the qubit and is medi¬ 
ated by all three types of backaction. To account for all 
these backaction effects, it is indispensable to keep the 
capacitive interaction (A ^ 0) when integrating out the 
electrodes. In this aspect, our work differs from the oth¬ 
erwise closely related approach of Ref. [88], which starts 
out from the assumption that the electrodes affect exclu¬ 
sively the SQD. There, all backaction effects derive only 
from the internal interaction, i.e., the stochastic backac¬ 


tion. 

A surprising finding of our analysis is that the total 
backaction exhibits a strong reduction when tuning the 
SQD towards the Coulomb blockade regime: we find that 
the coherent backaction actually cancels the cotunneling 
(“broadening”) corrections in the coupling of the quasis¬ 
tationary to the decay modes. This eliminates the naively 
expected leading power-law dependence oc l/{e — / i r ) of 
the backaction, affecting also the decoherence time scales. 
This indicates that a sensor with quantized orbital states 
can be switched off more efficiently by controlling its gate 
voltage than naively expected, the first important experi¬ 
mental implication of this article. This requires, however, 
to prepare the sensor state in a controlled way to avoid 
a slip of the qubit state (see part 2 below). 

It is important to emphasize already here that the 
coherent backaction, which is responsible for this mit¬ 
igation, is not an independent mechanism that can be 
“added” to counteract cotunneling noise. Instead, it 
arises together with cotunncling as an integral part of 
quantum fluctuation effects of the qubit-sensor system 
when consistently describing all types of backaction. No¬ 
tably, we show that this mitigation is not captured by 
widely-used classical stochastic fluctuator models and 
can also be easily overlooked in Born-Markov approaches 
that integrate out the entire environment of the qubit 
(i.e., including the SQD). The effect of the exchange of 
electrons between the SQD and the electrodes can thus 
not be fully captured by classical switching of the SQD 
charge state. We also review and compare in detail our 
results with earlier works and pinpoint a number of lim¬ 
itations of standard approaches. 

The prominent role of renormalization effects under¬ 
lying the coherent backaction distinguishes a QD sensor 
with few, discrete energy levels from what can be ex¬ 
pected for a sensor with a continuous energy spectrum. 
The recent study [89] showed that similar torque terms 
appearing in a spintronic context are much suppressed in 
single-electron transistors (continuous spectrum) as op¬ 
posed to QDs (discrete spectrum). In the former case, 
renormalization effects tend to nullify when averaging 
over their continuous energy spectrum. This motivates 
the extensive analysis of the detection of a qubit state by 
a sensor QD undertaken in this paper. Our work raises 
the interesting question to which extent backaction ef¬ 
fects due to renormalization effects are suppressed in a 
single-electron transistor. 

Part 2. One might think that following the above de¬ 
scription one can in a second step eliminate the sensor 
QD from the description to obtain an effective theory for 
the qubit only. However, already on general grounds, 
this is questionable: specific to our indirect detection 
problem is that the environment of the qubit is not sta¬ 
tionary. Moreover, initial correlations between SQD and 
qubit - both microscopic systems - might exist. When 
integrating out the environment, the factorizability and 
the stationarity of the environment are, however, often 
invoked to eliminate the so-called slip of the initial condi- 
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tion for the subsystem [90] coming from the initial state of 
the environment and its short-time transient evolution. 
To illustrate this point, we analyze in more detail the 
simpler high-temperature limit where the complications 
due to cotunneling and coherent backaction can be con¬ 
sistently ignored. Even in this high-temperature, weak- 
measurement limit the qubit develops a slip of order A/T 
on a time scale t < 1/T, which is beyond the control 
over the qubit system alone. The slip depends explicitly 
on the initial qubit-sensor state. The second experimen¬ 
tal implication of our work is that the dynamical state 
of the sensor and its correlations with the qubit cannot 
be ignored and must be brought under experimental con¬ 
trol. This slip effect is cumulative, e.g., it results in phase 
shifts that still affect the qubit on much longer time scales 
1/A^> 1/r relevant to the readout. By contrast, the 
relevant time scales Tj (relaxation) and T 2 (dephasing) 
of the transient qubit dynamics for times t 1/T do not 
depend on the initial sensor QD state. 

The precession axis of the qubit Bloch vector also turns 
out to be independent of the initial state. In the simple 
high-temperature limit, we furthermore identify an addi¬ 
tional effect of the (purely stochastic) backaction, which 
is to induce a tilt of the Bloch vector precession axis. We 
find that the circular isospin precession becomes slightly 
elliptical in the presence of the detector, adding as a 
fingerprint oscillations to the exponential decay of the 
qubit-state purity. This mixes the notions of relaxation 
and dephasing as we will see. It is an interesting question 
how these effects behave at low temperature and strong 
qubit-sensor coupling where they have received little at¬ 
tention so far. 

Outline. After this topical outline, we now present 
the organization of the sections of the paper and the key 
equations. In Sec. El we briefly review the generic indi¬ 
rect readout model of Ref. [32[ and discuss the dynami¬ 
cal variables that are needed to describe the mixed quan¬ 
tum state of the joint qubit-sensor system. This requires 
two isospins, which capture both the reduced qubit state 
as well as the correlations with the sensor QD. 

After this, we outline the key technical challenges of our 
approach in Sec. IIIII deferring details to the App. [A] and 
we present the time-local kinetic equation (l26l) for the 
coupled sensor QD plus qubit system. In App. [X] we fur¬ 
ther discuss the importance of including non-Markovian 
corrections to retain the positivity of the reduced den¬ 
sity operator. Without further approximations, we iden¬ 
tify the relevant unperturbed modes (A = 0) with the 
electrodes integrated out. From the representation of 
the kinetic equation (1561) in these modes we prove the 
exact cancellation between the coherent backaction and 
the cotunneling (“broadening”) noise [Eq. (1521) ]. We fur¬ 
thermore study the implications for the dependence of 
the total backaction on various experimentally relevant 
parameters (tunneling-rate asymmetries, bias, and gate 
voltages). From the formal solution of the effective qua¬ 
sistationary mode evolution [Eqs. (1711) and (1771) ]. with 
details given in Appendix El we infer that the qubit evo¬ 


lution is non-Markovian and exhibits a slip of the ini¬ 
tial condition that we characterize in Appendix EJ Ini¬ 
tial slips generally go h and in hand with non-Markovian 
dynamics as Refs. El l90l - l94| and the references therein 
point out. Initial entanglement between a qubit and its 
environment, one cause of initial slips, can drastically af¬ 
fect the qubit coherence [94| . 

In Sec. hyi we attempt to integrate out the sensor QD to 
derive an effective Liouvillian [Eq. (l80l) ] that effectively 
incorporates its fast switching dynamics. We then focus 
on the analytically tractable case of high temperature. 
This suffices to illustrate the general importance of initial 
slips in the context of detector backaction [Eq. (ISll) ]. the 
breakdown of orthogonality of relaxation and dephasing 
qubit modes, and the exponentially damped oscillatory 
but elliptical precession of the qubit Bloch vector. In this 
high-temperature limit, we obtain tangible expressions 
for the qubit relaxation and dephasing rates, expanded 
to leading order in A/T. In Sec. El we compare our results 
with semiclassical stochastic fluctuator models as well as 
Born-Markov and exact quantum approaches to provide 
further insight into the origin of the coherent backaction. 
In the accompanying Appendix El we show how the co¬ 
herent backaction affects the qubit phase evolution in a 
way that is not accounted for by semiclassical stochas¬ 
tic fluctuator approaches. We summarize our findings in 
Sec. ED 



FIG. 1: Sensor quantum dot (SQD) tunnel-coupled to source 
and drain electrodes and capacitively coupled to a qubit, 
whose different logical states involve two possible positions, 
left and right, in a double quantum dot. If the qubit electron 
is left (a) or right (b), the Coulomb repulsion to the SQD 
electron is larger or smaller, respectively, compared to the 
full delocalization of the qubit electron. 
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II. INDIRECT DETECTION 
A. Model 

We analyze the indirect detection setup sketched in 
Fig. [TJ the readout (Hi) of a double-quantum dot charge 
qubit ( Hq) by a proximal sensor quantum dot ( Hs ), 
which in turn is read out (Hr) by the conductance of 
the transport current in one of the attached electrode 
reservoirs (Hr). The model we use, discussed in detail 
in Ref. [32|, thus consists of three “layers” with their 
respective interactions: 

H = Hq + (H i + H s ) + (H t + H r ). (1) 

This models the essential physics found in many exper¬ 
iments on QD qubits and can be extended to supercon¬ 
ducting qubits [45[ as well as to spin qubits if measured 
by spin-to-charge conversion [3, [l!| ■ 

Qubit. The qubit is realized as a charge qubit, a single 
electron occupying a double quantum dot. This electron 
can reside either on the left dot, denoted by the state |L), 
or on the right dot, denoted by the state | R). The qubit 
state is represented by the ensemble average r = (r) 
of an operator f (corresponding to an isospin t/ 2) with 
components 

n = E (^vio<n- (2) 

l,l'=L,R 

Here, eq denotes the Pauli matrix for i = x,y, z. The av¬ 
erage z component r 2 quantifies the imbalance between 
the probabilities for finding the qubit electron in the left 
orbital rather than in the right orbital, while t x and t v 
quantify coherences between the left and right occupa¬ 
tion. The general form of the Hamiltonian of the isolated 
qubit is 

H q = n-f/2, (3) 

in which the qubit field Cl is used in applications to con¬ 
trol the qubit evolution: it induces coherent tunneling of 
the qubit electron between the two dots ( Cl x , Cl y ) com¬ 
bined with a detuning (Cl z ). In our analysis, fl is con¬ 
stant in time and later on, when we discuss tangible re¬ 
sults, we will chose Cl = Cle x . However, unless stated 
otherwise, we first keep Cl general. 

Sensor quantum dot. The sensor quantum dot (SQD 
or sensor QD) is modeled as a single, interacting, spin- 
degenerate orbital level with Hamiltonian Hs = efi + 
I/rhj-nj,. Here, h a = is the number operator for 

electrons with spin er =j\ J,, where d a denotes the cor¬ 
responding field operator, and n = + nj. is the total 

electron number operator. We take the Coulomb repul¬ 
sion energy U here to be the largest energy scale (except 
for the bandwidth 2 W of the electrodes), in accordance 
with typical experimental situations. We therefore ex¬ 
clude the double occupation of the sensor QD orbital in 
the following. This allows us to reduce the model to 

Hs = eP\ P 1 = E>>H (4) 


if we accordingly adjust the high-energy cutoffs in the 
electrodes [discussed below Eq. (13TT) ]. In the considered 
subspace, we can replace fi = P 1 . For the readout one 
tunes the level position e = —V g by a gate voltage V g 
close to the electrochemical potentials of one of the elec¬ 
trodes. While the spin of the qubit electron is irrelevant 
(it was not written above since the readout couples to the 
charge, see below), it is important to include the spin de¬ 
gree of freedom of the SQD because the spin degeneracy 
enters into the tunneling rates. 

Electrodes. The final stage of the readout involves the 
electrodes, treated as noninteracting reservoirs of elec¬ 
trons with spin: 

Hr = ’Yh r k< y Wrkac\.k a C r kai ( 5 ) 

with the field operators c r ka ( c tka) acting on the elec¬ 
trons in orbital k with spin a in the source (r = s) and 
the drain (r = d ), respectively. These are each held in 
equilibrium with a common temperature T, but at differ¬ 
ent electrochemical potentials y s =14/2 and yd = — 14/2 
by a applying a bias voltage 14- 

Readout. The indirect readout of the qubit state using 
the sensor QD involves two couplings: the first one is the 
capacitive interaction of the SQD electron charge h with 
the charge polarization t z of the qubit: 

Hi = hX-f/2. ( 6 ) 

The measurement vector A specifies both the basis in 
which one measures and the measurement strength A: 

A = Ae z . (7) 

Thus, depending on the qubit state, the sensor QD 
level experiences an energy/gate-voltage offset of at most 
±A/2 [see Figs. [T](a) and [11(b)] . This in turn affects the 
conductance measured in one of the electrodes due to 
tunneling to and from the SQD: 

H T = J2r,k,a trdlcrka + H.C. (8) 

The strength of this second coupling involved in the read¬ 

out is quantified by the tunnel rates = 27r|f T .| 2 ^ T ., 
where we take both the tunneling amplitude t r and the 
density of states v r to be spin (a)- and energy (k)- 
independent within the electrode bandwidth 2 W. 

The threefold layered structure of this indirect detec¬ 
tion (negligible direct coupling of the electrodes to the 
qubit) is reflected in our theoretical analysis of the mea¬ 
surement backaction. In Sec. IIII Cl we first eliminate the 
“outer” detection layer - the electrodes - in favor of effec¬ 
tive equations describing the joint SQD-qubit dynamics 
[99j . In Sec. hyi we then attempt to integrate out the 
“inner” detection layer - the SQD to find the effective 
qubit evolution. 

B. Weak-measurement and weak-tunneling limit 

We consider a sensor with a fast response, i.e., the in¬ 
ternal dynamics of qubit plus sensor QD is slow as com- 



6 


pared to the electron tunneling dynamics induced by the 
attached electrodes: 


A ~ \,fl < r. (9) 


This condition means physically that many electrons pass 
through the SQD during its interaction time with the 
qubit (A -C T) — a weak measurement is performed. 
Moreover, if fi lies in the x-y plane, the internal qubit 
evolution describes a coherent tunneling of the qubit elec¬ 
tron with dwell times of electrons in the SQD that are 
much smaller than the period of a qubit cycle (f2 <C T). 
Each electron sees a “snapshot” of the SQD-plus-qubit 
state. If we assumed T < Q instead (but still weak mea¬ 
surement A -C T), the readout would be too slow to re¬ 
solve any qubit evolution. 

As we see below in Sec. uni the leading-order response 
of the tunnel rates to a measurement-induced gate volt¬ 
age offset ~ A is given by TX/T. This, in return, induces 
a dissipative backaction affecting the polarization of the 
qubit that also scales as TX/T. Moreover, when condi¬ 
tion © holds, the tunneling also affects the isospin co¬ 
herences of the qubit-SQD state. This is well-known from 
the analysis of two-level systems coupled to a reservoir: 
Here, the density-operator coherences in the energy basis 
matter when the levels (here split by A) are degenerate 
on the scale of the coupling (here T) to the environment. 
In our case, the two levels correspond to the sensor QD 
and the qubit, each being a two-state system. One has 
to carefully identify which coherences are relevant, which 
is done below in Sec. Ill Cl These coherences are affected 
by tunneling processes: The simple physical intuition be¬ 
hind this is that if an electron on the sensor has time to 
interact with the qubit and change the current, it cer¬ 
tainly has time to fluctuate into the electrodes. This 
leads to a response of the level renormalization and re¬ 
sults in a coherent backaction which scales as TX/T , i.e., 
in the same way as the response of the sensor tunneling 
rates resulting in the dissipative backaction. This is a 
central result of Ref. [32] and here we explore its effect 
on transient dynamics. 

It should thus be noted that the energy scale for back- 
action on the qubit is not simply A (from the internal 
interaction Hj) but also AT/T, the scale of effective dis¬ 
sipative and coherent coupling between sensor QD and 
qubit, which are induced by tunneling processes. In a 
way, these couplings account for an indirect interaction 
of the qubit with the electrodes extending the approach of 
Ref. |88| • This is thus another relevant perturbative scale 
for a weak-measurement expansion, besides the scale A 
itself, see Ref. 32| for a detailed exposition. 

Another crucial point for this work is that if terms 
of order TX/T are taken into account and A < T, then 
at least cotunneling terms scaling as T 2 /T must also be 
accounted for (if not even higher-order tunneling terms). 
We can neglect higher-order tunneling processes beyond 
cotunneling if we restrict the temperatures such that 


T 2 /T < A. 


This condition means that the cotunneling-induced noise 
imposes only a weak perturbation of the qubit. Taken 
together, we employ here a weak-coupling limit in two 
ways, namely that of weak measurement and weak tun¬ 
neling: 


T/T < A/T < 1 . (11) 

Note that by Eq. m this imposes a stronger condi¬ 
tion than the usual weak-tunneling assumption T/T -C 1 
alone. We next discuss the dynamical variables needed 
to describe the measurement backaction. 


C. Charge-specific isospins and qubit decoherence 


To describe both the backaction of the sensor on the 
qubit as well as to compute the signal current through the 
sensor QD, one needs at least the reduced density opera¬ 
tor p(t) of the combined qubit plus SQD system obtained 
by tracing over the electrodes. Even though we do not 
analyze the sensor signal here, the signal is of course of 
high interest to gain insight into, e.g., the efficiency of 
the measurement [© - [l5j . Studying the backaction in a 
situation where the sensor signal current is not negligi¬ 
ble, is an experimentally highly relevant situation, which 
we pursue in this paper. 

The relevant part of the SQD-qubit density operator p{t) 
can be expanded as follows [321 ]: 


Pit) 


lY. pn ® 


i 


( 12 ) 


Here, P n denotes the projector onto the charge states n = 
0,1 of the sensor QD. The numbers p n (t) = ti(P n p(t)) 
give the probability for the sensor QD to be in the respec¬ 
tive charge state n = 0,1, which for any time t sum up 
to one due to the probability conservation: p° + p 1 = 1. 
The only irrelevant coherences (off-diagonal matrix ele¬ 
ments in the energy basis) of p are those involving differ¬ 
ent charge states on the sensor. These can be shown 
to decouple from the relevant part due to the charge 
conservation by the tunneling. However, all remaining 
qubit-SQD density matrix elements including their co¬ 
herences must be kept in (fl2l) . These are the six num¬ 
bers r"(t) = ti{P n Tip(t)), which are the averages of the 
isospin components i = x,y, and z for the two sensor QD 
charge states n = 0 or 1, respectively. To describe the 
correlated SQD-qubit system, one thus needs two charge- 
specific isospins t° and r 1 . Based on Eq. m, it is con¬ 
venient to introduce the following column representation 
of the density operator: 


P = 


/ P° 



(13) 


( 10 ) 
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In this form, one distinguishes the charge and isospin 
part, respectively, but the isospin part is still kept ba¬ 
sis independent. In other words, we may represent r° 
and x 1 in a different orbital basis than the one used in 
definition © if we transform the directions of A and ft 
accordingly. 

We now further explore the physical meaning and im¬ 
portance of the two charge-specific isospins. By construc¬ 
tion, they sum up to the total isospin, 

T = t q + t 1 , (14) 

which is often of main interest. This is the usual Bloch 
vector that describes the state of the qubit, i.e., its re¬ 
duced density operator. One can easily show that a single 
Bloch vector can describe the joint SQD-qubit state only 
if it is factorizable. In other words, if 

p = ps® pq= (^2p n P n ^j ® \ ( i + > 

(15) 

then the charge-specific isospins are given by r° = p°T 
and t 1 = p 1 r, respectively, by comparing Eq. IT5l) with 
Eq. (HU). The other combination of the two isospins, 

S := p°T l — pV 0 , (16) 

quantifies the nonfactorizability of the qubit-sensor den¬ 
sity operator: 

P = PS®PQ «=> $ = 0. (17) 

We emphasize that the nonfactorizability is crucial to 
describe the readout and its backaction on the qubit state 
t . We will see in Sec. IIII C 11 below Eq. (l54ll that the 
deviation d has no impact on the qubit evolution only if 
the qubit and sensor are strictly decoupled (A = 0). For 
the coupled case (A ^ 0), which is of interest, we have to 
keep the individual dynamics of r° and r 1 . 

The relevance of the two isospins for the decoherence can 
be seen explicitly from the equation of motion for |r(t)| 2 , 
which characterizes the purity of the isospin state: 

A[|t| 2 ] = 2r-r = —2 A-(t°xt 1 ), (18) 

where we inserted after the first equality 

t = SIxt + Axt 1 . (19) 

Equation m is an exact result which can be obtained 
from the Heisenberg equation of motion for r with re¬ 
spect to the full Hamiltonian ©. The only essential as¬ 
sumption on which Eq. (1181) relies is the indirect readout 
structure of our setup, i.e., [Ht,Hq\ = 0. It thus holds 
generally for any $7 and any tunneling T. Preservation 
of this exact equation imposes an exact isospin sum rule 
[discussed below Eq. (1261) ]. which any kinetic equation 
for p should satisfy: This was only realized recently, see 


Refs. and (32j, in particular, see Appendices D and 

E. 

Equation (fl8l) shows that the reduction of the purity of 
the qubit state appears only due to noncollinearities of r° 
and t 1 . These noncollinearities develop because of the 
readout: the isospins r° and t 1 are subject to different 
effective “magnetic” fields depending on the charge state 
n of the sensor QD. To see this, we rewrite Hq + Hj = 
B c ff(n) • t/ 2 with an effective field acting on the isospin 
T, 

B e ff = $7 + X5n, (20) 

where Sn = n — (n) and (h) = tr (P 1 p) = p 1 . Here, the 
first part is the mean field ,, 

17 = 17 + (n)A, (21) 

which the isospin experiences due to the internal isospin 
field 17 and the average field caused by the mean charge 
(n) = p 1 on the sensor QD with respect to the exact 
total density operator p. The mean-field contribution 
17 to B c ff is the same for both charge-specific isospins, 
t° and r 1 , and therefore not responsible for the qubit- 
state decay. The mean SQD occupation (n) = p 1 merely 
tilts the qubit precession axis and changes its frequency 
contributing to the detuning of the qubit. Note that 
the average is here an ensemble average but not a time 
average since p 1 = p : (t) can change in time with the 
state p(t). The qubit decay is induced by the second, 
fluctuating contribution \5h to Eq. (1201) . where Sn = 
h — (n) is the charge-state dependent deviation from the 
mean field. This generates a noncollinearity of t° and r 1 , 
which reduces the purity of the qubit state by Eq. (US- 

Even though our approach does not make use of the 
decomposition into a mean-field and a fluctuating part, 
we can identify both effects in our results in Sec. IIII El We 
will first identify p 1 = pl t with the state of the SQD in 
the stationary limit, in which the ensemble and the time 
average are equal. We further connect more precisely the 
decoherence rates to the components of the fluctuating 
part A Sh along and perpendicular to the mean field 17 
in Sec. IIV Al ( in accordance with the literature mmy 
This accounts for what we call stochastic backaction on 
the qubit by the sensor QD. This effect is also present 
for single-electron transistor sensors with a continuum of 
electronic levels, but (classically) quantized charge states. 

However, there are also a dissipative and a coherent 
backaction effect [12] (see Sec. IIIBI) . As we discuss below, 
they modify the relative orientations of r° and r 1 and 
therefore affect the qubit decay as well. This mechanism 

first noted in Ref. [32] - derives from a renormaliza¬ 

tion effect induced by the interplay of the readout inter¬ 
actions (A) and the tunneling on and off the sensor QD 
(T) as discussed in Sec. IIIII It results in isospin torques 
similar to those encountered in spintronic QD setups. As 
mentioned in the Introduction (Sec. [T|, the prominent 
role of renormalization effects distinguishes a QD sen¬ 
sor with few, discrete energy levels from sensors with a 
continuous energy spectrum. 










Finally, we note from Eq. m that the Bloch vector 
may not just shrink exponentially. We will find that r° 
and t 1 perform different precessional motions (due to 
both B c ff and the coherent backaction being dependent 
on the charge state of the sensor), which implies that 
the component of r° x x 1 along the measurement vec¬ 
tor A also oscillates in time. Thus, the rate of decay 
of the purity is not purely exponential but additionally 
oscillates in time as explicitly confirmed by our analy¬ 
sis in Sec. IIVD1 This illustrates that the motion of the 
charge-specific isospins is closely related to the qubit de¬ 
cay. Accounting for the interplay of their dynamics turns 
out to be the key to set up a correct description of the 
transient qubit dynamics that includes all the different 
types of backaction. 

III. QUBIT-SENSOR QUANTUM DOT 
DYNAMICS 

A. Outline 

The indirect measurement setup introduced in the pre¬ 
vious section poses several challenges for the theoreti¬ 
cal treatment of the measurement backaction. A central 
complication is that the environment of the qubit (the 
SQD plus the electrodes) is not in a simple equilibrium 
state since the detection is done by nonequilibrium trans¬ 
port. But even when specializing to near-equilibrium 
conditions, one has to treat the SQD as a strongly in¬ 
teracting quantum system with spin degeneracy. Both 
the nonequilibrium conditions and the interactions in the 
SQD prevent a simple direct approach where one aver¬ 
ages over the environmental degrees of freedom, leaving 
only the qubit degrees of freedom. Moreover, to obtain 
the sensor current, we need to retain the sensor degrees 
of freedom as well. As we discuss in Sec. IV B1 specifying 
the environmental state is the main difficulty when try¬ 
ing to directly calculate the evolution of the qubit density 
operator for an indirect detection setup. 

Therefore, we integrate out only the electrodes to ob¬ 
tain the density operator p for the joint qubit-SQD sys¬ 
tem. The resulting equation, which is of the form 

j t P{t) = ~iLp(t), ( 22 ) 

is given below [Eq. (1261) ] and is the first main equation 
of this work. Our main conclusion is that this provides 
a systematic description of the measurement backaction: 
in the weak tunneling, weak measurement limit T/T <C 
A/T -C 1, it does not get any simpler without making 
drastic concessions. Yet, a mere reformulation of Eq. (l26l) 
[see Eq. (1561) ] already provides important insights into 
the measurement backaction. 

Still, we will also describe an attempt to eliminate the 
SQD degrees of freedom in the high-temperature limit 
where corrections T /T can be dropped. This results in 
an effective Liouvillian L e g that reproduces the outcome 


of Eq. (l22l) for the evolution of total average isospin op¬ 
erator, r(f) = (r) = r°(t) Vr 1 (t), in the long-time limit, 

= -iL e ftr{t) (£-£ 0 >l/ r )> (23) 

at 

with initial time to- The preparations for this step pro¬ 
vide interesting insights into Eq. (1221) . However, Eq. C3l) 
turns out to be invalid for small times t — to < 1/T. The 
error made when still using Eq. (E3l) to compute r{t) 
starting from t = t 0 for t — to 1/T can be compensated 
by a correction to the initial condition r(fo), a so-called 
initial slip. This correction depends on the initial qubit¬ 
sensor state in an essential way, preventing the sensor 
from being integrated out completely. Before we discuss 
the details, let us first outline the further challenges posed 
in deriving the above two equations. 

The derivation of Eq. (1221) has to include various ef¬ 
fects: First, since we incorporate the measurement back- 
action terms ~ TX/T, we must also include next-to- 
leading order tunnel processes ~ T 2 /T into the kinetic 
equations for the SQD-qubit evolution (see Sec. Ill Bl) . We 
have given such a consistently expanded kinetic equation 
in Ref. (32j. However, there we employed an additional 
Markov approximation with respect to the electrodes, 
which is valid to obtain the stationary long-time limit 
studied in Ref. HI]. Here, by contrast, we are interested 
in the transient dynamics, where non-Markovian effects 
induced by the electrodes must be accounted for, as we 
explain in Sec. IIII B 11 We include the required leading 
non-Markovian correction perturbatively in the tunnel 
coupling T along the lines of Refs. [4(J [4l|, |U f78l - [80l . 
We present and explain the resulting time-local kinetic 
equations in Secs. IIII B 21 ITlI B 41 

We next analyze in Sec. IIII Cl how the qubit is affected 
by the measurement within the resulting description. For 
this purpose, we first solve in Sec. IIII CTl the kinetic equa¬ 
tions for zero capacitive interaction A = 0. An important 
step is to identify a set of quasistationary modes that 
contain the degrees of freedom of the qubit only, i.e., 
t = t° + t 1 . This identification remains valid also for 
nonzero capacitive interaction A ^ 0. The time scale 
~ I/O for the evolution of these modes - connected with 
the slow qubit dynamics - is well separated from that 
for the evolution of the residual decaying modes. Those 
are strongly damped on a short time ~ 1/T due to the 
fast tunneling dynamics of the SQD. We then introduce 
new dynamical variables to analyze the coupling between 
the quasistationary and decay modes in Sec. IIII C 21 for 
nonzero capacitive interaction A. This will reveal the 
mitigation of the measurement backaction by the coher¬ 
ent backaction, the first key result of the paper. Finally, 
we derive the evolution of the quasistationary modes by 
effectively incorporating the impact of the decay modes 
(see Sec. IIII El) . Importantly, the resulting equations are 
not independent of the detector evolution and even in 
the long-time limit an explicit dependence on the initial 
overlap with the decay modes remains as we will see in 

Sec. [TV] 
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B. Kinetic equation 

1. Integrating out the electrodes 

Whenever the time evolution of an open system is con¬ 
sidered, non-Markovian features arise from the memory 
of the environment, that is, its correlationjunctions de¬ 
cay within a nonzero correlation time tq HH, S3 • When 
integrating out the environment (here the electrodes, see 
Fig. HD, the time evolution of the reduced density opera¬ 
tor p(t) of the open system (here the SQD plus qubit) is 
governed by a time-nonlocal kinetic equation: 

p(t) = —HjQ S p(t) + f dt'W(t — t')p(t’). (24) 

Jt 0 

Here, Lqs» = [Hq + H$ + Hj,»] is the internal Liouvil- 
lian of the reduced system with denoting the operator 
the Liouvillian acts on. Moreover, all effects of the envi¬ 
ronment are contained in a kernel W that we compute by 
a real-time diagrammatic approach USES- If the initial 
value p(t 0 ) is specified, Eq. (l24l) can be used to compute 
p(t) without explicitly keeping track of the state of the 
electrodes. A key assumption enabling such a closed de¬ 
scription of p{t) is that the reservoir is stationary, i.e., 
[Hr,ph\ = 0, which is satisfied here because we assume 
the electrodes to be a thermal equilibrium state pR (see, 
e.g., Ref. f50jb 

When calculated to leading order in T, the kernel 
roughly decays as W(t — t 1 ) ~ Te _ ( t_i )/ TC with cor¬ 
relation time tc ~ l/T [49|. Within the Markov ap¬ 
proximation with respect to the electrodes, one replaces 
W(t — t') = W(iO)S(t — t'), where on the right-hand side 
W{z) = J 0 °° e lzt W(t) denotes the Laplace-transformed 
kernel. This yields a time-local kinetic equation dp/dt « 
(—iLtQs + W(i0))p(t) when inserted into Eq. ET1) . In 
general, (—iLqs + W(i0))p s t = 0 determines the ex¬ 
act stationary state, a fact which is often overlooked 
but easily shown [50j. Physically, this makes sense since 


a nearly constant state cannot ‘‘remember” much. As 
pit) approaches the constant stationary state p Bt , non- 
Markovian corrections in Eq. (1251) become weaker [for 
fixed if the memory kernel W(t—t 1 ) decays as t increases]. 

To go beyond this Markovian approximation to obtain 
the transient dynamics, we include the non-Markovian 
corrections induced by the electrodes perturbatively in 
T /T. To do so, we insert the Taylor expansion for the 
reduced density operator, 

p(t') = p(t) + p(t){t' - t) + ..., (25) 

recursively into Eq. (EH) , as explained in Refs. (40l. [4lt| . 
As we argue in Appendix I A 11 the derivatives d n p/dt n are 
on the order of T n , and we estimate (t' — t) n < rg ~ 1/T n 
within the correlation time of the kernel. Thus, higher- 
order terms in the expansion (1251) correspond to higher 
orders in the tunneling expansion in T /T. Truncating 
the expansion after the leading-order memory correction 
(n = 1), one can derive a time-local kinetic equation for 
p(t) as we show in Appendix IA II 

The above treatment is closely related to the tech¬ 
niques developed for full counting statistics f78l - l80| and 
to the recent study in Ref. j47j- There is also a concep¬ 
tual connection to time-convolutionless master equations 
[5345^ : In the latter approach, the full density opera¬ 
tor evaluated at time t’ is obtained by evolving the full 
density operator at time t backwards in time before inte¬ 
grating out the electrodes, resulting also in an effectively 
time-local kinetic equation. 


2. Kinetic equation 

Including all terms of order A, T, as well as T 2 /T and 
AT/T, where A ~ f2, A, as well as the leading memory 
corrections we obtain the kinetic equation expressed here 
in the representation (1131) of p (no time arguments writ¬ 
ten): 


d 

dt 


( p° \ 

/ 

-2 7 ° 

+7 1 

+2cA- 

+cA- \ 

/ P ( 

p 1 


+2 7 ° 

-7 1 

—2cA- 

-cA- 

P 

T ° " 


+2cA 

-|-cA 

—2 7 ° + (fi — kA) x 

+ 7 1 + kA/2x 

T 

Vr 1 / 

V 

—2cA 

—cA 

+2 7 ° + kAx 

— 7 1 + ($7 + A — kA/ 2) x / 

V T 


(26) 


When computing the matrix product with the col¬ 
umn vector in the above equation, the dot (cross 
“ x ”) in the entries of the matrix indicates that a three- 
dimensional scalar (vector) product is to be formed with 
the corresponding entries of p. The above equation is 
valid under the weak-coupling assumption T/T -C A/T -C 
1 introduced in Sec. Ill Bl such that corrections of order 
T 3 /T 2 , T 2 A/T 2 , and TA 2 /T 2 can be neglected. 


The above kinetic equation is the first central equation 
of this paper. It goes beyond a simple master equation 
by including all relevant coherences (see Sec. Ill Cl) and 
extends the kinetic equation of Ref. [33 , which is Marko¬ 
vian with respect to the electrodes, to access the transient 
dynamics by including the kernel frequency dependence. 
The kinetic equation (1261) respects the probability conser- 
vation, p° +P 1 = 0, and also the recently found [32| exact 
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isospin sum rule (1191) . f° + t 1 = fl x r + A x r 1 . The 
latter derives from the conservation of the total isospin, 
t = t° + t 1 , when electrons tunnel from the electrodes 
into the SQD and vice versa, a generic feature jig] of in¬ 
direct measurement models of type ©• We next discuss 
the expressions and physical significance of the four new 
coefficients 7 °, 7 1 , c, and k occurring in Eq. (f26]l : for the 
definition of A and L2 see Eqs. © and 13]), respectively. 


3. Stochastic, dissipative, and coherent backaction 

First, Eq. (l26l) incorporates the SQD switching rates 
7 0,1 = E r r )^ 1 with contributions from each junction 
r = s,d reading 


,,0,1 


= E 


r=s,d 


\fr + 


\ ' r r r„ 


/ y 2 T 
q—s,d 


(/r ± )>, 


t E ^;(2 /+ + /-) 


(27) 


where 0,1 corresponds to ±. Let us first focus on 
the meaning of the three different physical terms in 
Eq. 127]). The first term in the first line of Eq. m 
is the sequential tunneling contribution, whose depen¬ 
dence on the voltages is governed by the Fermi func¬ 
tions = / ± ((e — t l r)/T) for electrode r = s,d with 
f + (x ) = l/(e x — 1) and f~{x) = 1 — f + {x). We com¬ 
ment on the non-Markovian correction factor 77 [Eq. (1331) ] 
in Sec. IIII B dl The second term is a correction to the se¬ 
quential tunneling rate accounting for a renormalization 
of the level position e, incorporating the derivative of the 
Fermi function, 


(. fr)' = 


dfHx) 


dx 


x=(e— fi r )/T 

and the renormalization function, 


<Pr = (j>((e - Hr)/T), 


with 


/*+A 

= v / 

J -A 


dy f + (y) 

-A 7T x-y 

- Re HE^) +ln (^ 


(28) 

(29) 

(30) 

(31) 


Indeed, combining this second term with the first, 
/±(e) + (f±)'(e)(6/T) « f±(e + S), one identifies the 
shift S = E g IW 2 . The function (j>{x) is plotted in 
Fig. [2] and shows a maximum at x = 0 with logarithmic 
tails. In Eq. (l30l) . V denotes the principal value of the 
integral with a cutoff A = W/T, yielding the real part 
of the digamma function if) with a logarithmic correction. 
The latter depends on the electrode bandwidth W, which 


must be set to W ~ U, where U is the large but finite lo¬ 
cal Coulomb interaction energy of the SQD (we excluded 
the doubly occupied state from the SQD Hilbert space). 

The term in the second line of Eq. (03 relates to the 
cotunneling processes through the SQD, which incorpo¬ 
rates the derivative of the renormalization function, 


d(j>{x) 

dx 


c=(e-pr)/T 


(32) 


which is also plotted in Fig. [2] The contribution from 
each electrode r changes its sign close to the resonance 
e = fi r and takes its extremal values of 4> r r « =f 0.143 at 
e — y, r ~ ±1.911T. While the terms in the first line of 
Eq. (1271) depend exponentially on the distance to the res¬ 
onance |e — /i r |, th e cot unneling term is only algebraically 
suppressed, since |l 00 | 


<t>' r ~ — — ——— for |e — /j, r \ T. (33) 

7T e — n r 

and 2/++/“ = /++1 ^ 1 in Eq. (1771) . When these terms 
are added together, they result (for each electrode r) in a 
temperature-broadened step function, which approaches 
its asymptotes algebraically. Therefore, this must be ac¬ 
counted for when studying the qubit-sensor dynamics at 
the onset of Coulomb blockade where typically the read¬ 
out is performed. 



0.20 
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FIG. 2: Renormalization function 4>{x), Eq. (1301) . and its 
derivative 4>'{x ) = drf>(x)/dx. For illustration purposes we 
chose A = 15 to be rather small, noting that A only shifts the 
<j) vertically and drops out in (f>' [see Eq. 42D]. 


All the above-mentioned tunneling processes contribute 
to a stochastic switching of the SQD occupation n, which 
results via the capacitive interaction Hi = nX ■ t/2 in 
a fluctuation of the effective field B” ff = Cl + nX acting 
on the aubit as exDlained in Sec. Ill Cl [see Ea. ©]■ 
The importance of the capacitive interaction to produce 
this stochastic contribution to the total measurement 
backaction becomes apparent when rewriting the kinetic 
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equation (1261) in terms of quasistationary and decaying 
dynamical variables (see Sec. IIII C 21 ). It causes the 
decoherence of the qubit already in lowest order as we 
show in Sec. m 

Before we enter the detailed analysis of the next-to- 
leading order corrections, let us right away indicate 
their importance for the stochastic backaction on a more 
qualitative level. In Fig. [3l we compare the evolution 
of the x-conrponent r x (t) of the isospin, obtained by 
solving the kinetic equations m, when higher-order 
terms are included (red) or neglected by hand (green). 
The figure illustrates that noise from 0(r 2 ) terms 
indeed contributes to the qubit decoherence as naively 
expected. On a quantitative level, however, one would 
expect an algebraic suppression of the measurement 
backaction with e based on Eq. (1271) when entering the 
Coulomb blockade regime. It will turn out in Sec. IIII C 31 
that this expectation is incorrect, i.e., the backaction is 
weaker than expected. 

We further see from Fig. [3] that the oscillation period of 
the qubit is notably changed due to next-to-leading order 
corrections. This is due to the mean field, $7 = f2 + (n)A, 
acting on the qubit in the presence of the sensor QD. 
The average occupation (n) = pl t on the sensor is 
significantly modified by higher-order tunneling terms 
(see Sec. IIVB1I) . 


In addition to the stochastic backaction, there is also 
a dissipative backaction of the SQD on the qubit: These 
terms are related to the isospin-charge conversion rates 
~ cA with coefficient 

(34) 

r 

This coupling appears in two ways. The isospins influ¬ 
ence the SQD dynamics (allowing for the readout) and 
vice versa the SQD occupation probabilities directly in¬ 
fluence the isospins (backaction). This dissipative back- 
action drives the sensor QD and qubit into a correlated 
state in the stationary limit, in contrast to the stochastic 
backaction, which only changes the occupation probabil¬ 
ities of the qubit state. For the parameters chosen for 
Fig. [31 the dissipative backaction has a negligible impact 
on the qubit evolution and therefore we show no compar¬ 
ison. The reason for this suppression is that the SQD is 
already mildly Coulomb blockaded for these parameters 
and the dissipative backaction is exponentially peaked 
around the resonance as Eq. (I3H) shows. The dissipative 
backaction therefore only becomes relevant close to reso¬ 
nance. 

Finally, there is a third type of backaction: the tunneling 
gives rise to isospin torque terms ~ kA, where 

K = ( 35 ) 

r 



FIG. 3: Modification of the sensor backaction due to the 
next-to-leading order T 2 stochastic backaction and the co¬ 
herent backaction. Shown is the x-component of the isospin, 
T x {t) = r x (t) + ri(t) as a function of time. The isospin is 
obtained by solving the kinetic equations (1261) when all terms 
are taken into account (blue), when higher-order T 2 correc¬ 
tions are neglected (green), or when the coherent backaction 
oc k is neglected (red). One can see that the isospin precession 
period is larger if next-to-leading-order contributions are ne¬ 
glected since the isospin experiences a different mean field as a 
consequence of the differing average occupation (n). The time 
is given for all three curves in the same unit, the inverse of 
the natural frequency for the full result, fi = ^/Q 2 + ((n)A) 2 , 
which is the precession frequency including next-to-leading 
corrections. We insert here ( n) = p,[ t = 27 °/y from the sta¬ 
tionary solution of the full kinetic equation OH) at A = 0 [see 
Eq. (1491) 1. The parameters are lb = Td = 0.2 T, A = 0.02T, 
fi = 0.001T, V b = 0, Vg = —3T, W = 1000T. Initially, the 
sensor is empty, p°( 0 ) = 1 — p 1 ( 0 ) = 1 , and conditional upon 
this, the qubit isospin vector is prepared perpendicular to the 
measurement vector, r°( 0 ) = e x and t 1 (0) = 0 . 

incorporates the derivative of the renormalization func¬ 
tion shown in Fig. [2] This signals the coherent nature 
of this contribution to the backaction: it characterizes 
the response of the sensor QD level renormalization to a 
change in the state of the qubit. Similar to the cotunnel- 
ing corrections, the coherent backaction gains importance 
with the onset of Coulomb blockade [321 ]. 

The importance of coherent backaction for the qubit evo¬ 
lution stands out in Fig. [3] We can see that neglecting 
the coherent backaction (red) noticably (but artificially) 
enhances the qubit decoherence as compared to the full 
solution (blue). This points to a cancellation effect be¬ 
tween coherent backaction and cotunncling noise that we 
discuss in detail in Sec. nr C 31 We further note that the 
coherent backaction has a negligible effect on the qubit 
oscillation period [the period is the same for both the 
curves including 0(T 2 ) corrections]. The coherent back- 
action can thus not be interpreted as a simple correction 
to the qubit mean field (EH); it is the joint system of 
qubit and sensor QD that is renormalized and not just 
the qubit system. We finally emphasize that even though 
Fig. [3] shows theoretical results when different contribu- 





















12 


tions of the kinetic equation are neglected, they cannot 
be switched off individually in a real experiment - there 
they always appear together and have to be taken into 
account altogether. 


4- Impact of non-Markovian corrections 

It remains to discuss the three ways in which non- 
Markovian corrections induced by the electrodes are con¬ 
tained in Eq. (l26l) and how the latter differs from the 
Markovian kinetic equations (Eq. (18) of Ref. [US]). 
First, the non-Markovian corrections modify the leading- 
order SQD tunneling rates [see Eq. (12711 ] just by introduc¬ 
ing the prefactor 


77 = l + £) 


T 


(36) 


Since the correction 77 — 1, the cotunneling broadening 
term [see Eq. jm and the coherent backaction coef¬ 
ficient k [see Eq. (1551) ]. all depend on the same factor 
^ r r /T(f>' r with algebraic tails, the non- Markovian ef¬ 
fects should clearly be accounted for |lOl| . The correc¬ 
tion factor rj — 1 is an appreciable quantitative correc¬ 
tion that yields a contribution of 0(T 2 /T) to the switch¬ 
ing rates [4(j ■ However, in contrast to the cotunneling 
and coherent backaction, the correction r/ — 1 is mul¬ 
tiplied with the exponentially scaling SET contribution 
[see Eq. (1271) ] and therefore it has no qualitative impact 
here. 

By contrast, the second type of non-Markovian cor¬ 
rection affects the coherent backaction terms by qualita¬ 
tively changing them in general relative to the Markov 
approximation. The direction of the tunneling-induced 
isospin torque terms changes: In Ref. [32], we also found 
a contribution to the coherent backaction oc n fl. These 
terms are canceled out here up to 0(TQ/T) = 0(TX/T). 
This is expected on physical grounds as the backaction 
is mediated by the capacitive interaction A and there¬ 
fore we expect these terms to vanish when setting A = 0. 
We emphasize that the isospin torque terms oc ntl x r” 
do not affect the stationary state, which we studied in 
Ref. [13 , and therefore all the conclusions drawn in Ref. 
■32] remain valid. The third effect of the non-Markovian 
corrections is a sign change of the isospin torque terms 
~ k\ in the last column of the matrix (1261) as compared 
to Ref. [HI. 

Both the above modifications of the coherent backac¬ 
tion have important physical consequences illustrated in 
Appendix IA 21 If one naively computes the transient dy¬ 
namics of the SQD-qubit system using the equations of 
Ref. [13 , which neglect non-Markovian terms induced by 
the electrodes, one obtains exponentially increasing tran¬ 
sient modes leading to a violation of the positivity of the 
density operator. Moreover, within the Markovian ap¬ 
proximation the coherent backaction strongly enhances 
the measurement backaction in the Coulomb-blockade 


regime for a large parameter regime, while the coher¬ 
ent backaction suppresses the measurement backaction 
for nearly all parameter values when non-Markovian cor¬ 
rections are correctly accounted for (see Appendix I A 21) . 
This clearly illustrates that non-Markovian corrections 
go hand in hand with renormalization effects, which in 
an indirect measurement set up go hand in hand with the 
cotunneling effects of the sensor rates. All of these are 
of vital importance for describing the indirect measure¬ 
ment. 


C. Coupling of modes 


With the kinetic equations (1261) now in hand we can 
proceed to analyze the measurement backaction, but still 
without integrating out the sensor. To achieve this goal, 
we make use of the separation of different time scales 
in the coupled evolution of SQD and qubit in the weak- 
coupling, weak-measurement limit T/T <C A/T <C 1. To 
identify these time scales, we first solve in Sec. IIIIC li the 
unperturbed problem of the decoupled SQD-qubit system 
(A = 0) as described by Eq. (1251) . This produces eigen- 
modes which are well-separated in energy by T ( 3 > H and 
correspond to the wide-band limit for the sensor quan¬ 
tum dot. It turns out that one needs to compute the 
evolution of only a part of the modes - referred to as the 
quasistationary modes in the following - to construct the 
evolution of the total isospin t. 

In Sec. IIII C 21 we restore the coupling A and by sim¬ 
ply writing the kinetic equation (1261) in the basis of these 
eigenmodes, we can immediately extract A 2 /T as the rele¬ 
vant time scale for the qubit decoherence time. However, 
this is not the full story of the backaction: there is a 
prefactor which strongly affects this time scale. We ana¬ 
lytically identify a nontrivial competition of the coherent 
backaction and the cotunneling-induced stochastic back- 
action determining this prefactor. Finally, introducing 
an exact [relative to Eq. (l26l) ] projection of the dynam¬ 
ics onto the quasistationary modes in Sec. IIII El we gain 
further insight, still without integrating out the sensor. 
This projection incorporates the effect of the coupling 
between the modes and provides the starting point for 
deriving an effective equation for the isospin evolution in 
Sec. El 


1. Quasistationary and decay modes for A = 0 


We first solve the kinetic equation (l26l) for A = 0 in 
which case the dynamics of the occupation probabilities 
p° and p 1 decouples from the dynamics of the isospins t° 
and t 1 as shown by the “unperturbed” time-evolution 
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e° 



pit — quasi- 
p°t — stationary 
pO p i modes (q) 


P,t 

P°t 


decay 
modes ( d ) 


charae modes (c) isosDin modes 


FIG. 4: Sketch of the complex eigenvalues and the associated 
dynamics of the joint qubit plus sensor QD system without 
readout (A = 0). The four upper sketches depict the evolution 
of the total isospin r = r° + r 1 described by the coefficients of 
the four quasistationary modes in the expansion of p [see Eq. 

Moreover, the lower four sketches depict the evolution 
of the weighted difference <5 = Pst T ° — Pst " 1 " 1 of the charge- 
specific isospins associated with the coefficients of the decay¬ 
ing isospin modes [see Eq. (1451) ]. The indices a = +1,0,—1 
label the three different polarizations for the precessional mo¬ 
tion. The unit vector eo = f2 is the “bare” qubit field and 
the vectors ei and e 2 can be chosen arbitrarily in the trans¬ 
verse plane. The two sketches on the right show the different 
dynamics of the occupation probabilities p° and p 1 for the 
stationary and the decaying charge mode, respectively. Note 
that p° and p 1 are not the coefficients in Eq. (1451) but only 
their combinations p° + p 1 = 1 and 5 S t = pltP° — PstP 1 ■ If 
the readout is included (A ^ 0), the modes become coupled. 
Since L 7 ^ id, the transitions between the modes are not al¬ 
ways possible in both directions as indicated by the arrows. 

generator Lq := L|,\-o: 

/_2 7 ° + 7 i 0 

, r +2 7 ° - 7 ! 0 

lL ° 0 0 -2 7 ° + f2x 

\ 0 0 +2 7 ° 


This can be brought into diagonal form easily by noting 
that the cross product operation is diagonal in the basis 
of the complex unit vectors [Eq. (IB.301) ] 

eo = fl = fl/fl, (38) 

e± = (ei + ie 2 ) /a/2. (39) 


—y 1 + fix 


(37) 


we find 


- iL 0 


E [*n p v*v«t + (in p 


7 )VpVp' i ’j , (41) 


with p = c, 0 , ±, and eigenfrequencies 


n c = n 0 = o, n ± = ± n, (42) 

and V£t and V k are the left and right eigenvectors, re¬ 
spectively, using dyadic notation. The indices p = c, 0, ± 
and k = q,d label in total eight different modes. Be¬ 
fore we discuss the physical meaning of these modes, we 
first note the following general property: since Lq is di- 
agonalizable (although ^ Lq), all the left and right 
eigenvectors are mutually biorthonormal, 

v k ;-v k p : = s kk 's pp ,, (43) 

and they satisfy the completeness relation 

E E v p E f - = !• ( 44 ) 

k—q,d p=c,0,± 

This also follows explicitly by using Eq. (Ha and the 
expressions below. This can be exploited to expand the 
state vector p{t), defined by Eq. (fl3l) . as follows: 


p = (p 0 +P 1 )V« + ^(r°+ri)V« (45) 

a 

+(pltP° - PsV)E + E - p°sy a ) v d a 

a 

v* + E e<not «S(o)v« (46) 

a 

+ e~^{0)Vi + e (in ““ 7 )t ^(0)V[(, 

a. 

where r™ := • r [see Eqs. (l38l) and m] and the 

coefficients are given by 

*£(0) = V^-p(O), (47) 

taking the initial time to = 0 here. Importantly, equality 
m is generally valid for any state p(t), while the second 
equality (l46l) holds only if pit) = e~ lLot p{ 0). 

We now discuss the explicit form of the modes. The 
most fundamental one is the stationary charge mode 
with the conjugated left eigenvector and the right eigen¬ 
vector 


These are constructed from a right-handed orthonormal 
system (eo,ei,e 2 ) with eo = and an arbitrary choice 
of unit vectors eq and e 2 in the plane perpendicular to 
eo. The complex unit vectors satisfy the orthonormality 
and completeness relations (a = 0 , ±): 

e f a -e a > = 5 aa >, = 1, (40) 


V 9 

C 



V 9 

c 



(48) 


expressed in the occupation probabilities of the SQD in 
the stationary limit and for zero coupling A = 0, 


where here the dot • denotes the scalar product taken 
with the object to its right. Writing Lq in diagonal form, 


Pst 



Pst 


2 7 ° 

7 


(49) 
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introducing the often recurring rate combination 

7 := 2 7 °+7 1 - (50) 


as seen by inserting Eq. (l27l) into Eq. (l50l) and noting 
r/T -C 1. The conjugated left and right eigenvector of 
the charge decay mode read, respectively: 


The right zero eigenvector corresponds to a physical 
state, a valid density operator, which is factorizable, 
pi = ps , s t <8 |l. In this state, the SQD is stationary, 
ps, st = J^nPs- t-P”> and the qubit is in the completely 
mixed state with zero Bloch vector (r° = t 1 = 0, c.f. 
Fig. 01 upper left). Any valid solution of the A = 0 kinetic 
equation pit) = —iLop{t) always involves this stationary 
charge mode superposed with other modes. These addi¬ 
tional modes contain the isospin precession as we explain 
in the next paragraph. As one can see from expansion 
(SSJ), the coefficient v%(t) = 1 for all t and irrespective of 
the initial condition because the corresponding left zero 
eigenvector is just the trace operation, guaranteeing that 
pit) has unit trace for all times t: 


/+ X \ 

0 ■ (53) 

V 0 ) 

If only this mode is excited in addition to the funda¬ 
mental stationary charge mode, the SQD state devi¬ 
ates from the stationary state, i.e., the coefficient 6 s t := 
pl t p° — PstP 1 7 ^ 0 hr Eq. (H5l) . while the qubit remains in 
the completely mixed state (t° = r 1 =0). Clearly, such 
deviations from the stationary SQD state decay on the 
short time scale 1/7 set by the SQD tunneling dynam¬ 
ics (see Fig. H. Finally, there are three decaying modes 
with conjugated left and right eigenvectors, respectively, 
(a = 0 , ±) in which the isospins are polarized: 


/ +Plt 


V“ = 


~P°t 

0 

0 


V = 


tr (Pit)) = £>" = Vf -Pit) = 1. (51) 

n 

The remaining seven “modes” have zero trace (see 
Eq. (1431) with k = q, p = c) and therefore cannot repre¬ 
sent proper density operators on their own. These modes 
cannot be excited alone: They always appear in combi¬ 
nation with the stationary charge mode. In this respect 
these modes differ from modes encountered in, e.g., pure- 
state unitary evolution problems. 

There are three more quasistationary modes ik = q), 
for which the SQD remains in the stationary state ps , s t 
but the qubit state is not completely mixed, i.e., the 
isospin r is polarized (see upper right of Fig. H. The re¬ 
lated conjugated left and right eigenvectors, respectively, 
read: 


V« 

v a 



( 


0 

0 


Pst e « 

V Plt e a 


(52) 


The expansion fl5l) shows that the coefficients of these 
quasistationary isospin modes are connected with the to¬ 
tal isospin t = r° + T 1 . If the mode a = 0 is excited, the 
total isospin r points along the qubit axis eo and does 
not precess. If the other two modes a = +(—) are ex¬ 
cited, the total isospin r precesses (counter) clockwise in 
the plane perpendicular to eo with frequency Cl. Thus, if 
the isospin was nonzero initially, expansion (1451) involves 
at least one of these three modes in addition to the qua¬ 
sistationary charge mode. In the case A = 0, the former 
are not damped and the magnitude of the total isospin 
remains unchanged, reproducing exactly the free unitary 
evolution of the qubit. 

In addition to these four (quasi)stationary modes, 
there are four more decay modes ik = d) that are expo¬ 
nentially damped in time. As Fig. [I] illustrates, the eigen¬ 
values of these modes are well-separated from the quasis¬ 
tationary modes in the complex plane since Cl <C 7 ~ T 


v d 

y a 


+Ps t e a I ’ 
V -p° t e a ) 


\ d 
* a 


( 0 

0 

+e a 
V -e a 


■ (54) 


The coefficients of these modes characterize the decay of 
the weighted difference £ st = pl t T ° ~Pst Tl °f the charge- 
specific isospins [see Eq. (H5l) ]. The weighted difference 
5 s t and the sum t = t°-|-t 1 are linearly independent and 
together uniquely determine r° and r 1 . Note that <5 st in¬ 
volves the stationary occupation probabilities in contrast 
to 8 defined by Eq. m3. Similar to the total isospin for 
the quasistationary modes, the difference <5 st can point 
along eo = Cl without any precessional motion for a = 0 
or it can precess (counter)clockwise in the plane perpen¬ 
dicular to e 0 for a = +(—) (c.f. lower right of Fig. UJ). 
We see now that for A = 0, the total isospin r decouples 
from the motion of all other degrees of freedom and par¬ 
ticularly also from 8 S t. In contrast to r, the difference 
5 s t is strongly susceptible to the tunneling dynamics of 
the SQD, i.e., the switching between charge states n = 0 
and n = 1. This generates the noncollinearities of t° and 
t, which are responsible for the qubit decoherence for 
nonzero capacitive coupling [see discussion in Sec. Ill Cl 
Eq. (fl8l) ff.]. What is important for the following is 
that that the expansion (H5l) carries over to the case of 
nonzero coupling A / 0; even in this case, it suffices to 
compute the evolution of the quasistationary modes to 
obtain the evolution of the total isospin [ 102j |. This ob¬ 
servation provides the starting point of the subsequent 
discussion where we first investigate how the quasista¬ 
tionary modes - “containing” the qubit dynamics - are 
coupled to the decay modes and then even explicitly elim¬ 
inate the decay modes (except for their initial state). 


2. Coupling of quasistationary and decay modes 

We now turn the capacitive interaction back on, A 0, 
and investigate what we can say about the evolution us- 
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ing the above discussion of the eigenmodes of the Liou- 
villian L 0 . We can identify new dynamical variables that 
characterize the evolution in the quasistationary and de¬ 
cay subspace, respectively. Based on Eq. (H5l) . we intro¬ 
duce 


X 9 = [ P °+ P \ ), X d = 


pltP° ~ PstP 1 
Plt T ° ~ Pst T± 


, (55) 


and rewrite the kinetic equation as 

d_ f X 9 \ _ / L 99 + A" A 9d \f X 9 \ 

dt V X d ) 1 \ A d9 L dd + A dd . ) \X d ) ' 

(56) 


modes into the decay modes and back. These virtual 
processes are responsible for the qubit relaxation and de¬ 
coherence; the multiplet of quasistationary eigenvalues 
in Eq. m acquires an imaginary part fl a —> fl a + i^ a 
for a = 0 , ±, which induces a shrinking of the isospin 
Bloch vector. From Eq. (l5f>l) and Fig. [2 we expect that 
1/Tj,l/T 2 ~ A 9d r -1 A d9 oc A 2 /E, i.e., if the couplings 
A 9d and/or A d9 are small, then so will the backaction 
be. Before we make this more precise [see Eq. (l69l) ]. we 
investigate the detailed form of these couplings, which 
contains the first main result of the paper. 


3. Mitigation of cotunneling noise by coherent backaction 


This is the second main equation of our study. The A kk 
blocks are given and discussed below except for A dd , 
which is not needed here [it is given by Eq. (IB.34I) in 
Appendix [B]. First, the action of the unperturbed Liou- 
villian on these variables reads trivially 


The transition matrix from the quasistationary modes 
into the decay modes reads [see Appendix IB 3[ 


A d q 


( 0 (!+Pst) cA -\ 

V( 1 + P°t) cA ) 


(61) 


-*L 99 


0 0 \ 

0 fix J ’ 


—iL 


dd 

0 


-7 0 \ 

0 — 7 I + fix J ’ 


(57) with the transition factor 

(58) r = PstPlt - K f -jPlt ~P °t J > 


(62) 


and L qd = Lq 9 = 0 for the A = 0 solution. The perturba¬ 
tion A has two effects. It first introduces a direct action 
on the quasistationary variables: 


— iA qq 


(0 0 \ 
V 0 Pst Ax J 


(59) 


This produces the mean-field backaction, which amounts 
to a tilting of the internal qubit field fl as anticipated 
in Sec. IIVBI adding A" to Lg 9 gives the effective qubit 
field 


fl = fl+pl t A. (60) 

The term A", i.e., the term linear in A thus does not 
lead to dissipative dynamics of the isospin contrary to 
what one might naively expect. The isospin decoherence 
is at least quadratic in A (see below). We also note that 
the effective qubit field fl is modified by 0(T 2 /T) pro¬ 
cesses (cotunneling broadening, level shift), which affect 
the stationary occupation probability pl t [see Eq. (117)1) ]. 
Moreover, if one is close to resonance, the probability pl t 
may be a sizable fraction of 1 , which it approaches in the 
Coulomb blockade regime for e <C p s , Pd- Since we allow 
for A ~ fl, this implies that the correction pl t A leads to 
a large change in the qubit frequency and the direction 
of the qubit axis for a large range of parameters. 

The second effect of the perturbation is due to the 
coupling of the quasistationary variables to the decay¬ 
ing variables due to the off-diagonal blocks A qd and A dq . 
As a consequence, the decaying variables cannot be just 
ignored after a time ~ I /7 since they are permanently 
excited by virtual transitions from the quasistationary 


while the transition matrix back into quasistationary 
modes is given by 

A- = ( ” _“ x ) • (63) 

We first note that A, like Lo, is non-Hermitian, A^ ^ A, 
since the qubit-sensor evolution is nonunitary due the 
tunneling. As a consequence, the two transition matrices 
are markedly different: while transitions from the decay 
modes to quasistationary modes (A 9d ) do not depend on 
the parameters of the SQD (level position, bias voltage, 
and tunneling rates T), transitions from the quasistation¬ 
ary modes (A d9 ) exhibit a strong dependence on the sen¬ 
sor QD parameters that we discuss below. That A qd is 
entirely induced by the readout interaction Hi can be 
shown to be a consequence of the probability conserva¬ 
tion together with the conservation of the isospin during 
tunneling processes. The latter is specific to the indirect 
measurement setup. 

With Eq. (1611) in hand, we can now precisely pin¬ 
point what we mean by stochastic backaction: the di¬ 
agonal component of A dq can be split into a first term, 
p° t pl t A, associated with the stochastic backaction, and a 
second term, — n(p\ t /2 — p° t )A, associated with the co¬ 
herent backaction as signaled by the factor k. The com¬ 
bination p^tPlt appears as a simple consequence of the 
charge fluctuations of the SQD, which are characterized 
by (n 2 ) — (n) 2 = (1 — ( n))(n) = p® t pl t for a two-level 
system (see Sec. IV AD . The rates 7 0 ’ 1 determining p° t Ps t 
incorporate both the effect of the SET tunneling as well 
as that of next-to-leading order corrections. The stochas¬ 
tic term is multiplied with A because to act back on the 








16 


qubit, the “tunneling noise” has to act together with the 
internal interaction Hj = hX • t/2 to evoke a A-induced 
transition mediated by A dq . 

The most striking finding is that A dq is strongly sup¬ 
pressed when tuning the SQD towards Coulomb block¬ 
ade. To see this, we first note that Eq. EH) in¬ 
corporates the isospin-to-charge conversion rates EH. 
c ~ which depend on the derivative of 

the Fermi function. These rates are thus exponentially 
suppressed in the Coulomb blockade regime. We next 
inspect the diagonal element of A dq . It is useful to first 
consider the expansion of the transition factor r to zeroth 
order in T /T : 


2T+W /T\ 

( 2 r+ + r-) 2+ Viy ’ 


(64) 


with T± := £ r r r /± and 2T+ + T~ > F := £ r r r . The 
terms in Eq. (1641) derive only from the stochastic part 
in Eq. EH- Thus, in the single-electron tunneling 
approximation, the factor r is exponentially suppressed 
with gate voltage since either p® t or pl t becomes exponen¬ 
tially small when going off-resonance. One would expect 
that this exponential dependence is removed by includ¬ 
ing cotunneling corrections and the coherent backaction, 
which scale algebraically with e — p r and start to domi¬ 
nate over the single-electron tunneling rates as one moves 
into the Coulomb blockade regime. In our calculation, 
we include these terms as well, but still obtain an ex¬ 
ponential suppression of the transition factor. Indeed, 
expanding Eq. (1621) to the next order in T/T, we find 


r 


1 


2r+T“ 1+ 


Er r 


r'Pr 


( 2 T+ + r-) 2 

+ (r~ - 2 r + )(r~ + r + ) 


T 


+\/F- i P+'l ^ + ) / Er 


2 T+ 



(65) 


with r = dv+{e)/d£ = E r r r (/ r + )'/T. Clearly, T+, T~, 
and T' are determined by the Fermi functions and their 
derivatives. Thus, transitions from the quasistationary 
modes into the decay modes become exponentially small 
when tuning the SQD into the Coulomb blockade regime. 
What this implies is that any deviation from the expo¬ 
nential suppression must be due to even higher order 
tunneling contributions (i.e., beyond cotunneling) and 
thus must be a higher power law ~ l/(e — p r ) n with 
n > 2. The experimentally important conclusion that we 
can draw from this is that the sensor can be switched off 
better with the gate voltage than naively expected [see 
Eq. EH)]. 

What has happened in Eq. (1551) is that the coherent 
backaction k ~ J2 r ^rfir/T oc l/(e — p r ), which depends 
algebraically on e, has completely canceled out the alge¬ 
braically scaling cotunneling corrections to the stochas¬ 
tic backaction in the first term of Eq. EH- Hence, in 
Eq. Ml . the coherent backaction term (TA/T) counter¬ 
acts the change in the stochastic backaction term due 


to a change in the sensor QD rates by the cotunneling 
(A/T x T 2 /T ~ TA/T). This can be seen by explic¬ 
itly comparing Eq. (16-H) to Eq. (lt)5l) and is clearly visible 
in Fig. [3] We emphasize that for this cancellation also 
non-Markovian effects induced by the electrodes are im¬ 
portant, which modify the coherent backaction (see Sec. 
he B 41) . Without these, the transition factor r exhibits a 
different dependence on the level position e that can lead 
to a violation of positivity of the qubit-SQD density oper¬ 
ator as we discuss further in Appendix [A] Moreover, the 
cancellation does not imply that the backaction is pre¬ 
cisely the same as when only accounting for the lowest- 
order T approximation (see Fig. 0 ; the renormalization 
of the level position and non-Markovian contributions, 
both scaling exponentially, still modify the backaction. 

By formulating the problem in Eq. (1561) in terms of 
the A = 0 eigenmodes one most clearly sees how the 
cotunneling and coherent backaction, formally terms of 
different order, conspire to effectively cancel out. Note 
also that the dissipative backaction (through c) appears 
on its own. This highlights the importance of keeping 
track of all three types of backaction that are revealed 
only after integrating out the electrodes coupled to the 
sensor QD to obtain our central Eq. (l26l) . It should be 
noted that the dissipative backaction couples the qua¬ 
sistationary modes to the decay modes and therefore is 
not relevant for the leading-order A 2 /T dephasing times 
as Appendix IB 3 bl shows; see also the discussions after 
Eq. (RHl) and after Eq. (11111) below. Rather, the dissi¬ 
pative backaction must be kept to be able to calculate 
the response of the dissipative sensor current of which it 
represents the flip side, as explained in the introduction. 
As emphasized in Sec. IIII B 31 we were careful through¬ 
out our analysis to include all terms which depend on the 
function <f>' r with algebraic tails that could possibly can¬ 
cel out. In Appendix IB 41 we further discuss the cancel¬ 
lation in view of our weak-measurement, weak-coupling 
assumption (see also Sec. Ill Bl) . 

An important conclusion, which we draw in Sec. 13 is 
that this cancellation of cotunneling noise and coherent 
backaction cannot be understood within simple classi¬ 
cal fluctuator model. Although this approach could, in 
principle, be extended to account for the cotunncling- 
induced noise by modifying the switching rates, it seems 
not possible to include the coherent backaction. More¬ 
over, other approaches that aim at directly calculating 
the qubit Bloch vector t must make an assumption about 
the qubit environment, in particular the sensor QD. Here 
one is liable to miss the above cancellation as we also dis¬ 
cuss in Sec. IVB2I 

It is furthermore interesting to observe that this cancel¬ 
lation appears even though the coherent-backaction in¬ 
duced torque terms in the kinetic equations (02) scale 
with A, while the cotunneling corrections do not. How¬ 
ever, to affect the qubit, the “cotunneling noise” has to 
act together with the internal interaction Hj ~ hX ■ t 
to evoke a A-induced transition mediated by A dq . This is 
why they both affect the measurement backaction to first 
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order in A. Note also that the transition factor r is not 
only independent of the SQD-qubit interaction A (it ap¬ 
pears as a factor in A dq ) but also of the internal qubit held 
17. This means that the relative importance of coherent 
backaction over the stochastic backaction cannot be al¬ 
tered by measuring weaker or stronger (i.e., by changing 
A); in the weak measurement limit, these effects phys¬ 
ically come together and should be calculated together 
[103l |. Other experimental parameters alter this compe¬ 
tition and their effect is studied in the next section. 

In short, the delicate interplay of the qubit plus sensor 
renormalization oc T and sensor cotunneling rates oc T 2 
in indirect capacitive detection may be rationalized as 
follows. By keeping track of the sensor-quit coherences 
(since both are quantum systems), we find that coherent 
effects counteract decoherence, which is not really that 
unexpected. This may in fact present a key difference 
of a sensor with quantized levels from a single-electron 
transistor with a continuous spectrum. A comparison 
of both types of detectors regarding the importance of 
renormalization effects is therefore an interesting future 
task. Finally, we emphasize that for the above cancella¬ 
tion the modeling of the qubit as a charge qubit is not 
relevant as long as the isospin couples to the charge of 
the SQD. 


D. Experimental control over backaction strength 

Since Eq. (1651) is a key result of this paper, character¬ 
izing - together with c - the strength of the backaction 
[beyond the mean-field effect in Eq. (K>01) ]. we now in¬ 
vestigate its dependence on experimental parameters in 
some detail. 


1. Gate-voltage dependence of transition factor r. 


dephasing rates provided the sensor can be made the 
dominant environment (which should be the case for a 
good qubit, for which noise from manipulation “chan¬ 
nels” can be switched off). Measuring the qubit deco¬ 
herence rates could clearly distinguish between an al¬ 
gebraic and exponential dependence in an experiment. 
Here, we expect the measurement-induced decay rates 
to scale exponentially into the Coulomb blockade regime 
until higher-order tunneling processes at least of order 
T 3 /T 2 become relevant. They can lead to a crossover 
to an algebraic scaling oc l/|e — p, r \ n with n > 1 deep 
in the Coulomb blockade regime, see Appendix IB 41 In 
any case, our numerical examples illustrate that a QD 
detector can be switched off more efficiently with a gate 
voltage than naively expected. 


2. Bias-voltage dependence of r 

We continue with the discussion of the bias dependence 
of the transition factor r, which we show in Fig. [5(b). For 
small bias voltages, the SQD is Coulomb-blockaded and 
the coherent backaction strongly suppresses r. When the 
bias is increased, sequential tunneling sets in and when 
the level position £ ss 14/2 is resonant with the elec¬ 
trochemical potential of the drain, the transition factor 
saturates. Here, the correction due to the coherent back- 
action actually becomes positive , as shown in the inset 
in Fig. [5(b). Yet, one should note that in the sequen¬ 
tial tunneling regime the coherent backaction has only a 
small impact. For even larger bias voltages, the correc¬ 
tion from the coherent backaction drops to zero. This is 
in accordance with the general finding that renormaliza¬ 
tion effects can be neglected in the large-bias limit [3lil32 ! l] 
because re = J2 r ^r<P' r / T ~ X) r IV/(£—/v) oc 1/14 is sup¬ 
pressed [see Eq. @]. 


In Fig. [5(a), we plot the dependence of the transition 
factor r on gate voltage. The figure shows two curves, one 
including the coherent backaction (blue, n yf 0) and the 
other excluding the coherent backaction (green, k = 0). 
Clearly, the factor r is the largest in the sequential tun¬ 
neling regime of the SQD (| 14 1 < 14/2). Here, transitions 
from the quasistationary modes into the decay modes 
are induced by the fast succession of tunneling electrons, 
which impose a strong noise on the qubit. The coherent 
backaction and cotunneling are negligible in this regime. 
This drastically changes when tuning the SQD into the 
Coulomb blockade regime (|14| > 14/2). The full tran¬ 
sition factor r is actually exponentially suppressed with 
gate voltage [linear on the scale of Fig. [5(a)]. By con¬ 
trast, the gate-voltage dependence is markedly nonexpo¬ 
nential when neglecting the coherent backaction, charac¬ 
teristic of cotunneling noise, see also the discussion of Eq. 
(l65l) . 

Experimentally, we expect the cancellation to be reflected 
in the voltage dependencies of the qubit relaxation and 


3. Tunnel coupling dependence of transition factor 
transition factor r 


We finally discuss the impact of the coherent back- 
action when changing the tunnel couplings and their 
asymmetries. Since the coherent backaction is linear in 
k ~ T/T [see Eq. (1551) ], increasing the average tunnel cou¬ 
pling f = (r s +r d )/2 and lowering the temperature both 
increase renormalization effects in a trivial way (within 
the limit T/T -C A/T <C 1). By contrast, the asym¬ 
metry of the tunneling rates in the generic experimental 
situation, quantified by 


9 ■= 


14 -r d 
r s + r d ’ 


( 66 ) 


may have a nontrivial effect. Controlling the asymmetry 
has also been suggested [] 1 | as an experimental strategy 
for optimizing sensor efficiency in the limit r^/T,, > 1 . 
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K = ~e[T\ 



V„[T] 


FIG. 5: Voltage dependence of the transition factor r, 
Eq. (1621) , which determines the measurement-induced backac- 
tion of the sensor QD on the qubit (e.g., decoherence rates). 
We include the coherent backaction in the blue curves (k, ^ 0), 
while we exclude it by hand for the green curves (k = 0). We 
show the dependence of r in (a) on gate voltage V g for bias 
voltage Vb = 2.5 T and in (b) on bias voltage 14 for gate volt¬ 
age V g — 5 T. In all plots we use 14 = Yd = f = 0.02T and 
bandwidth W = 1000T. The inset in (b) shows the difference 
of the two curves for k = 0 and s ^ 0 and illustrates that 
the coherent backaction can also enhance the backaction for 
larger bias voltages. Note that r is independent of A and Yl 
[choosing, e.g., A ~ Yl ~ O.lf would be a parameter combina¬ 
tion consistent with the conditions for the kinetic equations 
(Infill to be valid: f/T <€. X/T •< 1]. 


In the stationary limit, we previously found [32J that 
for rrf/r s < 1 the impact of the coherent backaction is 
strongly enhanced. 

The effect of asymmetries on the transition factor r 
[Eq. (|521) ] strongly depends on the chosen bias and gate 
voltages. To illustrate this point, we plot in Fig. Ha) 
the relative change in r, 


K 


r K =o - r K ^o 
o 


(67) 


due to the coherent backaction as a function of gate volt¬ 
age V g for the different values of g as indicated. When 
introducing a nonzero junction asymmetry j ^ 0, the ex¬ 


ponential suppression with gate voltage effected by the 
coherent backaction, found in Fig. H a ), is simply rigidly 
shifted horizontally without changing its shape consid¬ 
erably. This can be understood from the fact that the 
maximum of r as a function of V g is shifted by an asym¬ 
metry g due to a basic effect of Coulomb interaction jl04j ] . 
In the vicinity of this maximum, leading-order processes 
dominate and cotunneling and the coherent backaction 
effects can be ignored. One then finds the rule that 
for large asymmetries (|<?| 1), the SQD level position 

e = —V g effectively lies closer to resonance with the elec¬ 
trochemical potential of the more strongly coupled elec¬ 
trode. Now, for positive bias 14, as in Fig. Ha), the 
electrochemical potential fi s = +14/2 (fid = —14/2) of 
the source (drain) is resonant with the QD level for neg¬ 
ative (positive) values of V g . Thus, when the coupling 
to the source is larger (smaller) than that to the drain, 
the maximal transition factor r is achieved for negative 
(positive) V g as Fig. H(a) demonstrates. 




V„[T\ 


FIG. 6: Relative change K in the transition factor r, Eq. (EU, 
by the coherent backaction as a function of (a) the gate voltage 
V g for Vi, = 2.5 T and (b) of the bias voltage 14 for V g = 1.5 T. 
The different curves correspond to the indicated values of the 
tunneling asymmetry g = (r s — Fd)/(r s +Fd) for fixed average 
tunneling rate F = (r s + Td)/2 = 0.02T and bandwidth W = 
1000T. 
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By contrast, Fig. 0b) shows that the impact of the asym¬ 
metry g on K as function of the bias voltage is more 
complicated close to resonance. In general, K decreases 
with 14, as Fig. 0b) already showed for g = 0. In the 
limit 14 —> 00 , one expects n ~ 1/14 —> 0 and there¬ 
fore K —> 0. Moreover, we find K > 0 for small bias 
[generally valid] and for large bias [specific to the param¬ 
eters chosen in Fig. 0b)]. Thus, both for large and small 
bias the coherent backaction suppresses the measurement 
backaction. However, for intermediate bias voltages, K 
shows strong drops as Fig. 0b) illustrates. This sup¬ 
pression appears since k and (pl t / 2 — p® t ) change their 
sign in the vicinity to resonance at e = p r and there¬ 
fore become both very small. This suppresses K and for 
a small intermediate bias regime K can become negative 
as pointed out already in Sec. lIIID~2l Figure0b) reveals 
that the position and even the existence of these drops 
and sign changes depends crucially on the asymmetry g 
and the gate and bias voltage polarity. In particular, 
for V g > 0, as assumed in Fig. 0b), the source cannot 
become resonant with the level for positive bias. Accord¬ 
ingly, the drop is absent for g = 0.5 in Fig. 0 b), i.e., 
when the source is coupled more strongly to the SQD 
than the drain. 

In summary, the coherent backaction starts to can¬ 
cel the cotunneling-induced stochastic backaction on the 
flanks of the SQD resonances where one enters Coulomb 
blockade - the optimal spot for sensing - with correc¬ 
tions depending on the junction asymmetry, see Fig. 0 
Junction asymmetries commonly found in experiments 
strongly affect the relative importance of the coherent 
backaction terms and thus also the qubit decay rates. 


E. Effective evolution of quasistationary modes 

Since the qubit isospin is contained in the quasista¬ 
tionary degrees of freedom [see Eq. (l55l) ] , it is of interest 
to find a description only for their evolution. Impor¬ 
tantly, we do this without making further approxima¬ 
tions, i.e., we reproduce exactly the result for X q (t) as 
one does when solving Eq. (l26l) for p. For all practical 
purposes, one should solve the latter equation [or equiva¬ 
lently Eq. (00)]: Eqs. (1711) and (1721) below merely serve to 
highlight some properties of this solution, in particular, 
its non-Markovian nature and the impracticability for re¬ 
ally getting rid of the sensor dynamics. It also serves as 
a basis for Sec. EY1 

Equation 02 is a coupled set of differential equa¬ 
tions for X q (t) and X d (t). The equations can be solved 
most conveniently by Laplace transformation, defined 
by f(z) = / 0 °° dte lzt f(t) for any (well-behaved) time- 
dependent function /(f). This yields 

( 68 ) 


with a Laplace-frequency dependent effective Liouvillian 

aw = n'+m 

and matrices given by Eqs. (1571) (1521) . (RJT1) . (02), 

and (IB.341) . The frequency-dependent initial condition 
is given by 

X’„(z,0) = X«(0) + A-—(70) 

where we take the initial time to be f = 0. The inter¬ 
mediate steps of the derivation of Eqs. (l 68 l) (1701) are 

given in Appendix [Bj There we show how to reproduce 
the above result in the general framework of a projec¬ 
tion appro ach following that of Nakajima and Zwanzig 
[48l l55l - [57| . The projection approach simply separates 
the dynamics in the complementary subspaces spanned 
by the quasistationary and decaying modes. Importantly, 
this projection technique treats all the different types of 
backaction we discussed so far on the same footin g, w hich 
allows us to go beyond the approaches of Refs. ( 88 | and 
& 

Transforming the solution (1681) back to time space and 
exploiting the convolution theorem yields [see Appendix 

0 


X q (t) = dt'U q e q (t-t')X q eS (t'), (71) 

-'O 

x e yt') = x q (0)6(t' - 0) + H q e d (t')X d (0), (72) 

where 5(t' — 0) indicates a (5-function with infinitesimal 
shift and 


/ OO 

-CO 

/ CO 

-CO 


Z2L p -izt ___ 

2^ e z-Ll q s ( Z y 

— e~ izt A qd _ - _ 

27t z — L dd — A dd 


(73) 

(74) 


Equations ED and (1771) are the third set of main equa¬ 
tions. We stress that they are generally valid in the 
sense that they do not involve any approximations be¬ 
yond those needed for the validity of the kinetic equa¬ 
tions (l26l) : Eq. (17U) exactly reproduces X q (t) as obtained 
from the solution of Eq. 02. The expression for X q (t) 
simplifies drastically for times 1 > I /7 if one accounts 
only for the leading-order contributions in A/ 7 . Such an 
expansion is valid only in the weak-measurement limit 
A /7 and complies with our kinetic equations in in the 
high-temperature limit T jT -C A/T as we discuss fur¬ 
ther below in Sec. 13 Before doing so, let us first note a 
couple of general properties of the above solution. 

Non-Markovian dynamics and time scales. The 
effective Liouvillian (1671) confirms the discussion in 
Sec. IIII C 21 the unperturbed evolution in the quasista¬ 
tionary subspace is perturbed directly by A qq , which can 
be absorbed into a redefinition Lq 9 = L qq + A qq that just 
leads to the mean-field tilting of the qubit axis. 
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Moreover, the third term in Eq. (1691) gives an explicit 
expression for the indirect and in general non-Markovian 
perturbation of the quasistationary evolution by virtual 
transitions via the decay modes. This interpretation fol¬ 
lows most clearly from an alternative derivation in time 
space given in Appendix IB 21 and illustrates that the time 
delay between these transitions — giving rise to non- 
Markovian effects — converts into the frequency depen¬ 
dence of this expression. This term entails the effect of 
the stochastically fluctuating deviations from the mean 
field (see Sec. Ill Cl) and separates it from the mean-field 
effect contained in A". Since the denominator in this ex¬ 
pression does not grow exponentially (it rather tends to 
be constant because \L^ q \ ~ y), the voltage dependence 
of this term is largely determined by that of A dq , which 
we discussed above. 

The third term in Eq. (l69l) has, in general, several ef¬ 
fects: the eigenvalues of Lq Q in general acquire (i) a real 
part leading to an additional shift of the qubit frequency, 
and (ii) an imaginary part, which corresponds to the re¬ 
laxation and dephasing rates of the qubit. To extract 
both effects, one has to inverse the Laplace-transformed 
function (1551) . which leads to an integral that should be 
computed by the applying the residue theorem. The 
residues are determined by the zeros of the denominator 
satisfying z p = L q ^{z p ) and determine the time scales of 
the qubit evolution, (iii) Moreover, the third term may 
not commute with L^ q . This induces transitions between 
the unperturbed eigenstates and leads to a rotation of 
the qubit eigenbasis. (iv) Finally, since the third term 
is not Hermitian due to the decoherence that it induces 
(in contrast to L g 9 ), the qubit eigenaxes may not be mu¬ 
tually orthogonal any more. This renders the circular 
precession induced by L^ q slightly elliptical as we illus¬ 
trate below in Sec. m 

The effective Liouvillian (l69l) also allows for a compar¬ 
ison with earlier results. We note that our approach is 
conceptually quite similar to that discussed in Ref. 8 Sj. 
However, Ref. [ 88 j employs the clearly stated additional 
assumption that the effect of the electrodes is only to 
modify the SQD dynamics (assumed to obey a Lindblad 
equation) without affecting the (effective) coupling to the 
qubit. In our formulation, this would mean that the cou¬ 
pling to between the quasistationary and decay modes 
is mediated only by the interaction Hi = hA ■ t/2, i.e., 
only the stochastic backaction. This means that in the 
approach of Ref. [gf| both the dissipative and coher¬ 
ent backaction are neglected. Within this approximation 
one can thus not calculate the experimentally measur¬ 
able signal current; it must be consistently set to zero. 
(The signal current was not calculated nor of interest in 
Ref. (gg|). Moreover, Ref. [ 88 ] focuses only on the time 
scales and does not consider initial-slip effects that we 
next turn to. 

Initial slip. Equations d) and m show that the 
decay modes affect the quasistationary modes not only 
through the effective Liouvillian but also through the 
effective initial state [79j, |90|-|93[ : the latter is not just 


given by the initial quasistationary variables contained 
in X 9 (0). In addition, there is a term in (1721) that ac¬ 
counts for a initial contribution from the decaying sub¬ 
space, X d (0), followed by a time integral over transitions 
into the quasistationary subspace. This leads to an ini¬ 
tial slip that affects the quasistationary modes. Like the 
time evolution m the slip (1721) of the initial state has 
a a time-nonlocal expression in terms of the initial state 
of the decay mode X d (0) [in the Laplace transform (l70l) 
the corresponding second slip term has frequency depen¬ 
dence] . 

An experimentally relevant question is how to elimi¬ 
nate or minimize the initial slip since it can induce errors 
even for a perfectly prepared initial qubit state. In Ap¬ 
pendix^ we show that initial qubit-SQD states p(0) that 
exhibit no initial slip form a subset of measure zero in the 
set of all valid initial density operators p(0). In general, 
a sufficient condition for zero slip is that X rf (0) = 0, i.e., 
using Eq. (l55l) : 

MO) = PstP°(°) -PstP 1 ^) = 0, (75) 

MO) = Pst T °(0)-Pst rl (0) = ( 76 ) 

To find all initial states with zero slip, one has to compute 
the kernel of matrix A qd (z — L dd — A dd )~ 1 . Assuming the 
frequency z does not hit a pole of the denominator (for 
example when considering the Markov approximation for 
z = 0 ), the inverse (z — L^ d — A dd ) _1 exists and has full 
rank. Thus, to determine the dimension of the kernel, it 
suffices to determine the rank of A qd , which is 2 since 



for any x,y £ IR and using Eq. (l 68 l) . Thus, the matrix 
A qd (z — Ly d — A dd ) _1 has a kernel of dimension 2 , which 
means that the set of initial states has zero measure since 
X d can be taken from a four-dimensional set. 

One way to eliminate the initial slip is to switch off 
the capacitive interaction before t = 0 , i.e., A (t) = 0 for 
— l/y -C t < 0. Then the detector can establish a sta¬ 
tionary state and the initial SQD-qubit state factorizes: 
X d (0) = 0 [see the A = 0 solution (gS] and see also Ap¬ 
pendix E]). By contrast, if one switches off the current 
through the sensor, T —>• 0, before t = 0, then one starts 
with a sensor in a definite charge state p 1 ( 0 ) = 0 or 1 , 
which is highly nonstationary for typical operation pa¬ 
rameters of the SQD (tuned close to resonance for high 
sensitivity, one finds usually ~ pit)- The backaction- 
induced initial slip thus leads to an essential difference 
between two ways of switching off a sensor, which should 
be taken into account in designing detection protocols. 

The magnitude of the slip in general depends on fre¬ 
quencies only for z > L^ d + A dd ~ 7. This means that 
the initial value X d (0) influences X 9 (t) for times t ~ I/7 
[through the integral (l72|) ]. In Sec. II V Cl we investigate 
the slip magnitude in more detail in a high-temperature 
limit to leading order in A/7; importantly, we find that 
even in this simple case there is a slip effect of order A/7 
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which affects the overall qubit dynamics by, e.g., phase- 
shifting the solution in a way depending on the sensor 
initial state. 

From the above we can generally conclude within the 
regime of validity of the kinetic equation (1761) that due to 
the backaction-induced initial slip additional errors are 
generated; since the zero-slip states are sets of measure 
zero, it is clear that preparation errors of the qubit-sensor 
state will invariably lead to an initial slip. Models of such 
possible errors are actually relevant for a different branch 
of quantum information. In quantum-error correction, 
one deals with decohcrence from the environment in a 
phenomenological way by introducing additional bit-flip 
errors [95}. This requires assumptions to be made about 
the type and statistics of the different possible errors. 
Our work thus provides in this context a possible scenario 
how such errors may arise and how they can be modeled, 
after, e.g., a measurement has been performed. This will 
become more concrete in the next section, where we dis¬ 
cuss simplified equations, which have, however, only a 
limited range of applicability. 


IV. HIGH-TEMPERATURE QUBIT DYNAMICS 

The kinetic equation, either in representation (1261) or 
(156|) . form central results of this paper. They fully suf¬ 
fice to compute the transient dynamics of the charge- 
specific isospins T°(f) and T 1 (t). From this result, the 
total isospin r(f) = r°(f) +x 1 (t), i.e., the reduced qubit 
state can be constructed. However, there are several rea¬ 
sons to attempt to obtain a closed description in terms 
of r(t) only. 

First, from Eq. m, it is not directly clear on which 
time scales r(t) evolves or decays. Second, an effective 
qubit description plays an important role, for example, 
in quantum error correction. It is an interesting question 
how far Eq. (17B1) can actually be reduced to a closed 
equation for the reduced density operator for the qubit 
alone. One indication that this requires additional as¬ 
sumptions is that © and (1561) can be solved only if 
the full initial state vector [either in form Eq. (fl3l) or 
Eq. ([55]) ] is specified and not just the sum t°(0) +t 1 (0). 
We have already seen that the initial values of the other 
degrees of freedom produce an initial slip [see Eq. m] 
and will see below that even in lowest nonvanishing order 
this slip cannot be avoided. The third consideration is 
related to this and concerns the minimization of backac- 
tion in quantum-information processing: One would like 
to know the effective qubit eigenmodes, e.g., to construct 
initial states that are least sensitive to backaction by set¬ 
ting experimental parameters. 

To investigate all this further systematically, an effec¬ 
tive theory for the qubit evolution in a simple limit is 
useful. In this section, we consider the regime where 
the coherent backaction and 0(T/T) corrections to the 
stochastic backaction can both be neglected. We stress 
that this is for the purpose of illustration mostly since 


the latter effects are relevant under typical experimen¬ 
tal conditions. This simplification allows us to perform a 
expansion of the effective Liouvillian (1691) and the initial 
slip (171)1) to leading order in A /7 and we investigate the 
resulting transient qubit evolution here in some detail. 
In Sec. IIIID II we give tangible analytical expressions for 
the relaxation and dephasing rates as well as for the qubit 
precession frequency. In Sec. IIVB1 we assess the accu¬ 
racy of the approximate theory by comparing with the 
numerical solution of the full kinetic equation (1261) . We 
further discuss the slip of the initial condition of the qubit 
isospin in Sec. IIV Cl which relates to a “kick” that the 
qubit experiences during the relaxation time of the sensor 
QD. Finally, we show in Sec. IIV Dl that the measurement 
backaction forces the isospin to precess about a tilted axis 
in an elliptical way. The eccentricity is connected to oscil¬ 
lations in the decay of the purity of the qubit state. This 
illustrates concretely that our density-operator approach 
goes beyond standard master-equation approaches as we 
discuss in Sec. m 

A. Effective Liouvillian, initial slip, and mode 
vectors 

The kinetic equation [Eq. (1761) or Eq. (1561) ] was de¬ 
rived using the weak-coupling, weak-measurement limit, 
T/T -C A/T -C 1 where A ~ A, fh This prevents one 
from just expanding in A/T since that would imply tak¬ 
ing A/T -C T/T. However, if we consistently neglect 
the corrections T /T (the cotunneling corrections to the 
stochastic backaction as well as dissipative and coherent 
backaction), then we can take A/T —> 0 and expand in 
this parameter. We refer to this as the high-temperature 
limit (since it is the large temperature that allows one to 
take the infinitely weak-measurement limit). It should be 
noted that in this approximation the current through the 
sensor QD is zero, i.e., at this level of the theory one is 
not accounting for the actual backaction effects due to the 
current measurement (rates ~ TX/T), but only for the 
leading effect of the tunnel coupling. Below, the station¬ 
ary occupations p 0,1 = 7 1,0 /7 are given by their leading 
order expressions (SET rates) 7 0 ' 1 = J2 r =s d^ r fr [ see 
Eq. (1771) ] and 7 = 27 ° + 7 1 . 

We thus simplify the isospin evolution obtained from 
the effective Liouvillian (1M1) and for concreteness assume 
from hereon 

Q = Tle x , (78) 

perpendicular to the capacitive interaction vector A = 
Xe z . This means that if we ignored the mean-field tilting 
$7 —> £7 (which we do not), the qubit would oscillate in 
the measurement basis. 

As we explain in detail in Appendix [B] the isospin evo¬ 
lution contained in Eqs. (l69l) and © can be simplified 
by performing a Markov approximation, i.e., by replac¬ 
ing z = 0. This Markov approximation with respect to 
memory induced by the sensor QD (after integrating out 
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the electrodes) is valid in the weak-measurement limit 
A/r < 1 (see Appendix IB 3 cl) . We note that in the 
high-temperature limit also non-Markovian corrections 
due to tunneling processes (r) are consistently neglected 
as next-to leading order T 2 /T corrections are not ac¬ 
counted for. With z = 0, the Laplace-transform inverse 
of Eq. (| 68 |) can be easily performed. Expanding the de¬ 
nominator in powers of A/r and extracting the isospin 
from X 9 , we find 


general case, Eq. (1721) . If the measurement is not weak 
any more, i.e., if A ~ 7 , one can expect the dynamics of 
the charge-specific isospins to become important on the 
entire time scale of the qubit decay and not just through 
an effective slip of the initial condition. Still, even in the 
weak-measurement limit studied here, Eq. (ED shows 
explicitly that the slip is of non-negligible order A/ 7 . 

2. Qubit time scales 


T {t) = e~ iLeltt r e ff(0) + 0 (A 2 / 7 2 ). (79) 

In this approximation the stationary state r(oo) = 0. 
The effective Liouvillian in Eq. © reads in diagonal 
form: 

— / ] iiSl a 7»)wff,ad e fj a , ($0) 

a= 0 ,± 


The simple formulas (17(71) - (15T1) for the relevant preces¬ 
sion, relaxation, and dephasing time scales form the third 
central set of equations of the paper. We now discuss 
their contents. The eigenvalues in Eq. (15(71) contain the 
effective qubit frequencies 

Sl a = att, a = 0, ±, (82) 


where the eigenvalues, the left and right eigenvectors, 
and the effective initial state T e a(0) are specified below 
by Eqs. ED ED- Equation m is valid for times 
I /7 <C t <C 7 2 /A 3 as Appendix [5] shows and also 
our numerical checks below confirm. The l ower limit 
indicates that we consider the wide-band limit |l05l ] with 
respect to the sensor QD band-width 7 (by setting 
7 > z -> iO above), whereas for times t 7 /A 2 , 

corrections of 0 (A 3 / 7 2 ) to the effective Liouvillian (l80l) 
accumulate and the error made for r(f) may become 
sizable. 


which up to 0(A 3 /y 2 ) read 


si 


si 


+ (rf t A) 2 . 


(83) 


Equation (l83l) is precisely the length of the mean isospin 
field announced earlier in Eq. CD. but now with ( n ) = 
Pst* 

$7 = SI + pit A = S7e||, (84) 

whose unit direction vector is relevant for the following: 


1. Initial slip 


= sisi 


Pst AA ) A 


(85) 


The effective initial condition appearing in Eq. CD 
reads (see Appendix 151): 

Teff(O) = t( 0) — — A X (Pg t T°(0) — Pg t r 1 (0)). (81) 
7 

This shows that even in this simple limit the qubit 
description is still not closed. Although the entire 
sensor variables (electrodes plus sensor QD) have been 
eliminated, the initial condition does not depend only 
on t( 0) (see Appendix [C]). Instead, it additionally 
requires the specification of the component of the 
initial qubit-sensor state p( 0 ) in the decaying subspace. 
Thus, both initial charge-specific isospins r"(0) are 
needed to compute r(i). In contrast to the general 
case discussed in Sec. IIII El Ea. ED only relies on 
<$st( 0 ) = 74 t °( 0 ) — Pst rl (°) and does not involve 

<5 st (0) = pltP°(0) — PstP 1 ^)- The reason is that the 
denominator in Eq. d is approximated here by —7 
and A qd does not act on the charge sector [see Eq. ©]. 
As mentioned in the general discussion of the initial slip, 
Eq. CD. initial qubit-sensor states p(0) with zero initial 
slip form a zero-measure subset of all possible initial 
states. As a result, the initial slip adds to preparation 
errors, an error that depends on the sensor dynamical 
state. In the present simple limit, the initial slip is 
a time-local expression, Eq. ED, in contrast to the 


Moreover, we will need the perpendicular unit vector, 

e_L = — pl t \SlJ/SI, ( 86 ) 

lying in the plane spanned by SI and A. Our calculation 
thus confirms the intuitive picture explained in Sec. Ill Cl 
The mean-field effect of the average SQD charge (n) = p\ t 
is just to tilt the qubit axis (see also below in Sec. IIVB II) 
and does not rely on the tunneling-induced fluctuations 
~ Sh , see Eq. ED- Tunneling influences the mean sen¬ 
sor charge only indirectly as noted in the discussion of 

Eq. ED- 

Compared to the first term of Eq. ED- the eigenvalues 
of the quasistationary modes have acquired small dissi¬ 
pative parts [see Eq. ED]: 

7o = 7± = (87) 

11 12 

Here, the relaxation rate is given up to 0 (A 4 / 7 3 ) by 


1 (A • ej_) 2 /DA 2 A 2 

T[ ~ T 7 - r UJ 7’ 


( 88 ) 


which is quadratic in the component of the measurement 
vector A perpendicular to the average isospin field SI with 
the transition factor r given up to zeroth order in T /T as 
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given by Eq. (l64l) . The dephasing rate 1/T 2 is expressed 
compactly in terms of the pure dephasing rate 1 /T^ = 
1/T 2 — l/(2Ti) (23 U P t° order 0( A 4 / 7 3 ) as 


1 

fl 


(A-e„) S 


x Plt 

n 


A 2 

7 


(89) 


which is quadratic in the projection of A on the unit 
vector ey = fl/fl along the average isospin field. In the 
following, we refer to both relaxation and dephasing as 
decoherence because both drive the qubit into a mixed 
state. 

Note that in both decay rates, the transition factor 
r appears, which links to the discussion of the previous 
sections. If higher-order A /7 terms are included into 
the relaxation and dephasing rate, additional terms ap¬ 
pear that depend on the dissipative backaction terms ~ c 
stemming from the off-diagonal elements in Eq. (Rail) . 

It is easy to see that the rel axati on and pure dephasing 
time Tj and T/, are positive [l06j | , since the transition 
factor is given to lowest order by r « 2r + r_/(2T + + 
r_) 2 > 0 [see Eq. (R>H) ]. Since > 0, the ratio T 2 /2Ti 
further satisfies the relation {27j[: 


T 2 1 

m = 1 + 2(pg t A/n) 2 < L (90) 

Equations (1551) and (1551) again confirm our intuitive ex¬ 
pectation from Sec. Ill Cl that only the fluctuating part 
of the SQD charge ~ A Sn (involving here virtual transi¬ 
tions into the decay modes) is responsible for the qubit 
decoherence. The energy scale A 2 /7 for the decoherence 
rates exhibits the expected quadratic scaling with the 
weak coupling A and inverse scaling with the large detec¬ 
tor band width 7, as discussed further below in Sec. IV Al 
The decay time is thus slow compared to the time scale 
of the relaxation ~ I /7 of the sensor QD and that of the 
intrinsic evolution of the SQD-qubit system ~ 1/A. We 
emphasize that this picture and our approach hold only 
in the limit A <C T: it breaks down if the tunneling be¬ 
comes to strong relative to eith er the measurement A or 
the qubit internal field f 1. [107] 

One should note that the corrections to the decoher¬ 
ence rates y a [see Eqs. (1551) and (1551) ] are of 0 (A 4 /7 3 ), 
while the corrections to the qubit frequency (l82l) already 
appear in lower 0(A 3 /'y 2 ). The reason is that these 
quantities behave differently under a simultaneous re¬ 
versal of the orientation of A and fl as one can see 
from a simple physical argument: mapping A —> — A 
and fl —> — $7 corresponds to spatially mirroring the de¬ 
tection setup about the vertical axis in Fig. [1] This 
clearly inverts the sense of the precessional motion, i.e., 
one has fl a (—A, — $1) = —fl a (A, fl) but the qubit de¬ 
cay cannot depend on mirroring the setup, i.e., we have 
7 a (—A, — fl) = 7 a (A, f2). This implies that corrections 
to 7 a must be of even order in A and therefore at least 
of fourth order in A, in agreement with our calculation. 

3. Mode vectors 


To complete the specification of the effective Liouvil- 


lian (1551) , we give the explicit formulas for the unit vectors 
e e ff jQ . Expressed in ey and ej_ given by Eq. (1551) and (1551) . 
respectively, they read: 
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(92) 

(93) 


a 2 fi 

+r^^ — 
fi 2 



fi 2 


-n 2 + nn/2 

-=-ej_ 

fi + fi 


where oc indicates that we suppress the normalization 
constants. As before [see Eqs. (1551) and (1551) ]. we define 
e eff,± = (e e ff,i T *e e ff j 2 ) /V2 and we note that ey x ej_ = 
f2 x A [see Eqs. (1551) and (1551) ]. 

In stark contrast to the unperturbed case, the real unit 
vectors e^o, e e ff.ij and e e ff j2 form a real nonorthogonal 
basis. This implies that to decompose a vector in the 
basis {e e fj iQ ,}, one needs to take the scalar product with 
the dual basis denoted {e e ff iQ ,}, see Fig. [7] The dual basis 
vectors e c ff jC[ are non-unit vectors orthogonal to the plane 
spanned by e e g^ and e e g 7 , 


®eff ,o: 


X 0eff,7 

|^eff,o: * (^eff,/3 X 6eff,7)| 


(94) 


where (a,/ 3 , 7 ) is a cyclic permutation of (0,1,2). We 
refer to this nonorthogonality of the eigenvectors in the 
following simply as a distortion of the isospin modes. 



FIG. 7: Distortion of isospin mode vectors: (a) Three- 
dimensional sketch of the effective qubit axes e e fj jQ! [a = 0 , ±, 
Eqs. dll]) ©] and its dual basis e e ff, a . Note that e e ff,o has 
a nonzero projection on the plane spanned by e e ff,i and e e ff,2, 
while e e ff : o is orthogonal to this plane. Likewise, a vector in 
the e e ff,i-e e ff,2 plane has a nonzero component along e e ff,o- 


The distortion of the mode vectors scales with the small 
ratio of the magnitude of the decoherence rates A 2 /y rel¬ 
ative to the effective qubit frequency fl ~ A and is fur¬ 
thermore suppressed by the transition factor r when go¬ 
ing off-resonance [see Eq. (1551) ]. Since we allow for A ~ fl, 
the remaining factors in Eqs. m ® can be of order 1 . 
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This distortion also has tangible physical consequences: 
Eq. (15T1) shows that the precessional motion of the qubit 
isospin t becomes tilted. It is also not circular any more, 
but becomes slightly elliptical instead as we investigate 
in more detail in Sec. Evd] 

f. Other approaches 

At this point, it is instructive to compare with some 
other approaches. At first sight, the mode distortion 
may appear peculiar and one may wonder why it does 
not show up in other approaches. In fact, the mode 
distortion disappears in the limit A/ST —> 0, i.e., when the 
capacitive coupling becomes smaller than all other en¬ 
ergy scales. In this limit, the modes become orthogonal 
with eo = n. This diminishes also the mean-field effect. 
This limit is equivalent to the frequently made secular 
approximation that Davies has shown to be exact in this 
strict weak-coupling limit 5^, [59j. Here, however, we 
consider the more general situation that the capacitive 
coupling A can be of the same magnitude as the internal 
qubit energy scale fl. The secular approximation is not 
applicable in this case as we showed previously [32] in 
accordance with other works [46]. This is furthermore 
signalled by the observation that the secular approx¬ 
imation conflicts with the isospin conservation when 
electrons tunnel between the electrodes and the sensor 
QD. The strength of the mode distortion thus reflects 
the importance of nonsecular corrections (coherences). 
We next compa re t o Ref. [33], whose approach is 
similar to ours [108 ] . This study applies a two-step 
procedure to derive a closed, effective description of 
the qubit dynamics, starting from a generalized master 
equation for qubit plus SQD, as we do. To treat 
the limit Cl A 2 / 7 , it is additionally assumed that 
|A ■ ejJ -C | A - e|| |, i.e., perpendicular fluctuations of 
B c fj along ej_ are small as compared to longitudinal 
fluctuations along the mean field f2. Perpendicular 
fluctuations are therefore treated there perturba- 
tively. In this way coherences between isospin states 
quantized along ey are only effectively included. By 
contrast, we allow for comparable fluctuations in the 
direction of both ey and ej_, which is a more general case. 

B. Accuracy of effective isospin dynamics 

Before we illustrate the effects of the initial slip and 
the mode-vector distortion, we discuss the accuracy of 
our Liouville-space perturbation theory as compared to 
the full solution of the kinetic equations (l26l) in the high- 
temperature limit and up to times times t~ 7 /A 2 ^>l/ 7 - 

1. Undistorted, mean-field qubit modes 

As a starting point for this discussion, we first discuss 
our result in view of the rough mean-field picture of the 


detector backaction. It is valid in zeroth order in A and 
B relative to T, where we neglect all relaxation and deco¬ 
herence rates A 2 /7 -C A. Expanding the expressions (l9ll) 
(1551) for the mode vectors to the corresponding zeroth 
order yields an orthonormal basis, 

^eff,o ~ ~ j_ > and o e ff 2 ~ ^|| x e j_, (95) 

with ey and e_L given by (1551) and (1551) . respectively. 
Thus, in this approximation the dual basis {e e fj a } coin¬ 
cides with {e e ff iQ ,}. In analogy to the unperturbed case 
discussed in Sec. m C 11 the isospin just precesses cir¬ 
cularly about ey ; however, the precession frequency - H 
instead of B - and the precession axis - along ey = Cl /Cl 
instead of fl are different. This rough mean-field pic¬ 
ture of the measurement is therefore to tilt the “bare” 
isospin field to the mean isospin field (1551) by an angle 
[see Eq. (12TT) ]: 

tan$ = (96) 

In our concrete charge-qubit model, this means that the 
capacitive readout simply detunes the charge qubit due 
the gating effect of the sensor QD with mean charge 
Pg t . Here the mean charge is identified with the ensem¬ 
ble averaged charge, see the related discussion in Sec. 
Ill Cl Since we only require A < 7 but impose no con¬ 
straint on the ratio A/D, this angle can be large. In Fig. 
[51 we illustrate this effect by showing the evolution of 
the three isospin components in the basis ( e x ,e y ,e ~) = 
(D,A x D,A) on a long time scale t I /7 when the 
isospin initially points into the direction of Cl, i.e., per¬ 
pendicular to A. If the coupling A was switched off, 
we would expect the isospin not to precess at all and 
to remain stable along D. By contrast, the oscillations 
of all components in Fig. [5] clearly demonstrates that 
the isospin revolves about a very different axis, roughly 
pointing in the direction of fl + A in line with Eq. (1551) 
for the parameters employed here. 


2. Accuracy of weak-measurement expansion 

The green curves in Fig. [5] depict the difference be¬ 
tween the isospin evolution obtained by solving the ki¬ 
netic equations (USD and the evolution computed from 
Eq. cna based on our perturbation theory, indicat¬ 
ing that both agree well (plotted is the error multiplied 
by 100). The remaining deviation is mostly due to a 
small phase shift between the full and perturbative so¬ 
lution that accumulates in time. The origin lies in ig¬ 
nored corrections of order of ~ A 3 /y 2 to the effective 
qubit frequency (1521) . In Fig. [51 this accumulates after 
the shown time At ~ 25 • (27 r/Cl) to a phase difference 
A ip < 25 • (A/4f ) 2 ~ 0.01, which is just visible. How¬ 
ever, we emphasize that the accuracy of the decay rates 
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FIG. 8: Time evolution of the isospin vector components (a) 
t x , (b) Ty, and (c) r z for times t > I/7. The blue curves show 
the analytic solution -r ana given by Eq. EH, whereas the green 
curves show 100 times the error with the r obtained from 
the solution of the full kinetic equations (1261) . The coordinate 
system is chosen as (e x ,e y ,e z ) = (fi,A x Cl, A). The initial 
state of sensor plus qubit is p°(0 ) = 0, t°(0) = 0, p 1 (0) = 1, 
and r 1 (0) = S2. The remaining parameters are T s = Vd = 
f = 10“ 3 T, A = fi = O.lf = 10“ 4 T, Vb = V g = 2.5 T, and 
W = 1000T. For these parameters, we find 7 ~ 3.75 • 10 _3 T 
and Cl ~ 1.37 • 10 _4 T and therefore I/7 rs (27r/f2)/200 -C 
t -C 7 2 /A 3 « 300(27r/f2) is well fulfilled for the times shown 
above. 

is higher as discussed below Eq. (1621) and therefore the 
exponentially decaying envelope of the isospin evolution 
agrees with much larger accuracy up to longer times. 

The rough mean-field picture introduced above com¬ 
plies with the physical picture developed in prior works 
1HH28], HI]. It also forms the basis of a simple classi¬ 
cal understanding of the qubit decoherence in fluctuator 
models, to which we compare our results in Sec. IV Al 
However, there are important corrections to this picture 
even in the high-temperature limit, which we next dis¬ 
cuss. 


C. Effect of initial slip 

A first illustration of the corrections to the mean- 
field picture is the effect of the slippage of the initial 
condition, Eq. m- To illustrate this effect, we com¬ 
pare in Fig. [9] the solution for the isospin for two dif¬ 
ferent initial states. We start from a factorizable ini¬ 
tial state pqs = Pq ® Ps with a fixed total isospin 
t(0) = Cl along the “bare” internal field [determining 
Pq = (1q + r(0) • t) /2] while changing the initial con¬ 
dition for the SQD through the sensor charge equal to 
p 4 (0) [determining ps = (1 — p 1 (0))P° +p 1 (0)P 1 ]. We 
show the outcome for r ana for the two cases of an ini¬ 
tially empty SQD (p 4 (0) = 0) and a SQD hosting an 
electron (p 1 (0) = 1). Figure [9](a) exhibits a phase shift 
between the two isospin evolutions that persists over an 
entire qubit cycle (and in fact for all future times, which 
are not shown here). 


The approximate analytical solution (1791) is valid over 
the entire time scale shown in Fig. [HXa) except for 
very small times t, < I/ 7 . In Fig. [Sib), we illustrate 
how the isospin computed from the full kinetic equations 
(l26l) (red) approaches the approximate analytical solu¬ 
tion (17(71) (blue): all curves for the full solution (red) 
start from the same value for t v ( 0 ) = 0 but immedi¬ 
ately develop differently depending on the initial SQD 
charge p 4 (0). On a time scale ~ I/ 7 , they approach the 
approximate analytical solutions r ana ( 0 ) (blue), which 
are offset by the initial slip (15TT) . Figure [9][b) confirms 
that precisely due to this slip the analytic solution ac¬ 
curately approximates the full numerical one for times 
t I/ 7 . The latter time scale is expected since the 
approximate curve relies on the SQD wide-band limit. 
(The analytic solution may even be unphysical for times 
t < I/ 7 : In certain cases, including Fig. [9l one may find 
l T off(0)| > 1.) 

If one, however, neglects the initial slip (1811) one ob¬ 
tains a curve similar to the blue one in Fig. [9] but with 
zero vertical offset for time t = 0. This clearly leads to 
a nonnegiigible deviation from the full solution for ini¬ 
tial conditions with nonstationary sensor. After a time 
t ~ I/ 7 , the qubit phase is advanced by D /7 which can 
be of the same order as the phase angle ~ A /7 of the 
initial slip (1811) (depending on the relation A, Cl, within 
the restriction A, Cl -C T). We stress that this slip leads 
to a cumulative effect: Even at long times f > I /7 the 
approximate solution without slip remains offset relative 
to the full solution. 

Altogether, this shows clearly that one cannot get rid 
of the detector completely - even though we describe 
the qubit state using only the Bloch vector r(f). It is 
difficult to eliminate the slip by a choice of the initial 
state of the coupled qubit-sensor system as mentioned 
(see Sec. lIIIEl) and further discussed in Appendix [Cl Im¬ 
portantly, one should note that it may not be possible to 
remove the initial-slip backaction by any qubit prepara¬ 
tion (i.e., just its reduced density operator r). For quan¬ 
tum error correction, it is thus important to model the 
failure of the preparation not only of the qubit dynamical 
state, but also the dynamical state of the sensor QD and 
their mutual correlations. 


D. Distortion of isospin mode vectors 

The rough mean-field picture also breaks down when 
accounting for the backaction effect on the isospin modes 
vectors : the eigenvectors e e g a are modified from Eq. 
(l95|) to Eqs. (IfJTl) (l93|l when taking into account the finite 
decoherence rate A 2 /q. This leads to both a tilting of the 
qubit axis and elliptical isospin precession. 
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tionally tilts the relaxation mode vector from ey —> eo^ff 
through the virtual-transition terms oc A 2 /y. 
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FIG. 9: Effect of initial slip: time evolution of the isospin com¬ 
ponent r ana computed from the analytical expression GS»- 

(a) Component r ana for a full qubit cycle. We assume a 
factorizable initial state of SQD plus qubit factorizes by 
taking t™(0) = p n ( 0)f2 for n = 0,1. Both solutions show 
a phase shift with respect to each other that does not die 
out, i.e., it persists even over many qubit cycles. In inset 

(b) , we compare r™ a (blue) with the component r y com¬ 
puted from the full kinetic equations (1261) (red). The high- 
temperature approximation r™ a approaches the full solution 
T y on the time scale ~ I /7 during which the SQD approaches 
stationarity. The inset (c) shows r^ na after one qubit cy¬ 
cle for the two initial conditions and illustrates that the off¬ 
set of both curves has not changed appreciably as compared 
to the evolution close to t = 0. The coordinate system is 
chosen as ( e x ,e y ,e z ) = (Cl, A x Cl, A). The parameters are 

r s = r d = f = io _3 t, a = ci = o.if = kt 4 t, v b = 3t, 

V g =T,W = 1000T, resulting in 7 = 27 ° + 7 1 « 3.30 • 10" 3 T 
[see Eq. (12711 ) and Cl « 1.27 ■ 10 _4 T [see Eq. (l83l) [. 


FIG. 10: Distortion of isospin modes: (a) Comparison of the 
components of the high-temperature approximation (1791) for 
the isospin r(f) along the mean-field qubit axis, ey = Cl /Cl, 
(black), along the tilted relaxation mode vector e e fj : o (red), 
and along its dual e e fr,o, (blue) as a function of time. Plotted 
are In (ey ■ r) + t/Ti (black), In (e e fj,o ■ r) + t/Ti (red), and 
In ((e e fr,o/ |e e ff,o|) • t) + t/T\ (blue), respectively. The initial 
state is given by p 4 ( 0 ) = 1 — p°( 0 ) = 0 , r( 0 ) = t°(0 ) = S7 
and all other parameters as in Fig. [9] Note that the decay 
of the precessional component is hardly visible here because 
T 2 ss T\ for the parameters chosen here and because of the 
short time window shown in this figure (see Fig. [HJ) (b) Sketch 
of the orientation of ey, e e fr,i, and e e fr,i relative to the e e ff,i- 
e e ff ,2 plane. The thick green line indicates the precession of 
the isospin in a plane parallel e e fF,i-e e fE ,2 plane and shifted 
along e e ff,o- The precessing part of the total isospin leads to 
an oscillation along the components along e^o (as indicated) 
and also along e||. Only the component along e e ff,o does not 
oscillate. 


1. Tilting of precession axis 


As a first consequence, the effective precession axis 
eefj ’0 acquires an additional tilting beyond the mean-field 
effect. This manifests as a nonzero component of e e ff o 
along e|Xei = Ox A, which is perpendicular to both 
the intrinsic qubit precession axis Cl and the measure¬ 
ment vector A and therefore perpendicular to Cl. By 
virtue of Eq. (15T1) this rotates the qubit axis relative to 
ey by an angle 


X 


r A 2 ClXpj 

jci n 2 


(97) 


plus higher-order corrections. This tilt becomes notice¬ 
able close to resonance, where detection is performed, as 
we illustrate in Fig. □Ufa), where we plot the projection 
of r(t) onto the mean-held axis ey = Cl/Cl. In addition to 
an exponential decay with the relaxation rate 1 /T \, this 
component acquires an additional oscillatory component 
for a general initial state. This simply indicates that we 
are looking at the component along a vector that is not 
the zero eigenmode of the qubit: the backaction addi- 


However, the component of r(f) along the tilted relax¬ 
ation mode vector e c ff,o shows also an oscillation (super¬ 
posed on an exponentially decaying contribution) as the 
red curve in Fig. [UK a) illustrates. This is an effect of the 
mode distortion: the plane spanned by the unit vectors 
e e ff,i and e e ff .2 is not anymore orthogonal to e e ff,o- Since 
the precessing part of the isospin lies in this plane, the 
projection of the isospin r(t) on e e ff,o becomes oscillatory 
as sketched in Fig. fTOl b). These oscillations are there¬ 
fore damped with the dephasing rate I/T 2 , which differs 
from the relaxation rate 1/Tj that sets the time scale for 
the exponential decay of the nonoscillatory contribution. 
Thus, the nonorthogonality of the isospin mode vectors 
mixes relaxation and dephasing in a nontrivial way. 

The only component of r(t) that does not oscillate for 
an arbitrary initial state is the one along e^o as the blue 
curve in Fig. flTH a) shows. The reason is that the dual 
vector e e ff,o oc e e ff,i x e e ff ,2 is normal to the e^i - e e ff ,2 
plane [see definition ©]■ This normal component is 
also independent of the dephasing time: it simply decays 
exponentially with the relaxation rate 1/Tj. 

However, the dual vector e e gf o should not be con¬ 
fused with the zero mode If we initially have 
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r e ff( 0 ) oc e e ff o, then the isospin still has a precessional 
component since e^o is not the relaxation mode vector. 
If one aims to prepare the qubit in a state whose Bloch 
vector direction is stable under the time evolution that 
includes the measurement backaction, one should take 

r eff (0) = Fe 0 , e ff, (98) 

where F is suitable real constant. [Il09 ] For the initial 
state (l98l) . we indeed find pure exponential decay to the 
stationary state r(oo) = 0: 

r(t) = Fe~ t/Tl e eSfi . (99) 

The above illustrates that the notion of “exciting a 
qubit mode” has to be treated with care due to both 
backaction-induced initial slip and due to mode vector 
distortion. Due to the initial slip discussed above, one has 
to prepare the qubit-sensor state [see Eq. d72l) . Eq. (18T1) . 
and Appendix [C] very carefully in order to achieve the 
initial condition ( 1 M 1 ) . 


2. Elliptical precession 


The second qualitative consequence of the distortion of 
the qubit modes due to the finite decoherence rate con¬ 
cerns the precessional motion in the e e g ! - e e g 2 plane 
with normal e^o. The trajectory of the isospin vector 
r(t) in this plane is changed from a circle to an ellipse. 
We illustrate this in Fig. filial for an effective initial 
isospin T c g(0) = Fe^i lying in this plane (again for a 
suitable real constant F). Applying Eq. (ITUl) . the evolu¬ 
tion of the isospin can be expressed as: 


r(f) = e~ t ^ T2 F cos(f 2 f)e e ff,i + sin(r 2 t)e e ff , 2 (100) 

_ e -t/T 2 p i^l — e/2 cos It + ^ vi 

+ \Jl + e/2 sin (fit + ^ v 2 ( 101 ) 

In the rewritten form, the part in the bracket describes 
an elliptical motion with linear eccentricity 

' = 2e *»-‘' e *<« = W, = (102> 


which is maximal near the resonance of the sensor QD 
[due to r, see Eq. (l62l) ] and proportional to the scale A 2 /q 
for the decoherence rate relative to the qubit frequency 
f l. The precession plane is spanned by two orthonormal 
vectors, 


Vl,2 


A G c ff 2 

v / 2"±7 


(103) 


which are at the same time unit vectors along the prin¬ 
cipal axes of the precession ellipse as sketched in Fig. 

Mb). 




FIG. 11: Elliptical precession, (a) Magnitude of the high- 
temperature approximation of the isospin |r(t)| (green). We 
take the initial condition p 1 ( 0 ) = 1 —p°( 0 ) = 0 , r e fr( 0 ) = e e fr,i 
here such that r(f) evolves in the e e fr,i - e e ff ,2 plane. From 
Eqs. (HU and m, it is easy to see that e e fi,i has a much 
larger component along A than e e ff !2 , leading to a stronger 
decay initially (see explanation in the text). The oscillating 
deviations of the isospin magnitude from the exponential de¬ 
phasing, |T e ff ( 0 )|e - A T2 (blue), reveal the elliptical precession. 
All other parameters are as in Fig. [Ij] Although the effects 
are weak in our controlled perturbative calculations, they in¬ 
dicate qualitatively new features that can be expected to grow 
for stronger readout couplings, (b) Two-dimensional sketch 
of the ellipse with principal axis vi and v 2 [see Eq. (11031) 1 de¬ 
scribed by the tip of the isospin for initial condition as in (a). 
The exponential shrinking with rate 1/T 2 is not indicated for 
simplicity. 


Equation (I101D shows that the magnitudes of both prin¬ 
cipal axes shrink exponentially with rate 1/T 2 . The mag¬ 
nitude of the isospin [green curve in Fig. fili al] thus os¬ 
cillates around a pure exponential decay with that rate 
[blue curve in Fig. II 1 T ad] . The oscillations are a signature 
of the elliptical distortion of the precessional motion. 

If there were no mode distortion, the above picture 
would be true also for slipped initial states out of the 
epeff - e e ff : 2 plane. However, due to the mode distor¬ 
tion arbitrary initial states must again considered with 
care because the projection of r(t) onto the e e ff.i - e e g !2 
precession plane has an additional contribution arising 
from the component of t along e^o [see Fig. »)]: 
this causes the center of the ellipse to be shifted away 
from origin. This shift of the center decays exponentially 
towards the origin with the relaxation rate 1/Ti, which 
is again different from the decay rate of the precession, 
which is damped with the dephasing rate 1 /T 2 . 

Relation to state purity. The discussed elliptical 
isospin motion reflects the exponentially damped but os¬ 
cillatory decrease of the qubit-state purity due to the 
readout by the sensor QD. This can be roughly under¬ 
stood as follows. The dephasing is the strongest when 
the isospin and the measurement vector A are perpendic¬ 
ular to each other because this corresponds to the charge 
qubit electron being delocalized between the two sites. 
Here, small fluctuations in the detuning, induced by the 
stochastic switching of the SQD, then introduce a strong 
dephasing. By contrast, no dephasing appears when the 
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isospin and the measurement vector A are collinear, i.e., 
when the charge qubit electron is localized in one of the 
QDs. This interpretation is consonant with the situation 
in Fig. HHa), which shows the magnitude of the qubit 
Bloch vector as a function of time. The qubit Bloch vec¬ 
tor has initially a large overlap with A, leading to a sup¬ 
pressed decay in the first quarter of a precession period 
2ir/Q. The oscillatory reduction of the qubit purity was 
anticipated also in our discussion of the exact relation 

m, 

j t [\r(t)\ 2 } = -2A ■[r°(f)xr 1 (t)]. (104) 

If r(f) precesses, its charge-specific components T n (t ) 
also precess and their components along A ch ange in 
time, resulting in a nonexponential purity decay |lldj |. 

Finally, we note that these oscillations of the purity 
decay are not an effect of non-Markovian corrections [see 
Eqs. mi) and (1721) ] induced by the sensor QD onto the 
qubit by capacitive interaction A: they may even be re¬ 
produced by modeling the SQD charge as a classical fluc- 
tuator (see, e.g., Refs. [60,163). I n fact, to obtain Eqs. 
m and (ED, we employ a Markovian approximation 
with respect to the sensor on the qubit (by setting z = 0 , 
see discussion in Appendix IB 3 cl) . The oscillations of 
the qubit decay should also not be mistaken for coher¬ 
ence revivals, i.e., an increase in the purity of the qubit 
state. Figure HTT a! clearly shows that the magnitude 
of the qubit Bloch vector decreases for all times, i.e., 
the information is permanently transferred to the envi¬ 
ronment during the measurement process. However, the 
rate of information loss is nonmonotonic, which simpler 
approaches might not predict (see, e.g., Ref. |Z3). This 
shows most clearly that our perturbation theory goes be¬ 
yond the simple “mean-field” detector picture discussed 
in Sec. IIVBI 

V. COMPARISON TO OTHER APPROACHES 

In this final Section we compare our results for the 
measurement backaction with the results of prior works. 
One of the central results that we obtained - the strong 
suppression of the measurement backaction when tun¬ 
ing the sensor into the Coulomb blockade regime - is 
surprising: One expects to underestimate the measure¬ 
ment backaction by a too simplified treatment which ig¬ 
nores cotunneling noise. We accounted for this cotun¬ 
neling noise contribution (“broadening” contribution to 
the rates) but our result (R)5l) revealed that it is canceled 
out by the coherent backaction, a renormalization effect. 
Such a cancellation should of course be viewed very crit¬ 
ically and we have carefully traced its origins. It also 
raises the question, which physical assumptions and ap¬ 
proximations may have caused it not to be noticed before. 
This is discussed here. 

We first show in Sec. IV Al that semiclassical approaches 
can only reproduce the contributions to the decay times 


induced by the stochastic backaction, but they fail to 
account for corrections from the coherent backaction. 
This also quite easily happens within a standard weak- 
coupling Bloch-Redfield approach aimed at finding a 
qubit-only description as we discuss in Sec. IV Bl Both 
these approaches are correct only under the assumption 
of the weak-coupling, high-temperature limit, which may 
experimentally be violated and does not allow the sen¬ 
sor QD current to be calculated. Finally, we show that 
our findings are in accordance with exact quantum treat¬ 
ments of models of decoherence due to noninteracting 
two-level fluctuators in equilibrium, insofar a compari¬ 
son is possible. 

A. Semiclassical stochastic approaches 

A popular way to study the decoherence a qubit suffers 
from the coupling to its environment is a semiclassical 
qubit-fluctuator model [24|, [2614281 Il62l ] . In this approach, 
the qubit is considered subjected to noise generated by a 
randomly switching two-level system. The basic idea is 
to replace the occupation number n of the SQD in the 
interaction Hamiltonian, Hj = h X ■ t/2, by a classical 
random process h —> £(t)/2 EH- Applied to our case, 
the SQD introduces a fluctuating effective magnetic field 
acting on the qubit [see Eq. ED]: 

H Q = + (105) 

Since our primary interest lies in the interplay of the 
coherent backaction with the cotunncling-induced back- 
action, we can take D = 0 and drop the spin degree of 
freedom of the SQD electrons and therefore the effect 
of the Coulomb interaction in the following. Neither of 
these assumptions are critical for reproducing these ef¬ 
fects (see Appendix iDl). The measurement then induces 
a dephasing of superpositions of the qubit states | L) and 
| R), physically related to the double-QD electron residing 
either in the left or right QD. This can be characterized 
by the decay of the off-diagonal elements of the qubit 
density matrix, the visibility [23j, |H, |6(j: 

(f+(f)> = ((L|e-^ dtlff e f(tl) |A) 

x {R\e +i ^ dtlH Q l{tl) \R)Y (106) 

Here, the bracket (...) denotes the average over many re¬ 
alizations of the random process £(t) with (£(i)) = 2p l st . 
Splitting the qubit Hamiltonian Hq = ( Hq ) + SHq 
according to Eq. ED into a mean-field part, (Hq) = 
\pl t X ■ t, and a fluctuating part, SHq = i<5£(t)A ■ r, 
with 6£(t) = £(t) — (£(t)) and (6£(t)) = 0, Eq. (11061) can 
be recast as 

(f+(t)) = fo * i, 5 5 ( ti )^ . (107) 

Hence, the decay of the coherences depends only on the 
“amplitude” of the fluctuations, while the average of the 
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fluctuations, (£(t)}, is irrelevant. This mean-field part 
just induces a constant shift of the qubit energy level, 
which has been absorbed into the average ( Hq ), see our 
discussion of the mean-field picture in Sec. IIVB II 
For simplicity, let us first take <5£(t) to be a Gaussian 
random process. In this case, Eq. (11071) can be rewritten 
as [m 


■ Jo* dti5£(ti) 


= g- V Jo dtl Jo <5^(i 1 )<5^(t 2 )) 


(108) 


resulting in an exponential decay of (f +(t)) 
The dephasing time, 


1 

T 2 


1 

Ta 


Ssd 0 ) 


= e-^. 


(109) 


can be related to the noise power spectrum Ss^(lo) of the 
Gaussian random process (f) 


/ +00 

dre 

-OO 


r (^( 0 )^(r)), ( 110 ) 


where we used that the correlator (<5£(fi)(5£(f 2 )) only 
depends on the time difference r = £2 — ti because a 
Gaussian random process is time-translational invariant. 
We emphasize that the relation (11081) holds exactly for a 
Gaussian process, whose entire statistics is fixed by two- 
point correlation functions. 

For a general random process, however, higher-order 
time correlators may contribute to Eq. (11081) [M . IM . [fll ] 
and therefore the spectral function is not sufficient to 
characterize a random process completely. A prominent 
example is a Poisson process inducing random telegraph 
noise. Such a process would actually be a better model 
for a capacitively coupled QD stochastically switching 
its occupation due to tunneling processes. Random tele¬ 
graph noise has also been extensively studied to explain 
the ori gin of flicker ( 1 //) noise in superconducting qubits 
, 28|,|62| and Gaussian noise in semiconductor QDs 
I . This noise results from an ensemble of fluctuat¬ 
ing background charges, each of which may be compared 
with our sensor QD, except that they do not carry a cur¬ 
rent on average. The non-Gaussian behavior of the Pois¬ 
son process comes to light only if a few fluctuators dom¬ 
inate the decoherence of the qubit state, which entails a 
nonexponential decay of the coherences [M, [M l60l l 6 ll 
and even coherence revivals have been predicted |60L 16 11 . 

The impact of charge fluctuations of an individual QD 
have also been modeled by classical random telegraph 
process Gum. In general, St; switches between two val¬ 
ues <5£ u (upper), and 5t;i (lower) with different switching 
rates I/77 for 6t; u 8£i, and 1 / t u for S£i 5£ u . This 
gives rise to different probabilities p u to find the value 
6£ u , and pi to find the value 5£i, respectively. In the 
weak-coupling limit, A -C 7 = 1 / t u + I/ 17 , it turns out 
that the deph asing is still exponential and the dephasing 
rate reads as Ill 2 f : 


1 

T 2 


[i-(p«-p 01 


2i (A/2) 2 


( 111 ) 


Identifying p u —> and pi —»■ p® t with the stationary oc¬ 
cupation probabilities of the SQD for our case and not¬ 
ing (A • e ||) 2 = 1 for Q = 0, Eq. (11111) reproduces the 
the dephasing rate (1551) . which we obtained in the high- 
temperature limit. 

Yet, in the classical fluctuator model, we have not ex¬ 
plicitly assumed high temperatures. Thus, if one naively 
extends the above result to include cotunneling correc¬ 
tions, one gets a faulty result. The cotunncling changes 
the occupation probabilities p!? t and p( t and the switch¬ 
ing rate 7 , but there is no way in which additional co¬ 
herent backaction terms could appear [see the transition 
factor r, Eq. As a consequence, one might overes¬ 

timate the qubit dephasing rate at the onset of Coulomb 
blockade. In addition, the classical approach does also 
not account for the dissipative backaction, which is less 
important for the d epha sing times as compared to the 
coherent backaction jll3f . 

The reason why the above classical approach is not able 
to reproduce the coherent backaction is that the electron 
number n is a classical variable with a definite value at 
each instant of time, meaning the SQD is fluctuating be¬ 
tween the states | 0 )( 0 | and | 1 )( 1 |, where |n) denotes the 
SQD state with n electrons on it. In this way, it can 
only produce the stochastic backaction. Quantum coher¬ 
ences | 0 )( 1 | and 11 )( 0 | involved in virtual processes of 
the SQD (quantum fluctuations), which play a role dur¬ 
ing the tunneling, are disregarded here. However, as we 
illustrate in Appendix [Dj these coherences - included in 
our calculation - are crucial for obtaining the coherent 
backaction. (Note that despite this, such coherences can 
not appear in the real quantum state, the relevant den¬ 
sity operator p [see Eq. (ED] due to charge conservation.) 
It turns out that the related quantum fluctuations induce 
an additional “phase kick” while a charge transition in 
the SQD takes place. This phase shift partially com¬ 
pensates the phase shift from the stochastic backaction 
while the SQD is in a specific charge state. This makes 
plausible why quite generally one may expect a coherent 
backaction that mitigates measurement backaction when 
employing a quantum sensor for the indirect detection 
of a qubit. Thus, a purely classical understanding of the 
indirect measurement backaction due to a sensor QD is 
incomplete. 


B. Bloch-Redfield approach 

Another frequently employed method to study deco¬ 
herence is the Bloch-Redfield approach (48[ Ez| • also in 
the context of qubits [M H3 ■ In the original presenta¬ 
tion (67] , this approach is first developed for a semiclassi- 
cal [114j and after this also for a quantum perturbation. 
In both cases, one considers the limit of weak coupling 
between a quantum system and its environment and an 
additional Born-Markov approximation is made. Impor¬ 
tantly, the approach generally predicts an exponential 
decay into a stationary state irrespective of the statis- 
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tics of the environmental fluctuations. The correspond¬ 
ing relaxation matrix of the reduced density matrix is 
furthermore related to the noise power spectrum of the 
perturbation. 


1. Stochastic backaction in spectral function 


For our case, the Bloch-Redfield approach confirms the 
dephasing rate (11091) also within a quantum treatment: 
one simply has to replace in Eq. (II101) the classical aver¬ 
age (<5£(t)<5£(t + t)) by the quantum-ensemble average: 


: — 


f + oo 


dre lu,T 4, (5n(t)6h(t + r)} . (112) 


For an indirect readout model of interest here, one treats 
the SQD as a noninteracting two-level quantum system 
in contact with a thermal reservoir. Applying the results 
of Ref . I~2~> for a single fluctuator, the spectral function 
reads |l!5| 

$sn(w) = [1 - (pit ~P°t) 2 ] 2 1 7 2 ( 113 ) 

7 ^ + oj z 

leading to Eq. (11111) when inserting Ssn(u = 0) into 
Eq. (II101) . Hence, also the Bloch-Redfield approach - 
even accounting for a quantum environment - does not 
reproduce the coherent backaction and may therefore be 
applied only in the high-temperature limit to derive de¬ 
phasing and relaxation times. 


2. Difficulty of capturing dissipative and coherent 
backaction 

The reason why the Bloch-Redfield approach fails to 
account for the dissipative and coherent backaction is not 
specific to this approach. It is rather a problem that ap¬ 
pears in principle for any procedure that tries to treat 
the sensor in an indirect measurement setup as a given 
environment without a nonequilibrium dynamics of its 
own. Such an approach always needs to make some as¬ 
sumptions about the sensor density operator, which, as 
we have seen, does not seem to have a simple structure 
allowing for an educated guess based on general physical 
principles. This means that one should simply calculate 
the joint quantum state of sensor QD and qubit, which is 
what we have done. It is, however, instructive to further 
understand the problems encountered when one tries to 
avoid this by making (too) simple approximations of this 
joint state. 

(i) The Bloch-Redfield approach (as well as many 
other density-operator approaches) applied in the weak- 
coupling limit involves a factorization assumption [68| for 
the state of the qubit and its environment (here the sen¬ 
sor QD plus the reservoir): ptotit) « PQ[t) ® psR{t)- The 
factorization assumption ignores the dissipative backac¬ 
tion valid only in the high-temperature limit. Otherwise, 


the dissipative backaction leads to a nonzero stationary 
qubit Bloch vector that reflects the nonfactorisability of 
the qubit-environment state even in the stationary limit. 
By contrast, to find the decoherence rates in the long¬ 
time limit t — to > 1/r, the factorization assumption 
may still work. 

(ii) The next critical point is then to find a proper de¬ 
scription for the unknown evolution psB.(t) of the envi¬ 
ronment. An assumption frequently made is that the 
qubit environment is in a stationary state, psFt(t) = 
Psr, st, for example, an equilibrium state. However, sta- 
tionarity of the entire qubit environment (SQD plus elec¬ 
trodes) is actually never reached whenever a measure¬ 
ment is performed, for which the electrodes must be held 
at finite bias to produce a nonzero measurement current. 
Here one should be careful to note that while the reduced 
sensor QD density operator may become stationary after 
some time, this is not true for the joint SQD-electrode 
state. In our approach, the reduced sensor-qubit system 
can become stationary after eliminating the electrodes, 
which are stationary. 

(iii) Even if we further simplify the problem and as¬ 
sume zero bias voltage (e.g., to compute the measure¬ 
ment backaction in linear response) and psR,e q is station¬ 
ary, then it is still difficult to compute the equilibrium 
state psR, e q since the sensor QD is a strongly interact¬ 
ing system. Naive assumptions made about psR, e q are 
prone to errors. Consider for instance the approximation 
PSR,e q = PS, st ® Pr, o, where ps,st denotes the stationary 
SQD state and pR$ is the grand-canonical equilibrium 
state of the reservoirs. (One may be inclined to make this 
approximation since weak coupling F often implies such a 
factorization.) If we use this state and average the two- 
point charge correlator in Eq. (I112D . this involves only 
charge-diagonal SQD states |0)(0| and |1)(1| - similar to 
the semiclassical approach of Sec. IV Al What goes wrong 
here is that ps, st ® Pr,o is not the correct equilibrium 
state if we go beyond the lowest-order T approximation 
in the tunneling. For larger T, the hybridization between 
both systems cannot be neglected any more [68| . During 
tunneling processes the total system can be in virtual in¬ 
termediate states involving sensor QD coherences |0)(1| 
and 11)(0| and corresponding charge coherences in the 
electrode. These intermediate virtual states explicitly 
appear in the calculation of the coherent backaction (see 
Appendix [Dj. Thus, even the quantum-ensemble averag¬ 
ing procedure that brings us from Eq. (11121) to Eq. (11131) 
misses the coherent backaction since it relies on a weak- 
coupling expansion between SQD and electrodes. 

This explains why for an indirect weak measurement 
setup, the procedure employed in this paper and also 
in Refs. [23, l30| seems unavoidable. By integrating out 
the electrodes first one incorporates their effect on the 
joint sensor-qubit system. As we have seen, this reveals 
that the qubit experiences a stochastic, dissipative, as 
well as a coherent backaction effect. This problem of 
describing the nonstationary environment is specific to 
indirect measurement setups and not encountered when 
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the qubit is directly coupled to, e.g., a stationary envi¬ 
ronment. This is the case, for example, for a bath of har¬ 
monic oscillators as in the spin-boson model [27, 69j, )70|, 
where the environment may indeed be described by a sta¬ 
tionary equilibrium state. Yet, previous studies for other 
types of environments also show that the approximations 
made to integrate out the entire qubit envronment can 
be too crude, including both a_spin-boson model [97] and 
a driven two-level fluctuator [96]. In particular, Ref. [96j 
explicitly compares results for different levels of approx¬ 
imations showing the breakdown of Markov and secular 
approximation beyond the weak qubit-environment cou¬ 
pling. 


C. Nonperturbative quantum solutions 

Our study indicates that renormalization effects, based 
on quantum coherences between the qubit and its detec¬ 
tor, are vital for the description of the measurement back- 
action even in the weak-coupling limit. We next compare 
our results to prior studies that treat the qubit decoher¬ 
ence arising from single fluctuators coupled to fermionic 
reservoirs fully quantum-mechanically. 

Such studies employ various approaches, such as an 
exact numerical evaluation of the visibility (11061) using 
electrodes of finite size ^flj, a Heisenberg equation-of- 
motion technique HU, or a Keldysh path-integral formal¬ 
ism HU HU . All these approaches are nonperturbative 
both in the measurement interaction A and the sensor 
tunneling T. However, in contrast to our model these 
studies are limited to noninteracting fluctuators (here: 
the SQD) in equilibrium with a single reservoir. Thus, 
they cannot access the situation of a (nonequilibrium) 
signal current through the fluctuator, an essential aspect 
of the indirect detection that we do consider. Moreover, 
they only consider qubit energy splittings along the mea¬ 
surement vector, in our notation Hq = £It z /2, that is, 
they only study pure dephasing in which the qubit Bloch 
vector has no precessional motion. This leads to drastic 
simplifications employed in the derivation of these ap¬ 
proaches which limit their applicability. 

References fiiol [HTj highlight coherence revivals for the 
short-time transients in the strong measurement regime 
1 > T. This nonexponential decay reflects the non- 
Gaussian statistics of the quantum telegraph process. 
Moreover, these studies also find oscillatory corrections 
to the dephasing even in the weak-measurement regime 
A < T in agreement with our work (see Fig. ED. The 
dependence of the dephasing rate on the level position is 
not reported and thus cannot be compared. 

In Ref. [Ill, a path integral method is used to study 
the qubit decoherence due to two-level fluctuators. Their 
dephasing rate also includes terms containing the real 
part of the digamma function if, which determines the 
renormalization function (f [Eq. (1301) ]. It appears both 
in the cotunneling rates [see in Eq. (ED] as well as in 
the coherent backaction [Eq. (l35l) ]. In contrast to Ref. 


HU, our expressions depend on the derivative <j>', which is 
arises from our expansions in A and T. However, a direct 
comparison is not possible because tangible results of Ref. 
HU are actually given only for the sequential tun nelin g 
regime where renormalization effects are neglected [l 16f | . 

Reference HU also uses a path-integral approach and 
considers also the case when the QD level e is tuned 
strongly away from resonance [e — n\/T 1. Expanding 
their result (Eq. (15) of Ref. HU) f° r th e long-time limit 
of the depha sing rate in the weak-measurement limit, we 
obtain |ll7[ : 


1 


T A 2 T 2 
47 t (e — p) 4 


+ 0((A/T) 4 ). 


(114) 


Thus, the dephasing rate drops algebraically in the fourth 
power with the level position. This is consistent with our 
results in the sense that also here the decoherence does 
not scale as expected from the cotunneling noise. In our 
case, we were only able to show that the algebraic scaling 
in l/(e— fj.) is suppressed at least to the second power (see 
Sec. lIIID II) . However, our analysis revealed the physical 
origin of this behavior by identifying what cancels the 
expected cotunneling noise contribution, namely the co¬ 
herent backaction. Note that the result (11141) does not 
rely on the high-temperature assumption T/T < A/T as 
in our case, which is why their result can be expanded 
in orders A/T while accounting also for higher-order tun¬ 
neling corrections [e.g., 0(T 3 /T 2 ) and higher, see Ap¬ 
pendix IB 4| . This implies that to compare concretely 
with their result (I114D , we would have to include higher- 
order tunneling terms, which, however, we do not expect 
to restore power laws with lower exponents. This is chal¬ 
lenging for our model since we account for Coulomb in¬ 
teractions and nonequilibrium conditions. This is beyond 
the scope of the present paper. 

In summary, various aspects of our assumptions and 
findings seem to be in accordance with previous approx¬ 
imate as well as exact quantum-mechanical treatments 
and shed new light on them. Our study extends these ap¬ 
proaches by simultaneously dealing with a nonstationary, 
nonequilibrium, Coulomb-blockaded sensor QD (fluctua¬ 
tor), which is fully quantum-correlated with the qubit 
(non-factorizing density operator p, virtual off-diagonal 
charge coherences during tunneling). The latter leads to 
the coherent backaction as an integral part of the total 
backaction together with stochastic and dissipative back- 
action, leading to the cancellation of cotunneling noise. 


VI. SUMMARY AND OUTLOOK 

We studied an indirect detection setup, in which a 
charge qubit is capacitively probed (A) by a sensor quan¬ 
tum dot (SQD). The SQD is in turn tunnel coupled (T) 
to electrodes in which the time-dependent conductance 
is measured. Electrons in the sensor QD occupy a single 
quantized orbital in which they strongly interact. 
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Kinetic equations. We considered the weak-tunneling, 
weak-measurement limit T/T -C A/r -C 1, in which 
quantum fluctuation effects are important on the time 
scale of the qubit-sensor interaction. We derived a ki¬ 
netic equation [Eq. (1261) ] by integrating out the current- 
carrying electrodes to obtain an effective theory for the 
composite qubit-SQD system. This revealed three types 
of backaction on the qubit: (i) a stochastic backaction 
due to random fluctuations of the qubit detuning, (ii) a 
dissipative backaction (coefficient c), the flip side of the 
modulation of the sensor tunnel current by the qubit, 
and (iii) a coherent backaction due to the level renormal¬ 
ization of the composite qubit-SQD system (coefficient 
«)• 

We showed the importance of the effects of single¬ 
electron tunneling (SET), as well as its cotunneling 
broadening and level-shift corrections. We also included 
the leading non-Markovian correction from the electrodes 
induced by the tunnel coupling (T) to the SQD (linear 
kernel frequency dependence). Moreover, our approach 
captures all non-Markovian effects introduced by the sen¬ 
sor QD on the qubit subsystem, which are mediated by 
the capacitive interaction. 

Suppression of cotunneling-induced backaction. By 
rewriting the kinetic equation in the basis of quasista¬ 
tionary and decay modes (defined by the A = 0 limit), 
we found that the interplay of these types of backaction 
leads to a nontrivial cancellation; whereas the dissipative 
backaction (c) independently couples these modes, the 
stochastic backaction and coherent backaction (re) par¬ 
tially cancel. In particular, the change in the stochastic 
backaction due to cotunneling broadening is canceled by 
the coherent backaction [Eq. (1521) ]. The expected alge¬ 
braic decay oc l/(e — p r ) of the backaction (determining 
the decoherence rates) is thus suppressed, implying that 
the actual power law must at least have a higher expo¬ 
nent. Experimentally, this is important since it indicates 
that a SQD can be switched off by applying a gate volt¬ 
age better than expected. By identifying the underlying 
physics, we suspect this to be a crucial difference with 
the backaction of sensors with dense level spectra, such as 
single-electron transistors. Thus, the less-than-expected 
backaction due to the coherent backaction is beneficial 
for switching the sensor on / off, provided one takes care 
to prepare sensor state to avoid initial slip errors (see 
next paragraph). 

Initial slip. We derived effective equations for the re¬ 
duced qubit density operator [Eq. (1711) ]. which are exact 
relative to the kinetic equations (1261) . In particular, we 
keep all non-Markovian effects induced by the sensor and 
account for a slip of the initial condition [Eq. (1721) ] . This 
slip depends on both the initial sensor QD state and 
the initial quantum correlations with the sensor QD. It 
is important for the long-time qubit evolution and thus 
the sensor QD needs also to be considered as part of the 
dynamical quantum circuit. The dynamical state of the 
sensor QD is relevant for qubit error propagation, e.g., 
the initial slip may introduce errors even for perfectly 


prepared qubit initial states. This provides a concrete ex¬ 
ample for errors usually phenomenologically introduced 
in quantum-error correction. Such a sensor QD cannot be 
considered (without further evidence) as a “black box” 
in a quantum circuit which is merely characterized by 
static parameters (e.g., relaxation and dephasing times). 
This is different from the treatment of the macroscopic 
electrodes coupled to the sensor QD. The electrodes, in 
which the current measurement is performed, can instead 
be assumed to be stationary, which eliminates initial-slip 
effects on qubit-SQD system (provided initial quantum 
correlations are neglected). 

High-temperature limit. Specializing to the high- 
temperature limit T/T —^ 0 and to times larger than 
the SQD relaxation time 1/T, we obtained the qubit evo¬ 
lution [Eq. (I8()[) ]. which neglects non-Markovian effects 
induced by the sensor QD. We connected dephasing and 
relaxation times oc A 2 /T with the component of the mea¬ 
surement vector A along the mean-held qubit axis f2 and 
perpendicular to it, respectively. We demonstrated the 
importance of the initial slip oc A/T [Eq. (1511) ] even in 
this simplest limit: the set of initial states of the qubit- 
SQD system without slip is only a subset of zero mea¬ 
sure. Generally, the magnitude of the slip increases the 
“less factorizable” the qubit-SQD state is and the “more 
nonstationary” the sensor QD is before the detection is 
started. Due to the latter, the initial slip should be reck¬ 
oned with in particular when the sensor QD is initial¬ 
ized in a fixed charge state. This happens, for example, 
when switching the sensor off by tuning the gate volt¬ 
age far away from resonance or by reducing the tunnel 
coupling T to the electrodes. By contrast, switching off 
the capacitive interaction A leads to an initially station¬ 
ary and factorizable state, which is favorable for avoiding 
an initial slip. This difference in the backaction between 
various parameters for switching off a sensor QD is an 
important experimental implication of our study. The 
need to control not only the qubit but also the readout 
device carefully may not only be a nuisance for engineer¬ 
ing quantum circuits but could also provide additional 
means of controlling qubits. For example, one could con¬ 
sider to switch the SQD during the readout to another 
readout point to compensate for manipulation errors de¬ 
tected in a weak measurement process. 

Mode distortion. The analysis of the isospin dynamics 
showed that additionally the qubit eigenmode vectors are 
distorted due to the small but finite value of the decoher¬ 
ence rates ~ A 2 /T. This corrects the simple mean-field 
picture, in which the qubit axis Tt = f2 + Pg t A is only 
influenced by the average occupation of the SQD. The 
distortion reflects the breakdown of the secular approxi¬ 
mation - often made in derivations - because the capac¬ 
itive coupling A can be of the same order as the internal 
qubit splitting Tl. Importantly, the distortion must be 
included to satisfy an isospin sum rule that follows from 
the conservation of the isospin by the tunneling of sen¬ 
sor electrons [H, Hgj. The distortion is thus enforced by 
a general principle. The experimentally relevant conse- 
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quence of this distortion for the qubit evolution is two¬ 
fold: first, the usual circular Bloch-vector precession be¬ 
comes slightly elliptical , with a shape that is not altered 
in time as the size shrinks with dephasing rate I/T 2 . Sec¬ 
ond, this “precession plane” is not orthogonal any more 
to the “relaxation axis,” along which the Bloch vector 
decays with the relaxation rate 1/Ti. The relaxation 
axis is moreover slightly tilted relative to the mean-field 
axis ft. Finally, the projections of the isospin on the 
mean-field axis, the relaxation axis, and the precession 
plane all show a superposition of relaxation and preces- 
sional dephasing motion. The measurement backaction 
thus mixes the effects of relaxation and dephasing, even 
in this simple Markovian limit. This effect may generally 
appear in indirect-coupling setups, which are typical for 
detection setups. 

The tilting of the qubit modes is related in a broader 
context to the concept of “gauge qubits” [98j: two phys¬ 
ically distinct qubits can be each subject to strong de¬ 
coherence but the joint Hilbert space formed by both 
qubits may contain a two-level subsystem with low deco¬ 
herence. Locating such subsystems is interesting for de¬ 
veloping strategies for quantum-error correction. In our 
case, the qubit mode tilting reflects an “admixing” of 
the sensor QD degrees of freedom to the low-decoherence 
subspace, which mostly overlap with the qubit degrees 
of freedom. This mixing can be strongly enhanced when 
the coupling of the quasistationary and decay modes be¬ 
comes stronger, which effectively happens, e.g., due to 
enhanced quantum-fluctuation effects at lower tempera¬ 
tures. Investigating this mixing further would therefore 
be interesting also in the context of quantum-error cor¬ 
rection strategies. 

Comparison with other approaches. We compared the 
above results at various points with existing approaches 
and explained why potential differences are expected 
within their validity: 

(a) Semiclassical stochastic approaches cannot capture 
the coherent backaction because the starting assumption 
of classical charge fluctuations on the sensor QD already 
excludes relevant quantum coherences of the qubit-sensor 
density operator. Including these into the description 
leads to additional “phase kicks” that counteract the 
stochastic “phase kicks.” Both together, when averaged, 
lead to a mitigated decoherence. It thus seems that one 
should quite generally reckon with coherent backaction, 
which can mitigate the measurement backaction when 
employing a quantum -dot sensor for the indirect detec¬ 
tion of a qubit. 

(b) Density operator approaches that try to integrate 
out the sensor together with the attached electrodes run 
into problems as well because one needs to “guess” the 
time-dependent, current-carrying sensor state as well as 
the quantum correlations with the qubit. These can lead 
to non-Markovian behavior and affect the initial slip of 
the qubit. All of this is systematically calculated in our 
approach. The more advanced approach of Ref. H), 
which first calculates the nonequilibrium SQD state in 


the absence of the qubit (an approximation carefully 
pointed out in Ref. (3), misses both the dissipative 
and coherent backaction. Extending such an approach, 
e.g., to include cotunneling broadening would thus lead 
to inconsistencies since only one part of two canceling ef¬ 
fects is taken into account. 

(c) Nonperturbative quantum solutions of related mod¬ 
els agree with the cancellation in the backaction we found 
here as does a separate calculation for a noninteracting 
limit of our model (U = 0). 

(d) It is also interesting to compare with prior studies 
not aiming at sensor backaction. In particular, Ref. [78| 
highlighted the importance of a competition of next- 
to-leading-order effects and non-Markovian corrections, 
which are closely related to the non-Markovian correc¬ 
tions of the first type affecting the dissipative switching 
rates © discussed in Sec. IIII B dl However, what is 
imporant here are non-Markovian corrections of the sec¬ 
ond type affecting the coherent backaction. We stress 
that although these non-Markovian corrections give an 
important part of the coherent backaction, they are not 
identical to it: the coherent backaction already arises in 
the stationary limit that we studied in Ref. [32) where non- 
Markovian effects can be neglected. In this way, we see 
that the general line of thought emphasized in Ref. [78| ex¬ 
tends to the coherent backaction, which is also a second- 
order effect, namely first order both in the tunneling (T) 
and the measurement interaction (A). 

Outlook. All this shows that the coherent backaction 
and other quantum fluctuation effects (cotunncling, level 
renormalization) are intrinsic effects of a quantum sen¬ 
sor: they can neither be “added” or “controlled” indepen¬ 
dently in an experiment, nor should they be neglected in 
a calculation. In view of the above, further studies that 
can address the experimentally relevant lower temper¬ 
ature dynamics of weak measurements using quantum- 
dot sensors are necessary. We expect the qualitative ef¬ 
fects that we identified to be present and quantitatively 
stronger under experimental conditions where, e.g., the 
capacitive interaction and tunnel coupling may not be 
weak anymore (e.g., because of the trade-off between a 
significant current signal and low backaction). Besides 
the commonly discussed relaxation and dephasing back- 
action, also the qubit initial slip, nonorthogonal mode 
distortion, and sensor-induced memory effects will be en¬ 
hanced, neither of which has received much attention so 
far. 
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Appendix A: Importance of non-Markovian 
corrections induced by the electrodes 

In this Appendix, we explain how to account for non- 
Markovian corrections arising from the memory induced 
by the electrodes on the qubit-sensor subsystem and il¬ 
lustrate their importance for a correct description of the 
detector backaction. One should clearly distinguish from 
this further non-Markovian behavior induced by the SQD 
- with the effect of the electrodes incorporated - on the 
qubit subsystem. This is discussed separately in Ap¬ 
pendix IB 3 cl In Appendix IA 11 we show how the non- 
Markovian corrections due to the electrodes can be in¬ 
corporated within the real-time diagrammatic formalism 
based on a perturbative weak-tunneling (T) expansion. 
Within the leading non-Markovian correction included in 
this paper, the kinetic equation remains a time-local and 
first-order differential equation for the density operator. 
Based on this, we then perform the weak-measurement 
expansion. In Appendix IA 21 we show that neglecting 
the non-Markovian correction leads to a violation of the 
positivity of the SQD-qubit density operator and to an 
overestimation of the measurement backaction. These 
unphysical features are removed when the leading non- 
Markovian correction is included. 


1. Incorporating non-Markovian effects in 
real-time diagrammatics 

In Ref. [32}, we started from a kinetic equation for the 
reduced density operator p(t) of the joint system of SQD 
plus qubit, obtained in the standard way by integrating 
out the electrodes’ degrees of freedom: 


with W = W(z = iO). 

Weak-tunneling expansion. Yet, in the present paper, 
we study the nonstationary time evolution of the qubit¬ 
sensor density operator p(t) and we show now that this 
implies that memory effects induced by the electrodes 
through the kernel W have to be included. We will refer 
to all the effects of the frequency dependence of W(z), 
Eq. (TOIl . as non-Markovian corrections induced by the 
electrodes [even when effectively a time-local equation re¬ 
sults, see Eq. (IA.71) ] . An approach to include such cor¬ 
rections within the real-time diagrammatic formalism has 
been given in Ref. jllj, and applied to study decay rates 
in Refs. [71, [z2]- The basic idea is to perform a Tay¬ 
lor expansion of p(t') in the integral in Eq. (IA.1I) around 
time t , 


P(t') 


OO 


E 

n=0 


1 d n p(t) , 
n\ dt n 


t) 


(A.4) 


and to subsequently perform the integration over t'. 
This results in a well-defined expression if the kernel 
decays faster than any polynomial in t — t', which is 
usually fulfilled because the kernel decays exponentially 
W(t — t') ~ Tc on the time scale of the inverse 

temperature tq ~ 1/T |l!8j ] . 

If we consider only times t r c , one may compute 
fg dt'Wit - t')[t' - t) n = f* drW (t)(— r) ra by replacing 
t —t 00 on the right-hand side since all contributions 
from t > t are negligibly small. Taking advantage of the 
Laplace transform (IA.2I) . we obtain from Eq. (IA.1I) 


dp(t) 

dt 


-HQsp(t ) 


V —d n w , 

' n\ dt 


d n p(f ) 


(A.5) 


P(t) 


-HjQspit) + f dt'Wit —t')p{t'). (A.l) 

Jo 


with the nth derivative of the Laplace-transformed kernel 
(IA.2I) at zero frequency with respect to — iz: 


Here, Lqs • = [Hq+Hs+Hi , •] is the internal Liouvillian 
of the reduced system with denoting the operator 
the Liouvillian acts on. The effects of the electrodes are 
incorporated in the kernel W, which we evaluated using 
the real-time diagrammatic approach [5(1 [Hi. In Ref. 
[32! , we were only interested in the stationary solution of 
Eq. (IA.1I) . p{t) = p st , which obeys p st = 0. This solution 
depends only on the time-integrated kernel f* dt'Wit — 
t') = fgdrW(T). In the long-time limit t —> 00, when 
the stationary state is approached, the time-integrated 
kernel is given by the zero-frequency limit of the Laplace 
transform 

pOO 

W{z) = / dre izT W(r). (A.2) 

Jo 


d n W 


d n W{z) 


d{—iz) n 


z—iO 


(A.6) 


To include the leading-order non-Markovian corrections, 
we truncate the sum on the right-hand side of Eq. (IA.5I) 
and keep only the terms with n ^ 1. This yields the 
following time-local differential equation of first order in 
time: 


dT = T^dw { ~ lLQS + w)p{t) = ~ iL P®^- 7 ) 


whose solution can be written as 


p[t) = e lLt pi 0) (A.8) 


Hence, the exact p s t is the same as the stationary so¬ 
lution of the approximate, time-local Markovian kinetic 
equation, 

= i~iL QS + W)pit), (A.3) 


for initial state p(0). When dropping the non-Markovian 
correction 1/(1 — dW) in Eq. (IA.7I) . one recovers the 
Markovian generator of Eq. (IA.3I) . 

We next show that this solution with L defined through 
Eq. (IA.7I) accounts for all non-Markovian effects up to 
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0{T 2 /T) provided one is seeking for exponential solu¬ 
tions of Eq. (1A.5I) . We do not discuss algebraic or loga¬ 
rithmic time dependencies here, which may also appear 
[73j . Thus, let us substitute the exponential ansatz 

Pit) = e At p{ 0), (A.9) 

into Eq. (1A.5I) . We obtain the following equation for A: 

A = —Hqs + W + (d 1 W)A + ... . (A.10) 

When all derivatives d n W are dropped from the above 
equation, we recover the Markovian generator A = 
— Hjqs + W. Non-Markovian corrections therefore en¬ 
ter through terms ~ d n W that capture the frequency 
dependence of the kernel [78l - [80j . 

Equation (IA.9D shows that the derivatives of the den¬ 
sity operator scale as d n p/dt n ~ A n ~ T” since Lqs ~ 
A -C T and W ~ T here. As a consequence, the 
expansion (IA.4D of the density operator in corrections 
from higher-order time derivatives is not independent 
of the perturbative expansion of the kernel W in pow¬ 
ers of r. When expanding A in powers of T, we find 
that the nth order derivative of the kernel scal es a s 
d n WA n ~ (T/T n )-r n plus higher-order corrections [l!9j . 
Thus, when expanding the kernel W = W r + ... only 
up to first order in T, Markovian corrections must be 
ignored, i.e., A r = — iL + W r . Yet, when expand¬ 
ing up to order T 2 /T as considered here, one is al¬ 
lowed to omit terms (d n W)A n ~ (T /T n ) ■ T" only for 
n ^ 2 from Eq. (IA.10D . while one has to keep the first 
derivative d 1 W. Solving for A, we obtain A = —iL = 
(1 - d^W^i-iL + W) + 0(T 3 /T 2 ) in agreement with 
Eq. (IA.7I) . Further expanding the inverse and the kernel, 
we find for the effective generator: 

-iL = (1 + d 1 W r )(-iL QS + W r ) + W r2/T 

+(T 3 /T 2 ). (A.ll) 

Weak-measurement expansion. The next step is the ex¬ 
pansion of —iL in A/T, where A ~ A ,fl denotes the 
small energy scale of the detection and intrinsic qubit 
frequency. While we have to keep the first-order terms 
in A /T for 

W r = w r ’° + W r ’ A + 0(TA 2 /T 2 ), (A.12) 

we can neglect the A dependence of 

d l W r = dW r ’° + 0(TA/T), (A.13) 

W r2/T = W r2/T ’° + 0(T 2 A/T 2 ), (A.14) 

and obtain for the effective, non-Markovian generator 

-iL = (1 + d 1 W r ’ 0 )(—iL QS + W T, °) 

+W r,A + W r2/T ’° 

+0(T 3 /T 2 ,T 2 A/T 2 ,rA 2 /T 2 ). (A.15) 

We now note that in Ref. [74j, we already computed the 
kernels W r '°, W r,A , and W r ,0 and therefore we have 


to merely evaluate dW T '° to obtain non-Markovian cor¬ 
rections up to the order considered here. The frequency 
derivative of <9W r ’° [see Eq. (IA.6I) ] can be easily con¬ 
verted into energy derivatives (d/de) [l20| . which results 
in the kinetic equation (1261) in the main text. It should be 
noted that dW T,t) is not simply the derivative of the SET 
contribution: the imaginary factor in Eq. (IA.6I) changes 
the role of imaginary and real parts, which are related 
to S functions and principal-value parts in the frequency 
integrals, respectively. While in the SET contributions 
only the S functions remain and the principal-value parts 
cancel out, this is opposite for the non-Markovian cor¬ 
rection dW T, °. These principal-value parts evaluate to 
the renormalization function <p(x) [see Eq. (l30l) ]. which is 
of central importance in our work and explains how the 
non-Markovian corrections can affect the coherent back- 
action as noted in Sec. IIII B 41 To evaluate dW T ’°, one 
thus must first compute W r, °, then apply d and only af¬ 
ter that take the relevant matrix elements (restricted by 
charge conservation). This completes the derivation of 
our non-Markovian Liouvillian, accounting consistently 
for all terms up to 0(A, T,T 2 /T, TA/T). 

To find the solution (IA.9I) of the kinetic equation 
p(t) = Ap(t ), we directly solve for p(t) without ex¬ 
panding p(t) i n T / T. Not solving the kinetic equation 
order-by-order 121| avoids well-known problems with ill- 
defined coherences UM and nonequilibrium occupations 
close to the Coulomb blockade regime [5lj| . As a con¬ 
sequence, our solution can comprise of terms of higher 
order in A/T and T/T. In particular, the stationary 
solution obtained from Ap st = 0 differs formally from 
that obtained by solving the corresponding Markovian 
equation Agw-oPst = 0. However, the deviations are of 
0(T 3 /T 2 ,T 2 A/T 2 ,rA 2 /T 2 ) and therefore consistently 
negligible in the perturbative limit considered here. 


2. Retaining positivity and effect on coupling of 
quasistationary and decay modes 


In this appendix, we elucidate the importance of non- 
Markovian corrections induced by the electrodes to de¬ 
scribe the indirect detection setup. We explicitly illus¬ 
trate that without these corrections the kinetic equations 
ESI) for the SQD-qubit system possess exponentially in¬ 
creasing solutions in time and the measurement backac- 
tion is qualitatively different. 

In Ref. flU, we reported that the Markovian kinetic 
equation when used to calculate a time-dependent solu¬ 
tion for the density operator p(t) for the joint SQD-qubit 
system violate the positivity condition, even though the 
stationary state showed no such problems. (It should 
be noted that for the stationary state to be positive, in 
addition to the dissipative and coherent backaction also 
0(T 2 /T) effects were required, already hinting at the can¬ 
cellation effect that we discuss in the present paper.) The 
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FIG. 12: (a) and (b): Two largest real parts bo,bi of the 
eigenvalues of the SQD-qubit generator — iL, as a function 
of gate voltage V g . We use the Markovian approximation 
(1A. 161) of the generator —iLgw =o in (a) and the full non- 
Markovian generator (IA.15I) in (b). (c) and (d): Component 
of r [Eq. (IA. 171) 1 along A. We show the projection tm = tm • 
A = p° t plt + (1 +Pst) [ see Eq. (IA. 181) 1 within the Markovian 

approximatio n in (c ) and r N M = r N M ■ A = p° t plt ~ «(Pst/2 - 
Pat) [see Eq. (IA. 191) 1 including non-Markovian corrections in 
(d). In (c) and (d), the blue curves include the coherent 
backaction (k ^ 0), while the green curves exclude them by 
hand ( k, = 0). The parameters in all plots are 14 = 0, 14 = 
F d = f = 10 _2 T, A = n = O.lf = l(r 3 T, and W = 1000T. 


positivity problems for time-dependent solutions arising 
in the Markovian approximation with respect to the elec¬ 
trodes can be readily inferred by diagonalizing the gen¬ 
erator (IA.15D in that approximation: 

-iL aw = o* = y^ibi + ia,i)Ai tr [Ai»], (A.16) 

S,Q 

i 

Positivity violation occurs when at least one eigenvalue 
exists with a positive real part, i.e., > 0. This leads 

to exponentially growing contributions p(t) ~ p s t + e bit Pi 
[32|. Figure fl2l al shows the two largest real parts of the 
eigenvalues of the Markovian generator (IA. 161) for typical 
parameters considered in this paper. It illustrates that a 
positivity violation may appear when the SQD is tuned 
deep into the Coulomb blockade regime for V g /T 1 
[see Fig. 7(b) in Ref. [32|]. By contrast, if we in¬ 
clude the leading non-Markovian correction and use our 


full equation (IA.15I) instead, all eigenvalues have a non¬ 
positive real part for all gate voltages as shown in Fig. 
fl2T b'). Thus, no positivity violation can occur here. We 
investigated this thoroughly by numerically exploring a 
large parameter regime. We note, as already pointed out 
in Ref. (32), that the inclusion of next-to-leading order 
r 2 /T corrections in the tunneling is also crucial to avoid 
such positivity problems. 

To assess the importance of non-Markovian corrections 
for the qubit backaction, we next consider the effective 
Liouvillian (1691) when starting from the Markovian ap¬ 
proximation of the generator (IA.16I) . As we discussed 
in Sec. nr C 31 the importance of the coherent backac¬ 
tion can be assessed from the transition matrix from the 
quasistationary into the decay modes. In both cases, the 
transition matrix takes the form 

A dq _ ( 0 (1+p° t )cA-e Q / \ 

V ( 1+ P°t)4 ’ ( cA ) -e£-(rxe a ,) J 

(A.17) 

where in the Markovian approximation the vector r reads 

r M = n{l+p° st )n + [p“ t Pst + k |(1 +P° t )] A (A.18) 

and when including the leading non-Markovian correc¬ 
tion [see Eq. the vector r reads 

tnm = [p%plt ~ k (|Pst - P°t)] A, (A. 19) 

which is the central result (1621) discussed in detail in the 
main text. The non-Markovian corrections lead to two 
important differences between Eqs. IA. 181) and (IA. 191) . 

First, the coherent backaction leads to a different gate- 
voltage dependence of the transition matrix when non- 
Markovian corrections are neglected. To illustrate this, 
we plot in Figs. &) and [Hd) the component of 
tm and tnm along the measurement vector A. Due to 
the cancellation effect explained in the main text, the 
coherent backaction suppresses the transition factor in 
the non-Markovian case, leading to an exponential gate- 
voltage dependence. This is different for the Markovian 
case; here, the term arising from the coherent backac¬ 
tion, k( 1 +Ps t ), changes sign at resonance and enhances 
the transition rate into the decay modes for positive gate 
voltages [see Fig. Had)]. Moreover, the transition factor 
scales there algebraically with V g in stark contrast to the 
exponential dependence in the non-Markovian case. 

The second difference between Eqs. (IA.18I) and (IA. 191) 
is that the vector tm in the Markovian case also has a 
component along $7. This component emerges because 
in the Markovian approximation the coherent backaction 
appears as torque terms in the kinetic equations that are 
proportional to k(£1 + A/2) [see Eqs. (25) and (26) in 
Ref. [32]]. By contrast, in the non-Markovian case the 
torque terms involve only the vector A [see Eq. (12511 ]. 

Both these differences clearly show that non- 
Markovian corrections induced by the electrodes have 
crucial consequences for the total backaction due to the 
capacitive interaction A and must be accounted for con¬ 
sistently. 
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Appendix B: Liouvillian perturbation theory 


a. Liouville-space projection technique 


In this appendix, we derive the effective Liouvillian 
(l69l) that we investigate in the main text to study the 
measurement backaction on the qubit evolution. We first 
keep our formulation as general as possible to bring out 
the generic features and to indicate that this procedure 
can be applied also to more complicated indirect detec¬ 
tion models. An important prerequisite is that one can 
identify a subspace of interest in Liouvillc space whose 
dynamics takes place on a time scale that is well sepa¬ 
rated from the dynamics in the complementary subspace. 
In our situation, this is related to the slow dynamics 
~ 1/A of the subspace of the quasistationary modes as 
compared to the fast dynamics ~ 1/7 in the subspace of 
the decay modes. Using the well-established projection 
technique 0, 57] along the lines of Chapter 17 in Ref. 
[55}, we outline how to obtain an effective Liouvillian 
that mediates the dynamics in the subspace of interest 
only. In contrast to most cases found in the literature, 
the unperturbed problem onto which we project already 
exhibits non-Hamiltonian dynamics, i.e., L* ^ [H, •], see 
discussion after Eq. (IB.3711 . Therefore, we review these 
steps here to highlight that the initial slip is a general 
phenomenon that appears when projecting out a comple¬ 
mentary subspace. This in fact prevents a full elimina¬ 
tion of these degrees of freedom unless special conditions 
apply. In the high-temperature limit, the projection can 
be analytically performed and we can obtain the effective 
Liouvillian and initial slip perturbatively in the coupling 
strength between the relevant subspace and its comple¬ 
ment. 

We give here two complementary approaches. The first 
one is given in frequency space in Appendix IB 11 which 
allows for a very compact procedure. However, to un¬ 
derstand the physical meaning of the involved approx¬ 
imations, we also show how to obtain our results in a 
time-space formulation in Appendix IB 21 In Appendix 
IB 31 we then apply our general Liouville-space projection 
technique to the detection setup considered in the main 
text. We finally derive the effective transient evolution of 
the qubit, Eqs. G9 m, by removing the SQD degrees 
of freedom as far as possible. 


1. Derivation in frequency space 


We start with the derivation in frequency space by ap¬ 
plying a general projection technique. This results in the 
most general expression for the time evolution, for which 
we give a perturbative expansion subsequently. 


Consider a system with density operator p whose time 
evolution is generated by a Liouvillian L, i.e., 

p{t) = - iLp{t ), (B.l) 

with initial condition p(0) given at time t = 0. We refer 
to L as the Liouvillian although it might not have a sim¬ 
ple commutator structure and can have dissipative, non- 
Hermitian parts. To solve the above equation, one can 
transform it to Laplace space. The Laplace transform of 
a function /(f) is defined as f(z) = / 0 °° dte lzt f (t). This 
yields using Eq. (IB. ID : 

— izp(z) — p(0) = —iLp(z). (B.2) 

Next, we are only interested in the evolution of p in a 
subspace a of the entire Liouville space, defined by a pro¬ 
jection superoperator P“, which satisfies ( P a ) 2 = P a . In 
contrast to Sec. IIII Cl our formulation here is completely 
basis independent to emphasize the generality. 

Projecting the kinetic equation (IB. 2D onto the subspace 
a and its complement b with projector P b = 1 — P a , we 
find 

~izp a (z)-p a ( 0) = ~iL aa p a (z)-iL ab p b (z){ B.3) 
-izp\z) - p\0) = -iL ba p a (z)~iL bb p b (z).(BA) 

Solving the second equation for p b (z), 


p\z ) = j^ I5E (p b (0)-iL ba p a (z)), (B.5) 


and inserting this into the first equation, we can formally 
write the exact solution as 


p a (z) = 


Z - L ts( Z ) 


Pes(z)- 


(B.6) 


This incorporates a frequency-dependent effective Liou¬ 
villian, 


L a £{z) = L aa + L ab —^L b \ (B.l) 
and the frequency-dependent initial condition 

P a eS (z) = p a ff)) + L°- b —±- rb p b ( 0). (B.8) 

To transform Eq. (IB.6D back to time space we apply the 
inverse Laplace transform: 

r+°° rl? 

p a it) = J ^ -e-^p a (z). (B.9) 

Identifying the subspace a ( b ) with the subspace of the 
quasistationary (decay) modes labelled by q ( d ), Eqs. 
(IB.6D (IB.SD yield Eqs. (I55D (T70l) of the main part. 

To compute the solution in time space, the integral (IB.9D 
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has to be solved by applying the residue theorem. Since 
L^q(z) can be represented by a finite matrix in our case, 
we only have isolated poles satisfying z p = L^g(z p ). If 
the coupling between the subspaces a and 6 is absent, 
these poles are given by the eigenvalues of L aa . 

Importantly, Eq. (IB. 81) shows that the component of 
the initial state in the complementary 6-subspace has 
been transformed into a correction to the initial state 
in the targeted a subspace. This is referred to as a 
slippage of the initial condition Q- In general, if the 
coupling of the 6-subspace to the a-subspace is nonzero 
( L ab ^ 0), this initial slip is present unless the initial 
state has no projection on the complementary space, 
p b ( 0) = P b p( 0) = 0. In our problem, P b projects onto 
density operators p 1 which are either not factorizable or 
for which the reduced state of the SQD is not station¬ 
ary (see Appendix [Cj. This has to be compared with the 
generalized master equation (EH) for the joint qubit-SQD 
state: As noted there, when deriving the general kinetic 
equation one in fact projects onto a stationary (equilib¬ 
rium) state of the reservoirs. One then often assumes (as 
we do) an initially stationary reservoir state that factor¬ 
izes with the SQD-qubit state. This eliminates the initial 
slip. 

b. Perturbative solution 

To next find a perturbative solution - required to ob¬ 
tain the high-temperature equations (15U1) and (ISlT) of the 
main text - we assume that the eigenvalues of L aa are 
well separated from those of L bb in the complex plane 
as compared to the coupling mediated by L ab and L ba . 
Thus, if 

l :=max(|L ab |,|L bQ |) < g := \\L bb \ ^ \L aa \\, (B.10) 

we can neglect the coupling-induced shift of the poles z p 
in the denominator of Eq. to lowest order in l/g. 

This is the basis of our perturbative expansion of L'/g in 
orders of l/g. 

For the measurement setup studied in this paper, the 
situation is even a bit simpler because the eigenvalues in 
the quasistationary subspace satisfy 

w := \L qq \ « 9, (B.ll) 

which allows us to insert z = 0 into Eqs. EZ) and (EH). 
This means that we carry out a Markovian approxima¬ 
tion (i) for the effective Liouvillian 

= L“%(z = 0) = L aa - L ab j^L ba (B.12) 
and (ii) for the effective initial condition: 

Peff(0) = P a { 0) - L ab ^P b ( 0). (B.13) 

One can next readily transform back to time space, which 
yields the effective time-evolution equation, 

p a (t) = -iL^p a (t), (B.14) 


with the “slipped initial condition” (IB. 131) . Equations 
(IB. 121) (IB. 141) reproduce Eqs. (1791) (1511) of the main 
text. 

We emphasize that this Markov approximation is pro¬ 
foundly different from the Markov approximation dis¬ 
cussed in Appendix [XJ The latter accounts for memory 
effects of the electrodes (kernel frequency dependence), 
which are already contained in L and therefore included 
into the projection approach from the start. In the 
present case we neglected the frequency dependence in 
Eqs. (IB.81) and (IB. 71) that is generated by considering a 
subsystem of the system described by L. We address this 
point further in Appendix IB 31 where we discuss this for 
our concrete detection problem. 

The Markovian approximation (IB. 141) neglects correc¬ 
tions to the effective Liouvillian of 0([max(Z, w)] 3 /7 2 ) 
and therefore the solution of Eq. (IB. 141) . p a {t) = 
e -* i e« t p“ ff (o) ; is limited to times t <C 7 2 /[max(Z,w)] 3 . 
Moreover, there is a restriction for Eq. (IB. 141) to be valid 
for small times: one has to require t > I/7 basically 
to ensure that the projection of the density operator on 
the decaying subspace, p b (i), becomes negligibly small 
and effectively contributes only through the slipped ini¬ 
tial condition Eq. (IB. 131) . This becomes clearer from our 
complementary time-space derivation below. 


2. Derivation in time space 

To gain further physical insight into the approxima¬ 
tions made in the above derivation, we re-derive the 
above results now in time space. We start here from 
decomposing Eq. (IB. II) into its components in the two 
complementary subspaces a and 6: 

p a (t) = -iL aa p a {t)-iL ab p b (t), (B.15) 

p b {t) = -iL ba p a (t)-iL bb p b (t). (B.16) 

The second equation is formally solved by 

p \t) = e~ iLbbt P b (0) 

nt 

-i / dt'e- iLbb{t - t ' ) L ba p a (t'), (B.17) 
J 0 

which inserted into the first equation yields the integro- 
differential equation 

p a (t) = —iL aa p a (t) 

—L ab [ dt'e- iLbb ( t - t ' ) L ba p a {t’) 

Jo 

-iL ab e~ iLibt p b { 0). (B.18) 

This equation incorporates three terms: The first line 
represents the internal evolution in subspace a, which re¬ 
produces the full time evolution of p a if subspaces a and 6 
are decoupled (L ab = L ba = 0). However, if the coupling 
is nonzero, the second line accounts for virtual transitions 
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from the subspace a into subspace b at time t' < t, fol¬ 
lowed by a period of internal evolution in the subspace b 
(mediated by Liouvillian L bb ) and a final transition back 
into subspace a again at time t. The third line is related 
to the initial slip: it accounts for the contribution to p a (t) 
that stems from the initial projection of the density op¬ 
erator p b (0) on subspace b, combined with a transition 
into subspace a at time t. 

To next derive the effective Liouvillian (IB.121) and ini¬ 
tial condition (IB. 131) from Eq. (IB.181) for the measure¬ 
ment setup under study, we make use of the separation of 
time scales governing the dynamics: This allows us to ne¬ 
glect the time evolution in the quasistationary subspace 
in Eq. (IB. 181) . i.e., we perform a Markov approximation 
with respect to the memory induced by the decaying sub¬ 
space (identifying now a = q and b = d), 

P q (t') « P q (t)- (B.19) 


start. Here these terms must be kept and the solution of 
Eq. (IB.201) can therefore be approximated as 


P q {t) = 


e - iL £ * 


P q ( 0) 

-i /Ve + ^V^V(0) 

J 0 

+ o(£). 


(B.22) 


Again, one can exploit that the exponentials in the sec¬ 
ond line of Eq. (IB.221) vary on a different time scale: 
While the factor e~ lL 4 is nonzero only on a short time 
scale ~ 1/g, the factor e +li eff* changes on a much longer 
time scale max(w, l 2 /g), and we may therefore expand 
e+^ltP' « 1 + 0{u/g,l 2 /g 2 ) in Eq. (lR22j) . In the long¬ 
time limit t —> oo, we can then approximate the integral 
well by setting its upper bound t —> oo, resulting in 


Inserting Eq. (IB.191) into Eq. (IB.181) yields 


p q (t) = -i 


p —iL t 

L q % +L <ld^_^ 

eff Ijdd. 


-L dq 


P q (t) 


-iL qd e~ iLddt p d (<S), 


(B.20) 


with the effective Liouvillian 

L q e q s = L qq -L qd ^L qd + 0 (#,£). (B.21) 

The higher-order terms in Eq. (IB.21[) are due to non- 
Markovian corrections A p q {t') = p q (t') — p q (t) to 
Eq. (IB. 191) . They can be estimated as follows: Since the 
exponential in the integral Eq. (IB.181) decays on a time 
scale 1/g, we only need to account for corrections A p q (t') 
for times t' satisfying t—t' < 1/g. Integrating Eq. (IB.181) . 
this yields corrections A p q ~ uj/g,l/g. Multiplied with 
the order l 2 /g of the corrections from the decaying sub¬ 
space, we obtain the estimate of the higher-order terms 
in Eq. (IB. 211) in accordance with the result obtained in 
frequency space. 

The time-local, Eq. (IB.201) is not yet fully Markovian in 
the sense that it still contains explicitly time-dependent 
terms in the time-evolution generator, which there fore 
becomes frequency dependent in Laplace space Il22j . In 
the long-time limit t 1/g this time dependence drops 
out: As we now argue, one can omit the second term 
in the bracket in Eq. (IB.201) while the second line of 
Eq. (IB. 201) must not be neglected. One may drop the first 
term since it gives a correction to the derivative ~ l 2 /g 
on a time scale 1/g, i.e., they result in an accumulated 
correction A p q [t) on the order of ~ l 2 /g 2 , which can be 
neglected. This is usually achieved in standard deriva¬ 
tions of master equations by setting t —> oo in the inte¬ 
gral in Eq. (IB. 181) [48]. By contrast, the corrections from 
the second line of Eq. (IB.201) are of lower order l/g as in¬ 
tegrating Eq. (IB. 201) shows. In many cases, these terms 
do not appear as one often assumes /O d (0) = 0 from the 


p q it) = e-^-VeV 0) + o(#,f), (B.23) 

with the slipped initial state 

PeV 0) = p q (0)~L qd ^ n p d (0). (B.24) 

We have thus arrived again at Eq. (IB.12[1 and (IB. 131) . re¬ 
spectively, particularly emphasizing that Eq. (IB.231) and 
()B.24I) describe the time evolution correctly as long as 
1/g -C t -C g 2 /[max(Z,w)] 3 ; otherwise, the correction 
terms to the effective Liouvillian (IB.211) can accumulate 
to a large error in Eq. (IB.231) . This is also borne out by 
numerical checks that we performed. 


3. Effective Liouvillian for indirect detection 

We next apply the above Liouville-space projection 
technique to the indirect detection setup of the SQD- 
qubit system studied in the Sec. Mil and HVl of the main 
text. In Appendix IB 3 al we first make the connection 
between the projections just discussed and the dynami¬ 
cal variables considered in the main text. After this, we 
provide in Appendix IB 3 bl some important intermediate 
steps in the derivation of the effective Liouvillian (1801) . 
We comment in Appendix IB 3 cl on the Markov approxi¬ 
mation with respect to the memory induced by the sensor 
QD on the qubit. Finally, we discuss the validity of our 
effective Liouvillian in view of the perturbative expansion 
of the kinetic equations in Appendix IB 41 and comment 
on the U = 0 limit of our problem. 


a. Definition of projections and representing matrices of 
the Liouvillian 

To make a connection between the projection p q = P q p 
of the SQD-qubit density operator on the quasistation¬ 
ary subspace and the dynamical variables introduced in 
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Sec. IIII C 21 we exploit the eigenbasis of the Liouvillian 
L 0 [Eq. (1371) ] with left and right eigenvectors V' p and V p , 
respectively. In the following, “=” denotes that we rep¬ 
resent a basis-independent object on the left hand side 
by its components on the right hand side with respect 
to the eigenbasis of Lq. Representing, for example, the 
density operator p in this basis, we obtain 


P = 



T® T 1 

' rv 1 rv 


pltP° 

V plt T ° 


■ PstP 1 

■ Pst T l 


(B.25) 


with X 9 and X d given by Eq. (l55l) in the main text. In 
contrast to the representation (l56l) in Sec. IIII C 21 work¬ 
ing in the eigenbasis of La also fixes a particular ba¬ 
sis for the isospin part, namely the polarization basis 
which we order as e_ = (ei + ie 2 ) /\/2, eo = 17/fi, and 
e + = (ei — ie 2 ) /y/2. The latter vectors are constructed 
from the right-handed orthonormal basis eo, ei = A/A, 
and e 2 = eo x ei = $7 x A/I2A. This yields the three com¬ 
ponents r” = e]j • r n (a = — ,0, +) in Eq. (IB.25|) that 
make up the isospins r" = J2 a T a e a- The projection of 
p on the quasistationary and the decaying subspace are 
next represented as 


p q = P q p = ( ) , (B.26) 

P d = P d P = ( x d ) ■ (B - 27) 

The matrix (Lg k ) pp > = V p t • L 0 ■ V*, of the unperturbed 
Liouvillian [see Eq. (|41jl ] expressed in terms of its left and 
right eigenvectors V p and V p . respectively, is specified 
completely by the two diagonal blocks 

~ iL ” = ( 2 JL)’ < B ' 28 > 

-< d = <R29> 

whence the off-diagonal blocks vanish: L/ = L^ q = 0. 
Here and below we use the shorthand notation M aoi ' for 
a 3 x 3 matrix M in the above-mentioned polarization 
basis. In this polarization basis the cross product Jlx is 
diagonal: 


= Eaa' [ e a ' X e a')] e„ej,, 

/ -in 0 0 \ (B.30) 

= 0 0 0 , 

y 0 0 +in J 


Next expressing the four blocks of the perturbation A = 
L — Lq in this basis, one finds 


A qq 

A qd 

A dq 

j^dd 


0 0 

0 \p st s aa > 

0 0 

0 \s aa ' 


0 

(1 +Pst)^C c 

0 Ac*, 

Ac a (p° t - 3 k/2)A s, 


(1 + Pst) Ac*, 
'T’XSaa' 


(B.31) 

(B.32) 

(B.33) 

(B.34) 


The components of the matrix s aa ' are proportional to 
the (signed) volume of the parallel epiped spanned by the 
vectors e Q , e a >, and A = A/A, 


Sac*' 
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0 +1 0 \ 
+1 0 -1 , 
0 - 10 / 


(B. 35 ) 


and the remaining factor c a is given by the projected 
isospin-to-charge conversion vector 



(B.36) 


with c given by Eq. (l34l) . Finally, Eq. (IB.331) incorporates 
the transition factor 


r = pltPst ~ « (jbst - P° s t) > ( B - 37 ) 

whose dependence on the SQD parameters we thoroughly 
discussed in Sec. nr C 31 in the main text. 

Equation (IB.311) together with Eq. (IB.351) reveals that 
the “direct” perturbation of the quasistationary modes 
for A ^ 0 is not diagonal in the unperturbed eigenbasis. 
This expresses the fact that the mean-field f2 = r2+Pg t A 
is tilted with respect to $7. Moreover, since the Liou¬ 
villian is not Hermitian, L ^ L\ the transition matrices 
are also not simply related by Hermitian conjugation, i.e., 
A dq ^ (A qd y. This means that transitions between qua¬ 
sistationary and decay modes are not always possible in 
both directions (see Fig. a, in contrast to Hamiltonian 
dynamics. In particular, all transitions into the station¬ 
ary charge modes ( q , c) are forbidden (the first rows of 
A qq and A qd are zero). This is a consequence of the sum 
rules guaranteeing probability and isospin conservation 
by tunneling (see Ref. Hi). Physically, the eigenvalue 
of the stationary charge mode must stay pinned to zero. 

We further find from the above equations [see Eqs. 
(IbTToI) - (IbTTTD 1 : 


\\L bb \-\L aa \\ 

- 7 , 

(B.38) 

max(|L afe |, \L ba \) r 

- A, 

(B.39) 

\L qq \ 

- A - 

- A, 17. (B.40) 
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Thus, the time scales on the quasistationary subspace 
and the decaying subspace separate as required by 
Eq. (IB. Ill) and satisfied by our assumption A <C 7. 
Moreover, the weak-coupling assumption (IB.101) is ful¬ 
filled in the weak-measurement limit A <C 7. For further 
discussion of the consistency see also Appendix IB 41 

b. Computation of effective Liouvillian 

Exploiting Eqs. (IB.26I) and (IB.141) . the effective evolu¬ 
tion of the quasistationary modes can be expressed as 

Uw) = - iL £ x '< f > = L ’*U\ t) ) 

(B.41) 

Consistently accounting for terms up to second order in 
A in the high-temperature limit, the effective Liouvillian 
reads 

L q e | = L^ q + A qq — iA qd ^-A qd + O ■ (B.42) 
Inserting Eqs. (IB. 281) (IB.34[) . we obtain 

- iL % ^ (£ < b -«> 

with 

U = o x E) (b.44) 

\ 2 r 

2 Ol q , > ^PstSotat' (s * s)aa' 

7 

TO (^) . (B.45) 

The injection term I a of Eq. (IB.44[) is thus negli gible 
in the high-temperature limit considered here 123]. In 
Eq. (IB. 451) the first two terms are responsible for the 
qubit precession with frequency ~ A, while the third 
term induces the isospin decay with the a rate ~ A 2 /7, 
which coincides with Eq. (1621) in the main text. 

Inserting Eq. (IB. 431) back into Eq. (IB.411) . we obtain 
for the total isospin evolution: 

7"(t) ^ ] fqBo; ^ ( ^-b e ff )acK / T Q ./(t) 

a ot 

= - iL eS r(t ) (B.46) 

with 

Teff — (I J eff)aa'®a® a ' 

= ^ 7a) e eff,Q6 c g' ct ) (B.47) 

a=0,± 

with Cl a , 7 a , and e e fj a , thus establishing Eqs. 

(1891) in the main part of the paper. 


Finally, the slipped initial state following from 

Eq. (IB.241) reads 

= ( WO)) <B - 48 > 

= X«(0)+A 9d ixf ff (0) + 0(^) (B.49) 

( (T° + 'T 1 )a(0) - jSoca 1 (Pst^ ~ Pst T °)a(0) ) 

+ o(^) (B.50) 

with t™ = el ■ r n . Computing r cff (0) = J2 a T"eff,a(0) e a 
from Eq. (IB.501) . we arrive at Eq. (1511) from the main 
text. The set of states with zero slip is discussed in de¬ 
tail in Appendix [C] and it is shown that the “size” of 
this set grows with the “distance” of the sensor from the 
stationary state. 

This completes the derivation of the effective isospin 
evolution discussed in Sec. IIV Al 


c. Markov approximation relative to the sensor QD 

Regarding the discussion of non-Markovian effects, we 
have to distinguish between two different memory ef¬ 
fects, namely those arising from the electrodes (imposed 
on the qubit-SQD system) and those arising from the 
sensor quantum dot (imposed on the qubit only). As 
explained in Appendix IA 21 we account for the leading 
non-Markovian effect on the composite qubit-sensor QD 
system from the electrodes , which is induced by the tun¬ 
nel coupling; these effects are contained in the effective 
Liouvillian L in the kinetic equations (1261) . which forms 
the starting point of the above analysis. Their effect is 
to modify the coupling of the quasistationary and decay 
modes through the coherent backaction terms in the ki¬ 
netic equation (1551) as explored in Appendix [A2] and Sec. 
IIII B 41 

By contrast, the frequency dependence of the effective 
Liouvillian L e g for the reduced qubit system, Eq. (IB. 71) . 
incorporates non-Markovian effects on the qubit due to 
memory of the sensor quantum dot after the electrodes 
were integrated out (thus effectively of SQD plus elec¬ 
trodes) . These are accounted for by our general Laplace- 
space approach given in Appendix IB II in an exact way 
relative to Eq. (1261) . In principle, these non-Markovian 
effects can be already studied based on the solutions of 
our full kinetic equations (1261) . However, the expressions 
CD and (1721) provide a more convenient starting point 
to gain further insight on an analytical level. 

In the main text, we focused for further illustration on 
the high-temperature limit that allows us to make the 
Markov approximation (IB.191) when expanding the effec¬ 
tive Liouvillian (I8U1) to lowest order in A/7. In this limit, 
SQD-induced memory effects on the qubit are thus ne¬ 
glected. Moreover, we neglect 0(T 2 /T) corrections to the 
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SQD tunneling rates and, thus, non-Markovian correc¬ 
tions from the electrodes are also consistently neglected. 
This is reasonable because if memory effects of the SQD 
are not accounted for, then memory effects from the elec¬ 
trodes should have no effect a fortiori. 

4. Validity of perturbation expansion 

In this part of the Appendix, we collect various re¬ 
marks on the validity and the limitations of our pertur¬ 
bation theory. 

a. Effective Liouvillian [M 

As stated in Sec. IIIBI our kinetic equations are appli¬ 
cable as long as T /T -C A/T -C 1. One may now wonder 
whether it is permissible to expand the denominator in 
Eq. (1691) in orders of A /T and to truncate this expansion. 
In general, the answer is no, for the following reason: To 
lowest order in A/T, the effective Liouvillian scales as 
A 2 /r. The coherent backaction and and the cotunneling 
terms then yield corrections of order A 2 /r x T/T; how¬ 
ever, higher-order corrections in A/T are at least of order 
A 2 /r x A/T and are therefore not neg ligible once we in¬ 
clude terms of order A 2 /T x T/T |l24l |. Therefore, even 
though we start from the weak-measurement limit, one 
must not expand the denominator in Eq. (ESD from the 
start. 

However, if we first neglect all T/T corrections (dis¬ 
sipative and coherent backaction and T 2 /T corrections 
to the SQD rates), then one can consistently expand the 
effective Liouvillian in A/T. The lowest-order approxi¬ 
mation to this is investigated in Sec. HYI and requires high 
temperatures so that the T/T corrections are sufficiently 
small. 


b. Kinetic equation (M3) 

The reader may also wonder whether one should not 
include terms of order A 2 into the kinetic equation (1551) 
since the lowest order to the effective Liouvillian scales 
at least quadratically in A. However, such terms appear 
only in combination with additional tunneling processes 
as the internal interaction L ~ A is treated without ap¬ 
proximation in the kinetic equations. Thus, terms in¬ 
cluding A 2 must be at least of order TA 2 /T 2 and are 
therefore strongly suppressed in the high-temperature 
limit that we consider in this paper. For example, if 
these terms appear in A qq , they lead to corrections of or¬ 
der TA 2 /T 2 -C A 2 /r x T/T since T/T -C 1. They are 
even less important if they appear in A qd , A dq , A qq , where 
they would lead to terms of order A 2 /T{AT/T 2 ) when 
inserted into Eq. (R7J1) . 

Moreover, one may also wonder whether terms of or¬ 
der r 3 /T 2 should be included into the kinetic equation 


(l26l) since they can also modify the stochastic backac¬ 
tion. However, in the weak-coupling, weak-measurement 
limit considered in this paper, those result in corrections 
of even higher-order A 2 /TfT 2 /T 2 ) -C A 2 /T. Yet, deep 
in the Coulomb blockade regime (which is beyond the 
present scope), these terms should play an important role 
since the transition factor r [Eq. (1621) ] - which controls 
the strength of the total backaction together with c — is 
exponentially suppressed by virtue of the effective cancel¬ 
lation of T 2 /T and TX/T terms explained in Sec. IIII 0~3l 
We expect higher-order corrections to cause deviations 
from this exponential suppression of the transition fac¬ 
tor r, resulting in an algebraic scaling ~ l/(e — /J, r ) n 
but with an exponent of higher than that of cotunneling 
broadening, n > 1, in agreement with other theoretical 
work (see Sec. IV Cl ) . We emphasize that the main point 
of our work is just to explain physically the cancellation 
of the expected leading power law n = 1. A much more 
elaborate analysis is needed to find the actual power law 
(including the complications of a nonstationary detector 
with strong local Coulomb interaction, etc.). 

c. Exact result for (7 = 0 

Finally, to further support the cancellation in Eq. (1551) 
that we found here for large Coulomb interaction U and 
perturbatively in F, we computed in a separate calcula¬ 
tion the effective Liouvillian for a noninteracting sensor 
QD (U = 0) nonperturbatively in T but only to lead¬ 
ing order in A. This releases the condition T/T -C A/T 
and can be used to investigate also the opposite regime 
A/r < F/T(-C 1). These results actually confirm the 
mitigation of the cotunneling-induced stochastic backac¬ 
tion by the coherent backaction also in this case. Yet, it 
remains an interesting future question to understand the 
role of the coherent backaction in the low-temperature, 
strong qubit-sensor coupling regime when also the strong 
interaction U is accounted for. 

Appendix C: Initial states without slip 

In this Appendix, we characterize the set of initial 
qubit-sensor QD density operators with zero initial slip 
[see Eqs. (17X1) and (1511) ] in the reduced dynamics of the 
qubit in the high-temperature limit. We show that this 
imposes a strong condition: Initial states with zero slip 
form a set of measure zero, thus making the initial slip a 
relevant source of errors in indirect detection unless the 
qubit-SQD quantum state ( not just the qubit state) is 
under accurate control. In addition, increasing the non- 
stationarity of the sensor reduces these sets of zero slip. 
We also formulate the constraint in terms of relations be¬ 
tween the reduced qubit state and the composite qubit- 
SQD state and discuss how it relates to the factorizability 
of the qubit-SQD state. 
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FIG. 13: Geometric restrictions on the qubit-sensor density 
operator imposed by requiring zero initial slip. We show a 
2D cross section of the 3D construction described in the text, 
i.e., the circles represent spheres and the dotted lines indicate 
the boundaries of a cylinder around the line xX. 

Kinematic restrictions. It is most convenient for the 
following considerations to work with the representation 
of the initial qubit-SQD state, 

p(0) = ^P"®i(p"(0)i Q +r"(0)-f), (C.l) 

n 

in terms of the occupation probabilities p n ( 0) and the 
charge-specific isospins r n (0). We recall that for this 
state to be a valid quantum state, p n ( 0) and r”(0) have 
be real (for pi = p), )T) n p n (0) = 1 (for trp = 1), and the 
magnitudes 

|r"(0)| < p n { 0) (C.2) 

have to be restricted (for positivity, p > 0). We note that 
the reduced state of the SQD and the qubit, 

p s := trp(0) = 5> re ( 0)P", (C.3) 

n 

pq ■= trp(o) = §(iQ+T(o)-f), ( c - 4 ) 

are completely specified by the occupation probabilities 
p n { 0) (in fact by only one) and the Bloch vector t(0) = 
t°(0) + t 1 (0), respectively. When p n ( 0) = p" t , the SQD 
is stationary, but this does not imply that the qubit-SQD 
state state p(0) is factorisable (see below). 

Charge-specific isospins. We first investigate the slip 
based on the charge-specific isospins T n (0) since these 
have simple kinematic constraints, allowing for an easy, 
complete characterization. The zero-slip condition ob¬ 
tained from Eq. ©, 

0 = A x (r° + r 1 ) , (C.5) 

is expressed using rescaled isospin vectors r° = Pg t r°(0) 
and r 1 = —Pg t r 1 (0). For an arbitrary initial state, the 
two three-dimensional vectors r° and r 1 are taken from a 
six-dimensional set that is constrained only by the pos¬ 
itivity of p(0): the vectors r° and r 1 have to lie within 
spheres of different radii given by |r°| ^ R° := p\ t p°{ 0) 


and | r 1 1 ^ R 1 := p° t p 1 (0), respectively. The radii are 
equal if and only if the reduced sensor state p$ is sta¬ 
tionary: the condition Pg t p 1 (0) = Pg t p°(0) is equivalent 
to p n ( 0) = p" t (n = 0,1) due to the normalization condi¬ 
tion = J2nPst = !• 

From this set of valid initial states let us now construct 
those which have zero slip. According to Eq. (1C.51) . this 
requires the sum of the vectors r° and r 1 to lie on the line 
defined by the measurement vector: r° + r 1 = xX with 
any x £ R. For the construction, sketched in Fig. [T3l 
first draw a sphere with radius R° (blue) and draw the 
line a; A through its origin (black). For each vector r° in 
this sphere draw a second sphere of radius R 1 centered 
at its tip (red). The set of vectors r 1 that give zero- 
slip state are just given by the intersection of this second 
sphere with the line xX. This is a set of measure zero. 
Moreover, from the figure it is clear that the construction 
is possible only if r° is inside a cylinder of radius R 1 with 
the line xX as its axis. Whenever r° lies outside this 
cylinder, it is not possible at all to construct a zero-slip 
initial state. 

The radii R n are controlled by the initial reduced sen¬ 
sor state through p n (0) and there are two extreme limits: 

(i) Stationary state p n ( 0) = p" t : R° = R 1 . In this 
case for every r° one can find an r 1 giving a zero-slip 
state. Still, the subset has measure zero in the total set 
of possible states. 

(ii) Integral charge state p o,1 (0) = 0,p 1,o (0) = 1: 
R 0 ’ 1 = 0, R 1,0 = pl t . In these two cases, the charge- 
specific isospin t o,1 (0) = t(0) coincides with the total 
isospin, while r 1,0 (0) = 0. The slip is r c ff(0) — t(0) = 
=FA x Ps t °r( 0)/7 and can be avoided only if r°’ 1 (0) is 
collinear to A so that the tip of r 0,1 lies on the line xX 
to allow for r 1,0 = 0 as sketched in Fig. [Tdl a'). 

As the sensor deviates from the stationary state, the 
radii R 1 and R° differ and the above possibilities for con¬ 
structing initial zero-slip states are reduced. Since in any 
case the zero-slip states are sets of measure zero, it is clear 
that most preparation errors of the sensor QD will lead 
to an initial slip on the qubit , i.e., the backaction-induced 
initial slip generates additional errors beyond the control 
over the qubit. 

Restriction on the qubit isospin. It is instructive to 
further clarify the restrictions that the above imposes 
on the reduced qubit state, a density operator completely 
characterized by r(0), relative to the composite qubit¬ 
sensor state p(0). To this end we change the variables 
to 


r(0) 

= T°(0) + T 1 (0) 

(C.6) 

M0) 

= pl t p 0 (0)-p 0 st p\ 0) 

(C.7) 

m 

= p 1 (0)x°(0) — p°(0)x 1 (0) 

(C.8) 


The conditions imposed by the positivity p(0) on r( 0) 
and <5(0) are not easy to formulate and will be ignored 
in the following (note that |t(0)| ^ 1 is only necessary, 
not sufficient). Rewriting the zero-slip condition (1811) by 
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inserting r°’ 1 (0) = p°’ 1 (0)t(0) ± 5(0) yields 

0 = Ax [5(0) + 5 st (0)r(0)] (C.9) 

Noting that 5(0) =0 corresponds to a factorizable state 
[see Eq. and 5 st (0) = 0 to a stationary reduced 

sensor state (see above), we have four cases in which there 
is zero initial slip: 

(i) nonfactorizable, nonstationary initial state 5(0) ^ 
0, 5 s t(0) ^ 0: The qubit state and the qubit-SQD correla¬ 
tions must be fine tuned such that 5(0) +5 st (0)r(0) = xX 
or some x € R. 

(ii) nonfactorizable, stationary initial state 5(0) ^ 0 , 
5 s t(0) = 0: The qubit state can be arbitrary, but the 
qubit-SQD state must have have very special correlations 
such that 5(0) oc A. 

(iii) factorizable, nonstationary initial state 5(0) = 0, 
5 s t(0) ^ 0: The qubit state must be prepared in a 
measurement-basis state t(0) oc A. 

(iv) factorizable, stationary initial state 5(0) = 0, 
5 s t(0) = 0: no conditions, there is no slip. 

Equation (1C.91) emphasizes that the set of states with¬ 
out slip does not simply coincide with factorizable initial 
states (“uncorrelated states”). It further shows that the 
more the sensor deviates from the stationary state, now 
quantified by 5 s t(0), the larger the initial slip becomes. 
This concludes our discussion of the states with zero slip 
and, as the previous analysis showed, the zero-slip states 
are subsets of measure zero of the set of possible initial 
states. The experimental possibilities to avoid the slip 
[case (iv)] are discussed in Sec. IIIIEl and Sec. IVII 


Appendix D: Renormalization-induced qubit phase 
kicks 

In this final Appendix, we explain how the coherent 
backaction appears in the visibility (11061) as an addi¬ 
tional contribution that cannot be understood in a semi- 
classical stochastic picture discussed in Sec. IV Al We 
show that, loosely speaking, the coherent backaction in¬ 
duces additional “phase kicks” on the qubit that par¬ 
tially “undo” the phase kicks induced by the stochastic 
backaction that result in decoherence. As a result, the 
decoherence is mitigated. We emphasize from the start 
that these phase kicks are completely unrelated to the 
initial slip: They do not lead to a phase shift of coherent 
isospin precession but instead only affect the net qubit 
decoherence. 

Our objective in the following is to merely further il¬ 
lustrate the physical origin of the coherent backaction by 
calculations that are certainly not rigorous. The proper 
treatment is achieved by our kinetic equation (1261) which 
is based on the systematic real-time diagrammatic ap¬ 
proach. Still, we believe the following may be instructive. 

We start from the expression (11061) for the visibility, 

D{t) = tr[(|+)(—| <g> l Si j) 

xe-««(|+><—| ® PS R)e im ], (D.l) 


assuming the initial qubit-environment state factorizes 
and H is given by Eq. ©■ We furthermore neglect the 
spin degree of freedom for the SQD and thereby also 
the Coulomb interaction effect, i.e., the SQD takes two 
charge states |0) and |1). As noted in Sec. IV Al and 
Appendix IB 41 neither is essential for the coherent back- 
action. We next split up the Hamiltonian H = Ho + Ht 
and expand Eq. (ID. II) in the tunneling Ht of electrons 
between the SQD and the reservoir and vice versa. In 
O(T) ~ O(Ht), the following term contributes (there 
are more terms; we just consider a relevant one): 


D(t) ~ f dhdt 2 tr [{+\e- lH °( t - tl) H T e- lHotl \+) 

J S,R 

0<ti 

® PSR(-\e iHo{t - t 2 ) H T e iH ^\-)} (D.2) 


Importantly, the tunneling process for the ket-evolution 
(+| <r- (+| may happen at a time t\ different from the 
time t 2 for the tunneling process of the bra evolution 
(—| -A |—). The coherent evolution in between is re¬ 
sponsible for an additional phase shift. To see this, 
we assume the SQD to be initially unoccupied, i.e., 
PSR = |0)(0| ® PR.fi, and to end up in a singly occu¬ 
pied state |1)(1|. We assume again for simplicity J7 = 0. 
Evaluating the trace over the electrodes, we get 


m 


J dtidt 2 ^2, J du)f+(u)r r (u)e tl ' > 


X 


e -i(e+^-oi)(ti-t 2 ) _|_ e -i(-e+!+w)(ti-t 2 ) 


(D.3) 


where the two terms in the bracket relate to the two 
cases, t± > t 2 , and t\ < t 2 in Eq. (ID. 21) . respectively. It 
is easy to show that the terms in the second line arise 
from SQD-qubit coherences of the type |l,+)(0, — | and 
|0, V)(l, — |, respectively. Finally, the SQD reaches the 
singly occupied real state |1)(1|, giving rise to the phase 
shift g-i’Ht-ti)' The latter is related to the stochastic 
backaction (since the time t\ is random), while the for¬ 
mer phase shift contains the coherent backaction as we 
explain next. For this purpose, we redefine r := ti — t 2l 
assume large t, and introduce a small imaginary part into 
the exponentials to ensure convergence. Since r can take 
nearly all positive values for large f, the coherent phase 
shift becomes approximately 


dr 


= H) 


g— i(e+-|— ui— i0)r . i( — e+^+u— i0)r 

1 1 


s V A/2 — to — 10 s — A/2 — to iO 


1 d 1 

—2 Im-— — i\— Re ■ 


e — to + i0 de £ — to + i 0 ’ 


(D.4) 


where the last step holds in first order in A. We ignore 
the first part, which is independent of A, and insert the 
second part into Eq. (ID.3[) . which yields a term propor¬ 
tional to the renormalization function. When taking the 
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time derivative of Eq. (ID.31) , it reproduces the coherent 
backaction: 

D ~ ( £ = ~ iK ( D ' 5 ) 

which coincides with the effect of the term r 1 = n\ x r° 
in the kinetic equation (1261) . 

What our heuristic argument clarifies is why a semiclas- 
sical approach as discussed in Sec. IV Al is not capable of 
reproducing the coherent backaction: this approach cru¬ 
cially relies on the assumption that the charge state 
of the SQD is a classical variable , i.e. n(t) is just 
a fixed (but random) function of t. Thus, the transi¬ 


tion 0—^1 happens on both the bra- and the ket-part 
of the evolution at the same time t\ = t 2 (while the 
time t\ is random). In this case, the coherent phase 
shift just vanishes. This amounts formally to replacing 
1/x —> —inS(x) and “by hand” dropping the principal 
value integral term P(l/x) that contains renormaliza¬ 
tion effects. In the single-step Born-Markov approach 
discussed in Sec. IV Bl the problem is to find a proper 
description of the environmental state that contains the 
intermediate coherences |1)(0| and |0)(1|. Our above non¬ 
rigorous derivation of the coherent backaction terms thus 
illustrates that virtual fluctuations in nanoscale sensors 
can have a real impact on the measurement backaction. 
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[99] The first step of our approach is in principle generally 
possible at least for weak tunnel coupling T, i.e., not 
limited to the weak-measurement situation we study in 
this paper. 

[100] Here we use Re ip( 1/2 + ix/2x) ~ In |a:| for |*| 1. 

[101] Note that we investigate the impact of the coherent 
backaction by setting here and below “by hand” k = 
0 in our results. Although one can express the non- 
Markovian correction factor as g = 1 + k, we do not 
set r) = 1 since this would affect the stochastic backac¬ 
tion and would not lead to the comparison we intend to 
make. 

[102] Note that the dynamics of the quasistationary modes 
does not contain information about the response of the 
sensor QD to the qubit since the coefficient of the sta¬ 
tionary charge mode V/ trivially equals 1. To describe 
the sensor response, one must compute the dynamics 
of <5 s t = PstP° — PstP 1 : see Eq. (1551) (while the isospin 
degrees of freedom may be projected out for that pur¬ 
pose). 

[103] The coherent backaction does not only affect the transi¬ 
tions into the decaying subspace but also affects the evo¬ 
lution in the decay space by modifying A dd in Eq. (1561) . 
However, these terms yields only a small correction as 
compared to the large free evolution term \LQ q \ ~ 7 . 

[104] It is known that at finite temperature the condition for 
a resonance is not e = p r , but there is an offset linear 
in T with a coefficient that grows with junction asym¬ 
metry. This is related to Coulomb interaction effects on 
the QD and appears already in O(T) [76j. This causes 
the maximum of the decoherence rates to lie at nonzero 
V g (see Fig. [5]) and also causes the shift with changing 
junction asymmetry g in Fig. [G] 

[105] The kinetic equations for sensor QD plus qubit are con¬ 
sidered here in the wide-band limit with respect to the 
bandwidth W of the electrodes to which the sensor is 
coupled. 

[106] The transition factor satisfies r > 0 even when cotun¬ 
neling and coherent backaction terms are included, see 
Eq. (1651) . 

[107] Note that the decay rates would diverge in the limit 
A/r —> 00 . 

[108] We stress that Ref. [rf] is based on a kinetic equation to 
leading order in F, i.e., cotunneling and renormalization 
effects are excluded from the start. Thus, the authors 
do not address the same questions as we do here. 

[109] Note that T e ff(0) differs from the initial isospin r(0) due 
to the initial slip m and therefore F is not necessarily 
limited to the range [— 1 , + 1 ]. 

[110] We note that applying Eq. (11041) requires a calculation 
of the charge-specific isospins from the full kinetic equa¬ 
tions (|26|. Using only the projections quasistationary 
subspace, X 9 (f) = (1 r(t))^, see Eq. (IB. 261) . would just 
yield zero when inserted into Eq. (fT04l) . 

[111] We introduce a factor 1/2 for convenience to adjust the 
amplitude of the fluctuations to A£ = 2 , as usually used 
for a random telegraph process. 
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We adjusted the notation in Ref. [ 66 | to our notation by 
replacing 1/to = 1 /t u + 1 /n = 7 and Jt — > A/ 2 . 

The reason is that the dissipative backaction does not 
affect the coupling between the quasistationary and de¬ 
cay modes [see Eqs. (l 6 Tll (l63ll l. 

In Redfield’s notation in Ref. | 6 jfl, the random pertur¬ 
bation is for our model given by G(t) = |£(f)A ■ r. 

To compare with Eq. (7) of Ref. [25|, we note that for 
our model there is no direct hybridization Aj = 0 be¬ 
tween the orbital levels of the SQD, i.e., cos ( 6 j) = 1, 
and sin ( 9 j) = 0. Furthermore, the relaxation rate Fry 
to the stationary state is identified with 7 here. 

This is a consequence of the approximation K 2 ,j( A) « 
K 2J (0) employed below Eq. (3) in Ref. [*22^. 

We express the parameters g and e in Ref. [23|] in terms 
of our notation as g = A/T and i = 2e/T. 

The correlation time is set by the time dependence of the 
contraction functions in the diagrammatic expansion 
(reservoir correlation functions), which drop on the time 
scale of the inverse temperature 1 /T, see Eqs. (93) and 
(100) in Ref. Note that the contribution from the 
stationary state p s t just vanishes, i.e., W(t — t')p s t = 0 
for any time t'. 

This can be understood as follows: the finite fre¬ 
quency kernel is given by Eq. (A 6 ) in Ref. [32| 
by replacing iO —>• z/T. Each frequency deriva¬ 
tive in d n W acts on a propagator, i.e., it involves 
expressions such as 
1 0 1 


a 


dz z/T+Xi-^LQs-^/T 

t d(L QS -vi)/T z/T+x i -(L QS - l j.i)/T > which can 

be rewritten as a derivative of dimensionless 
energy ratios and a factor 1 JT for each z- 
derivative. This yields schematically bd n W ~ 


(n 


Ei' 

T , 


d n 


(Lq S ~^\ 

v T )’ 


9[(Lqs-G/ t ] 
dimensionless function that contains 
integrals. Thus, d n W ~ V/T n since 
higher-order corrections). 

The first-order kernel reads 


where I is a 
the frequency 
r (plus 


w 


W r (z) 
d/dL f dxf 


f dxf ± {x)/[z + x 
{x)/[z + x — L\. 


~ L], 


schematically 
9 1 W r (2) ~ 


To obtain a solution p(t) that is also consistently ex¬ 
panded to all orders of A and T, one should insert the 
expanded density operator p(t) = p 0,0 (t) + p F,0 (t) + 
p°’ A + ... together with the expanded kernels into the 
kinetic equation p(t) = —iLp(t), and compare left- and 
rightOhand side in all orders. 

See also the discussion of microscopic derivations of 
master equations in Ref. [i3 |. 

We checked that when accounting for the leading-order 
expression in Eq. (IB.441) for lower temperatures, one 
reproduces the stationary state of Ref. [33] with nonzero 
stationary isospin r s t. 

We note that we still consider the weak-measurement 
regime (A/T <S 1). The validity of our kinetic equations 
is just limited to high temperatures implying T /T <C 
A/T. 










